снартев 30
ANISOTROPIC AND INHOMOGENEOUS COSMOLOGIES
§30.1. WHY IS THE UNIVERSE SO HOMOGENEOUS AND ISOTROPIC?
This chapter is entirely Track 2.
The main text requires no special preparation, although Chapters 27-29 would be helpful.
Box 30.1 contains more technical sections: ideal preparation for it would be Chapters 4,9-14, 21, and 27-29, plus §25.2; minimal preparation would be exercises 9.13,9.149.13,9.14, and 25.2, Chapter 21 through §21.8, and §§ 27.8,27.1127.8,27.11, and 29.2.
Chapter 30 is not needed as preparation for any later chapter.
Motivation for studying inhomogeneous and anisotropic cosmologies: Why is universe so uniform?
The last three chapters studied the Friedmann cosmological models and the relatively satisfactory picture they give of the universe and its evolution. This chapter describes less simplified cosmological models, and uses them to begin answering the question, "Why are the very simple Friedmann models satisfactory?" This question is intended to probe more deeply than the first, obvious answer-namely, that the models are satisfactory because they do not contradict observations. Accepting the agreement with observations, we want to understand why the laws of physics should demand (rather than merely permit) a universe that is homogeneous and isotropic to high accuracy on large scales. Because this question cannot be answered definitively in 1972, many readers will prefer to omit this chapter on the first reading and return to it only after they have surveyed the major results in other areas such as black holes (Chapter 33), gravitational waves (Chapters 35-37), and solar-system experiments (Chapter 40).
The approach described here to the question "Why is the universe so highly symmetric?" is to ask Einstein's equations to describe what would have happened if the universe had started out highly irregular.
The first step in this approach is to ask what would have happened if the universe had started a little bit irregular. This problem can be tackled by analyzing small perturbations away from the high symmetry of the Friedmann models. Such an analysis is most fruitful in its discussion of the beginnings of galaxy formation, and
in its ability to relate small upper limits on the present-day anisotropy of the microwave background radiation to limits on density and temperature irregularities that might have existed ten billion years ago, when the radiation was emitted. These studies are described so well in the book by Zel'dovich and Novikov (1974) [see also Field (1973), Peebles (1969), Peebles and Yu (1970), Jones and Peebles (1972), and references cited therein] that we omit them here.
Another approach is to allow large deviations from the symmetry of the Friedmann universes, but to put the asymmetries into only a few degrees of freedom.
§30.2. THE KASNER MODEL FOR AN ANISOTROPIC UNIVERSE
The prototype for cosmological models with great asymmetry in a few degrees of freedom is the Kasner (1921a) metric,
{:(30.1)ds^(2)=-dt^(2)+t^(2p_(1))dx^(2)+t^(2p_(2))dy^(2)+t^(2p_(3))dz^(2):}\begin{equation*}
d s^{2}=-d t^{2}+t^{2 p_{1}} d x^{2}+t^{2 p_{2}} d y^{2}+t^{2 p_{3}} d z^{2} \tag{30.1}
\end{equation*}
which was first studied as a cosmological model by Schücking and Heckmann (1958). In this metric the p_(i)p_{i} are constants satisfying
Each t=t= constant hypersurface of this cosmological model is a flat three-dimensional space. The world lines of constant x,y,zx, y, z are timelike geodesics along which galaxies or other matter, treated as test particles, can be imagined to move. This model represents an expanding universe, since the volume element
sqrt(-g)=sqrt((3)g)=t\sqrt{-g}=\sqrt{(3) g}=t
is constantly increasing. But it is an anisotropically expanding universe. The separation between two standard (constant x,y,zx, y, z ) observers is t^(p_(1))Delta xt^{p_{1}} \Delta x if only their xx-coordinates differ. Thus, distances parallel to the xx-axis expand at one rate, ℓ_(1)propt^(p_(1))\ell_{1} \propto t^{p_{1}}, while those along the yy-axis can expand at a different rate, ℓ_(2)propt^(p_(2))\ell_{2} \propto t^{p_{2}}. Most remarkable perhaps is the fact that along one of the axes distances contract rather than expand. This contraction shows up mathematically in the fact that equations (30.2) require one of the pp 's, say p_(1)p_{1}, to be nonpositive:
As a consequence, in a universe of this sort, if black-body radiation were emitted at one time tt and never subsequently scattered, later observers would see blue shifts near one pair of antipodes on the sky and red shifts in most other directions. In terms of this example, the fundamental cosmological question is why the Friedmann metrics should be a more accurate approximation to the real universe than is this Kasner metric.
Kasner metric: an example of an anisotropic model universe
§30.3. ADIABATIC COOLING OF ANISOTROPY
In seeking an answer, ask a question. Ask, in particular, what would become of a universe that starts out near t=0t=0 with a form described by the Kasner metric of equation (30.1). This metric is an exact solution of the vacuum Einstein equation G=0\boldsymbol{G}=0. It approximates a situation where the matter terms in the Einstein equations are negligible by comparison with typical non-zero components of the Riemann tensor. Schücking and Heckmann (1958) give solutions with matter included as a pressureless fluid. In this situation, the curvature of empty spacetime dominates both the geometry and the expansion rate at early times, t longrightarrow0t \longrightarrow 0; but after some characteristic time t_(m)t_{m} the matter terms become more important, and the metric reduces asymptotically to the homogeneous, isotropic model with k=0k=0.
This example illustrates the possibility that the universe might achieve a measure of isotropy and homogeneity in old age, even if it were born in a highly irregular state. Whether the symmetry of our universe can be explained along these lines is not yet clear in 1972. The model universe just mentioned is only a hint, especially since the critical parameter t_(m)t_{m} can be given any value whatsoever.
The standard Einstein general-relativity physics of this model can be described in other language (Misner, 1968) by ascribing to the anisotropic motions of empty
Anisotropy energy
Adiabatic cooling of anisotropy spacetime an "effective energy density" rho_("aniso ")\rho_{\text {aniso }}, which enters the G_(00)G_{00} component of the Einstein equation on an equal footing with the matter-energy density, and thereby helps to account for the expansion of the universe:
For pressureless matter gamma=1\gamma=1; for a radiation fluid gamma=4//3\gamma=4 / 3; for a nonrelativistic ideal gas gamma=5//3\gamma=5 / 3 ).
This arrangement of the Einstein equation allows one to think of the anisotropy motions as being adiabatically cooled by the expansion of the universe, just as the thermal motions of an ideal gas would be. Since the adiabatic index for homogeneous anisotropy is gamma=2\gamma=2, the anisotropy will be the dominant source of "effective energy" in a highly compressed state, whereas the matter will dominate in an expanded state.
§30.4. VISCOUS DISSIPATION OF ANISTROPY
The model universe sketched above can be further elaborated by introducing dissipative mechanisms that convert anisotropy energy into thermal energy. Suppose that
such an anisotropic universe were filled at one time with thermal radiation. If the radiation were collisionless or nearly so, the quanta moving parallel to the contracting xx-axis would get blueshifted and would develop an energy distribution corresponding to a high temperature. The quanta moving parallel to the other (expanding) axes would be redshifted to an energy distribution corresponding to a low temperature. Any collisions taking place between these two systems of particles would introduce a "thermal contact" between them, and would transfer energy from the hot system to the cold one, with a corresponding large production of entropy. This provides an irreversible dissipative process, which decreases rho_("aniso ")\rho_{\text {aniso }} and increases rho_("radiation ")\rho_{\text {radiation }} relative to the values they would have had under conditions of adiabatic expansion. [For further details, see, e.g., Matzner and Misner (1972).]
It is possible that both the adiabatic cooling of anisotropy and the dissipation of anisotropy by its action on a gas of almost collisionless quanta have played significant roles in the evolution of our universe. In particular, neutrinos above 10^(10)K10^{10} \mathrm{~K} may have undergone sufficient nu\nu-e scattering to have provided strong dissipation during the first few seconds of the life of the universe.
§30.5. PARTICLE CREATION IN AN ANISOTROPIC UNIVERSE
Adiabatic cooling and viscous dissipation might not be the chief destroyers of anisotropy in an expanding universe. More powerful still might be another highly dissipative process, which might occur at still earlier times, very near the initial "singularity." This is a process of particle creation which was first treated by DeWitt (1953), then explored by Parker (1966 and 1969) for isotropic cosmologies and finally by Zel'dovich (1970) in the present context of anisotropic cosmologies. In this process one again turns to the Kasner metric for the simplest example, but now quantummechanical considerations enter the picture. One realizes that not only would real quanta propagating in different directions be subject to red shifts and blue shifts, but that virtual quanta must be considered as well. Vacuum fluctuations (zero-point oscillations) entail a certain minimum number of virtual quanta, which are subject to the redshifting and blueshifting action of the strong gravitational fields. Virtual quanta that are blueshifted sufficiently violently can materialize as real particles, thanks to their energy gain. In this context "sufficiently violently" means not adiabatically.
In an adiabatic expansion, the number of particles does not change, although the energy of each one does. This adiabatic limit is just the geometric-optics approximation to wave equations, which was discussed in §22.5\S 22.5§. There one saw that, if spacetime were not flat on the scale of a wavelength, then the wave equation could not be replaced by a particle description with conserved particle numbers. Thus, the adiabatic limit (geometric-optics approximation) is violated in the conditions of high curvature near the singularity at the beginning of the universe.
By studying wave equations in the Kasner background metric, Zel'dovich and Starobinsky (1971) find quantitatively the consequences of the failure of the adia-
Creation of particles by anisotropy of expansion
Anisotropy might have created the matter content of our universe, damping itself out in the process
Inhomogeneous cosmological models:
(1) with spherical symmetry
(2) with (rather symmetric) gravitational waves
(3) near a singularity, with few or no symmetries
batic approximation near the singularity. Classically, the amplitudes of waves at frequencies comparable to the Hubble constant for any given epoch increase faster than a simple blue-shift calculation would imply (amplification through parametric resonance). Quantum-mechanically, the same amplification, applied to zero-point oscillations, leads to the creation of particle-antiparticle pairs. The calculations indicate that this effect is very strong at the characteristic time t_(q)=sqrt(Gℏ//c^(5))≃10^(-43)t_{q}=\sqrt{G \hbar / c^{5}} \simeq 10^{-43} sec. (All calculations performed thus far are inadequate when the effect becomes strong, thus for t <= t_(q)t \leqq t_{q}.
For the creation of massless particles, it is essential that an anisotropically expanding universe be postulated (except for scalar particles, for which particle creation occurs already in the Friedmann universe, unless the particle satisfies the conformalinvariant wave equation). The isotropic Friedmann universes are all conformally flat, so that solutions of the wave equation for a field of zero rest mass can be given in terms of solutions for flat-space wave equations where there is no particle creation. There is some particle creation even in the isotropic Friedmann universe when the particle has finite rest mass and low energy. However, the particle-creation process normally uses anisotropy energy as the energy supply that it converts into radiation energy.
The pioneering work by Parker and Zel'dovich suggests that one should study in detail cosmological models in which the initial conditions are a singularity, and in which quantum effects near the time t=t_(q)t=t_{q} dissipate all anisotropies and simultaneously give rise to the matter content of the model. This program of research, which is in its infancy, seems to require extrapolating laws of physics down to the very natural looking but preposterously small dimension sqrt(Gℏ//c^(5))≃10^(-43)\sqrt{G \hbar / c^{5}} \simeq 10^{-43} sec, or equivalently sqrt(Gℏ//c^(3))∼10^(-33)cm\sqrt{G \hbar / c^{3}} \sim 10^{-33} \mathrm{~cm}.
§30.6. INHOMOGENEOUS COSMOLOGIES
The model universes considered above were all homogeneous although anisotropic. It is also crucial to study inhomogeneous cosmological models, in which the metric has a nontrivial dependence on the space coordinates. One class of such models is spherically symmetric universes, where the matter density, expansion rate, and all other locally measurable physical quantities have spherical symmetry about some preferred origin. Models of this sort were first considered by Lemaitre (1933a,b), Tolman (1934b), and Datt (1938), and were also treated by Bondi in 1947. These models provide a means for studying density perturbations of large amplitude.
A recent tool is making it possible to study large-amplitude, spatially varying curvature perturbations of other symmetries; this tool is the Gowdy (1971, 1973) metrics. These metrics, which are exact solutions of the Einstein equations, represent closed universes with various topologies ( S^(3),S^(1)xxS^(2),T^(3)S^{3}, S^{1} \times S^{2}, T^{3} ) containing gravitational waves. The wave form in these solutions is essentially arbitrary, but all the waves propagate along a single preferred direction and have a common polarization.
A rather different approach to understanding the behavior of inhomogeneous and anisotropic solutions of the Einstein equations has been developed by Khalatnikov,
Lifshitz, and their colleagues. Rather than truncate the Einstein theory by limiting attention to specialized situations where exact solutions can be obtained, they have sought to study the widest possible class of solutions, but to describe their behavior only in the immediate neighborhood of the singularity. These studies give a greatly enhanced significance to some of the exact solutions, by showing that phenomena found in them are in fact typical of much broader classes of solutions.
Thus, in the first large class of solutions studied [Lifshitz and Khalatnikov (1963)], it was found that near the singularity solutions containing matter showed no features not already found in the vacuum solutions. Furthermore, space derivatives in the Einstein equations became negligible near the singularity in these solutions, with the consequence that a metric of the Kasner form [equation (30.1)] described the local behavior of spacetime near the singularity, but with a different set of p_(i)p_{i} values possible at each point of the singular hypersurface. Subsequently, broadened studies of solutions near a singularity [Belinsky and Khalatnikov (1970)] showed that the mixmaster universe [Misner (1969b); Belinsky, Khalatnikov, and Lifshitz (1970)] is a still better homogeneous prototype for singularity behavior than the Kasner metric.
§30.7. THE MIXMASTER UNIVERSE
The simplest example of a mixmaster universe is described in Box 30.1. It shows how, near the singularity, the Kasner exponents p_(i)p_{i} can become functions of time. The result is most simply described in terms of the Khalatnikov-Lifshitz parameter uu :
As one extrapolates backward in time toward the singularity, one finds that the expansion rates in the three principal directions correspond to those of the Kasner metric of equation (30.1), with p_(i)p_{i} values corresponding to some fixed uu parameter. In these mixmaster models, however, the metric is not independent of the space coordinates (the spacelike hypersurfaces can, for instance, have the same 3 -sphere topology as the closed Friedmann universes).
The Kasner-like behavior at fixed uu can persist through many decades of volume expansion before effects of the spatial derivatives of the metric come into play. The role then played by the space curvature is brief and decisive. The expansion is converted from a type corresponding to a parameter value u=u_(0)u=u_{0} to a type corresponding to the value u=-u_(0)u=-u_{0} (which is equivalent, under a relabeling of the axes, to the value u=u_(0)-1u=u_{0}-1 ). Extrapolating still farther back toward the singularity, one finds a previous period with u=u_(0)-2u=u_{0}-2. Throughout an entire sequence u=u_(0)u=u_{0}, u_(0)-1,u_(0)-2,u_(0)-3,dotsu_{0}-1, u_{0}-2, u_{0}-3, \ldots, with u_(0)≫1u_{0} \gg 1, nearly the entire volume expansion is due to expansion in the 3 -direction, whereas the 1 - and 2 -directions change very little, alternating at each step between expansion and contraction. Sufficiently far in the past, however, such a sequence leads to a value of uu between 0 and 1 . This value
Mixmaster universe:
(1) "anisotropy oscillations" explained in terms of Kasner model
(2) as a prototype for generic behavior near singularities
Are there any other generic types of behavior near singularities?
can be interpreted as the starting point for another, similar sequence, through the transformation u longrightarrow1//uu \longrightarrow 1 / u, which interchanges the names of axes 2 and 3 .
The extrapolation of the universe's evolution back toward the singularity at t=0t=0 therefore shows an extraordinarily complex behavior, in which similar but not precisely identical sequences of behavior are repeated infinitely many times. In terms of a time variable which is approximately log(log t^(-1))\log \left(\log t^{-1}\right), these behaviors are quasiperiodic. In the generic example to which the Khalatnikov-Lifshitz methods lead, one has a metric whose asymptotic behavior near the singularity is at each point of the singular hypersurface described by a mixmaster-type behavior, but with the principal axes of expansion changing their directions as well as their roles (as characterized by the uu parameter) at each step, and with the mixmaster parameters spatially variable. [For more details see Belinsky, Lifshitz, and Khalatnikov (1971), and Ryan (1971, 1972).]
It is not yet (1972) known whether there are important solutions or classes of solutions relevant to the cosmological problem, with asymptotic singularity behavior not described by the Khalatnikov-Lifshitz generic case. The difficulty in reaching a definitive assessment here is that Khalatnikov and Lifshitz use essentially local methods, confined to a single coordinate patch, whereas the desired assessment poses an essentially global question. The global approaches (described in Chapter 34) have not, however, provided any comparable description of the nature of the singularity whose necessity they prove. One attempt to bridge these differences in technique and content is the work by Eardley, Liang, and Sachs (1972).
(continued on page 815)
Box 30.1 THE MIXMASTER COSMOLOGY
The Mixmaster Cosmology is a valuable example. As described in §30.7\S 30.7§, it shows a singularity behavior which illustrates most of the features of the most general examples known. In particular, it shows how properties of empty space reminiscent of an elastic solid become evident near the cosmological singularity.
The mathematical path to this example, as given in this box, also illustrates several important techniques in using the variational principles for the Einstein equations to elucidate the solution of these equations. The Mixmaster example can also be used to provide simple examples of superspace ideas and of quantum formulations of the laws of gravity [Misner (1972a)].
A Generalized Kasner Model
Two generalizations must be implemented in order to progress from the Kasner example (30.1) of a cosmological singularity to the Mixmaster example. The first is to allow a more general timedependence while preserving some of the simplicity of the conditions (30.2) on the exponents p_(i)p_{i}. Note that these exponents satisfy, e.g., p_(2)-=p_{2} \equivd ln g_(22)//d ln gd \ln g_{22} / d \ln g. Therefore one is led to parametrize the 3xx33 \times 3 spatial metric as
or equivalently, (ln g)_(ij)=2alphadelta_(ij)+2beta_(ij)(\ln g)_{i j}=2 \alpha \delta_{i j}+2 \beta_{i j}, where beta_(ij)\beta_{i j} is a traceless 3xx33 \times 3 symmetric matrix, and the exponential is a matrix power series, so dete^(2beta)=1\operatorname{det} e^{2 \beta}=1 and
For the purposes of this paragraph only, define p_(ij)=d(ln g)_(ij)//d ln det gp_{i j}=d(\ln g)_{i j} / d \ln \operatorname{det} g. Then from equations (1) and (2), one computes
is an identity in view of trace beta_(ij)=0\beta_{i j}=0. The second condition on the Kasner exponents is trace (p^(2))=\left(p^{2}\right)= 1 , and becomes (dbeta_(ij)//d alpha)^(2)=6\left(d \beta_{i j} / d \alpha\right)^{2}=6 by equation (3). This is not an identity, but a consequence of the Einstein equations in empty space. For the (Bianchi Type I) metric
{:(4)ds^(2)=-dt^(2)+e^(2alpha)(e^(2beta))_(ij)dx^(i)dx^(j):}\begin{equation*}
d s^{2}=-d t^{2}+e^{2 \alpha}\left(e^{2 \beta}\right)_{i j} d x^{i} d x^{j} \tag{4}
\end{equation*}
and in the case when beta_(ij)\beta_{i j} is diagonal, the Einstein equations are,
together with a redundant equation involving T_(kk)T_{k k} and the equation T_(0k)=0T_{0 k}=0. [The stress components here refer to an orthonormal frame with basis 1 -forms omega^(i)=e^(alpha)(e^(beta))_(ij)dx^(j)\omega^{i}=e^{\alpha}\left(e^{\beta}\right)_{i j} d x^{j}.] From equation (5) one immediately derives
{:(7)rho_("aniso(t) ")=(c^(2)//16 pi G)(dbeta_(ij)//dt)^(2):}\begin{equation*}
\rho_{\text {aniso(t) }}=\left(c^{2} / 16 \pi G\right)\left(d \beta_{i j} / d t\right)^{2} \tag{7}
\end{equation*}
as a formula for the effectiveness of Type I anisotropy in contributing to the Hubble constant H=H=d alpha//dtd \alpha / d t on a basis comparable to matter energy, as in equation (30.4). Similarly, for equation (6) in the case of fluid matter (isotropic pressures), the stress terms vanish, and one obtains rho_("aniso(I) ")e^(6alpha)=\rho_{\text {aniso(I) }} e^{6 \alpha}= const., as in the equation following (30.4). The Kasner condition Sigmap_(i)^(2)=1\Sigma p_{i}{ }^{2}=1 or (dbeta_(ij)//d alpha)^(2)=6\left(d \beta_{i j} / d \alpha\right)^{2}=6 follows from equation (5) whenever T^(00)≪rho_("aniso ")T^{00} \ll \rho_{\text {aniso }}.
In the diagonal case, beta_(ij)\beta_{i j} has only two independ-
ent components, and it is convenient at times to define them explicitly by the parameterization
For these the Kasner condition (dbeta_(ij)//d alpha)^(2)=6\left(d \beta_{i j} / d \alpha\right)^{2}=6 becomes
{:(9)(dbeta_(+)//d alpha)^(2)+(dbeta_(-)//d alpha)^(2)=1:}\begin{equation*}
\left(d \beta_{+} / d \alpha\right)^{2}+\left(d \beta_{-} / d \alpha\right)^{2}=1 \tag{9}
\end{equation*}
The beta_(+-)\beta_{ \pm}are related to the Kasner exponents p_(i)p_{i} or the uu parameter of equations (30.5) by
{:[dbeta_(+)//d alpha=(1)/(2)(1-3p_(3))],[=-1+(3//2)(1+u+u^(2))^(-1)],[(10)dbeta_(-)//d alpha=(1)/(2)sqrt3(p_(1)-p_(2))],[=-(1)/(2)sqrt3(1+2u)(1+u+u^(2))^(-1)]:}\begin{align*}
d \beta_{+} / d \alpha & =\frac{1}{2}\left(1-3 p_{3}\right) \\
& =-1+(3 / 2)\left(1+u+u^{2}\right)^{-1} \\
d \beta_{-} / d \alpha & =\frac{1}{2} \sqrt{3}\left(p_{1}-p_{2}\right) \tag{10}\\
& =-\frac{1}{2} \sqrt{3}(1+2 u)\left(1+u+u^{2}\right)^{-1}
\end{align*}
Introducing Space Curvature
The first step in generalizing the Kasner metric has focused attention on the "velocity" beta^(')-=(dbeta_(+)//d alpha:}\boldsymbol{\beta}^{\prime} \equiv\left(d \beta_{+} / d \alpha\right., dbeta_(-)//d alphad \beta_{-} / d \alpha ) which is a derivative of anisotropy with respect to expansion. The effects of matter or, as will soon appear, space curvature can change the magnitude ||beta^(')||\left\|\boldsymbol{\beta}^{\prime}\right\| from the Kasner value of unity. The second step of generalization is to introduce space curvature. This one achieves in a simple example by retaining the metric components of equation (1), but employing them in a non-holonomic basis. Use the basis vectors introduced in exercises 9.13 and 9.14 on the rotation group SO(3)S O(3), whose dual 1-forms are
{:[sigma^(1)=cos psi d theta+sin psi sin theta d phi","],[(11)sigma^(2)=sin psi d theta-cos psi sin theta d phi","],[sigma^(3)=d psi+cos theta d phi]:}\begin{align*}
& \sigma^{1}=\cos \psi d \theta+\sin \psi \sin \theta d \phi, \\
& \sigma^{2}=\sin \psi d \theta-\cos \psi \sin \theta d \phi, \tag{11}\\
& \sigma^{3}=d \psi+\cos \theta d \phi
\end{align*}
to form the metric
{:(12)ds^(2)=-N^(2)dt^(2)+e^(2alpha)(e^(2beta))_(ij)sigma^(i)sigma^(j):}\begin{equation*}
d s^{2}=-N^{2} d t^{2}+e^{2 \alpha}\left(e^{2 \beta}\right)_{i j} \sigma^{i} \sigma^{j} \tag{12}
\end{equation*}
where N,alphaN, \alpha, and beta_(ij)\beta_{i j} are functions of tt only. When
Box 30.1 (continued)
alpha=0=beta_(ij)\alpha=0=\beta_{i j}, the three-dimensional space metric here reduces to the one studied in exercise 13.15, which is the metric of highest symmetry on the group space SO(3)S O(3). The simply connected covering space has the 3 -sphere topology, and is obtained by extending the range of the Euler angle psi\psi to give it a 4pi4 \pi period [SU(2):}\left[S U(2)\right. or spin (1)/(2)\frac{1}{2} covering of the rotation group]. With N=1,(1)/(2)a=e^(alpha)N=1, \frac{1}{2} a=e^{\alpha}, and beta_(ij)=0\beta_{i j}=0, one obtains from equation (12) the same metric (in different coordinates) as that treated in exercise 14.4 and in Chapter 27 in discussions of the closed Friedmann cosmological model. A non-zero value for beta_(ij)\beta_{i j} allows the 3 -sphere to have a different circumference on great circles in each of 3 mutually orthogonal principal directions, thus destroying its isotropy but not its homogeneity.
Let us consider only the case with beta_(ij)\beta_{i j} diagonal, as in equation (8). Then the T^(00)T^{00} Einstein equation becomes (with N=1N=1 as a time-coordinate condition)
is different from equation (5). This term [see equation (21.92)] is the scalar curvature of a threedimensional slice, t=t= const [which has symmetry properties known as "Bianchi Type IX" for the metric of equations (11) and (12)]. If equation (13) is interpreted in terms of an anisotropy energy density contributing, with T^(00)T^{00}, to the volume expansion alpha^(˙)^(2)\dot{\alpha}^{2}, then there are not only kinetic energy terms beta^(˙)^(2)\dot{\beta}^{2} [as in equations (5) and (7)], but also a potential energy term. This term shows that negative scalar curvature, which can be produced by anisotropy ( beta!=0\beta \neq 0 ), is equivalent to a positive potential (or "internal") energy, and suggests that empty space has properties with analogies to an elastic solid and resists shear strains. The more detailed analysis which follows shows that, near
the singularity, the scalar curvature is always negligible when positive.
Negative curvatures, however, arise in this closed universe from large shear (beta)(\beta) deformations near the singularity and become large enough to reverse one Kasner shear motion [ uu-value, etc.; equation (10)] and change it to another.
These conclusions and further details of the time-evolution of the "Mixmaster" metric (11,12)(11,12) require, in principle, the study of all the Einstein equations, not just equation (13) for T^(00)T^{00}. As described in Chapter 21, however, this T^(00)T^{00} constraint equation is central, and actually contains implicitly the full content of the Einstein equations when formulated properly.
Variational Principles
One adequate formulation, adopted here, involves treating equation (13) not as an energy equation (involving velocities), but as a Hamiltonian (involving momenta). Take the Einstein variational principle (21.15) in ADM form (21.95) and carry out the space integration, using intsigma^(1)^^sigma^(2)^^sigma^(3)=int sin theta d phi^^d theta^^d psi=(4pi)^(2)\int \sigma^{1} \wedge \sigma^{2} \wedge \sigma^{3}=\int \sin \theta d \phi \wedge d \theta \wedge d \psi=(4 \pi)^{2}, to obtain the action integral in the form
When introducing the specific form (1) and (8) for g_(ij)g_{i j}, it is convenient also to parameterize the diagonal matrix pi^(i)_(k)\pi^{i}{ }_{k} as follows:
so V(0)=0V(0)=0; and adjust the zero of alpha(alpha longrightarrow\alpha(\alpha \longrightarrowalpha-alpha_(0)\alpha-\alpha_{0} ) so that e^(2alpha)longrightarrow(6pi)^(-1)e^(2alpha)e^{2 \alpha} \longrightarrow(6 \pi)^{-1} e^{2 \alpha}. Then the metric is
{:(20)ds^(2)=-N^(2)dt^(2)+(6pi)^(-1)e^(2alpha)(e^(2beta))_(ij)sigma^(i)sigma^(j):}\begin{equation*}
d s^{2}=-N^{2} d t^{2}+(6 \pi)^{-1} e^{2 \alpha}\left(e^{2 \beta}\right)_{i j} \sigma^{i} \sigma^{j} \tag{20}
\end{equation*}
and the variational integral is
{:[(21)I= intp_(+)dbeta_(+)+p_(-)dbeta_(-)+p_(alpha)d alpha],[-(3pi//2)^(1//2)Ne^(-3alpha)Hdt","]:}\begin{align*}
I= & \int p_{+} d \beta_{+}+p_{-} d \beta_{-}+p_{\alpha} d \alpha \tag{21}\\
& -(3 \pi / 2)^{1 / 2} N e^{-3 \alpha} \mathscr{H} d t,
\end{align*}
One demands delta I=0\delta I=0 for arbitrary independent variations of p_(+-),p_(alpha),beta_(+-),alpha,Np_{ \pm}, p_{\alpha}, \beta_{ \pm}, \alpha, N to obtain the Einstein equations. From varying NN, one obtains the fundamental constraint equation K=0\mathscr{K}=0 [which would reduce to the vacuum version of equation (13) when the momenta are replaced by velocities (via equations obtained by varying the pp 's) if the coordinate condition N=1N=1 were imposed.]
ADM Hamiltonian
The standard ADM prescription for reducing this variational principle to canonical (Hamiltonian) form is to choose one of the field variables or momenta as a time-coordinate, and solve the con-
straint for its conjugate Hamiltonian. Here an obvious and satisfactory choice is to set t=alphat=\alpha, and solve K=0\mathscr{K}=0 for H_(ADM)=-p_(alpha)=[p_(+)^(2)+p_(-)^(2)+e^(4alpha)(V-1)]^(1//2)H_{\mathrm{ADM}}=-p_{\alpha}=\left[p_{+}^{2}+p_{-}^{2}+e^{4 \alpha}(V-1)\right]^{1 / 2}.
The alpha^(˙)\dot{\alpha} equation [vary p_(alpha)p_{\alpha} in equation (21)] is
The reduced, canonical, variational principle which results when equation (23) is used to eliminate p_(alpha)p_{\alpha} reads deltaI_("red ")=0\delta I_{\text {red }}=0 with
{:(26)I_(red)=intp_(+)dbeta_(+)+p_(-)dbeta_(-)-H_(ADM)d alpha:}\begin{equation*}
I_{\mathrm{red}}=\int p_{+} d \beta_{+}+p_{-} d \beta_{-}-H_{\mathrm{ADM}} d \alpha \tag{26}
\end{equation*}
and must be supplemented by equation (25).
Super-Hamiltonian
A more convenient approach here is one more closely related to the Dirac Hamiltonian methods than those of ADM. Note, however, that one does not remove the arbitrariness in the lapse function by taking it to be some specified function N(t)N(t) of the coordinates. Instead the procedure adopted here is to eliminate NN from the variational principle (21) by choosing it (coordinate condition!) to be some chosen function of the field variables and momenta, N=N(alpha,beta_(+-),p_(alpha),p_(+-))N=N\left(\alpha, \beta_{ \pm}, p_{\alpha}, p_{ \pm}\right). Any such choice, inserted in equation (21), leaves a variational integral in canonical Hamiltonian form. The content of this new variational principle becomes equivalent to the original one only when supplemented by the constraint
which can no longer be derived from the variational principle. [The other Euler-Lagrange equations for these two principles differ only by terms proportional to K\mathscr{K}, and thus are equivalent when
Box 30.1 (continued)
H=0\mathscr{H}=0 is imposed on the initial conditions.] The choice
is obvious and convenient. It makes H\mathscr{H} become a super-Hamiltonian in the resulting variational principle
{:(29)I=intp_(+)dbeta_(+)+p_(-)dbeta_(-)+p_(alpha)d alpha-Hd lambda",":}\begin{equation*}
I=\int p_{+} d \beta_{+}+p_{-} d \beta_{-}+p_{\alpha} d \alpha-\mathscr{H} d \lambda, \tag{29}
\end{equation*}
where t-=lambdat \equiv \lambda has been written to label the specific time-coordinate choice that equation (28) implies.
Mixmaster Dynamics
If matter terms with no additional degrees of freedom are included, the super-Hamiltonian in equation (29) is modified simply. For an example, choose
for the energy density of matter in a frame with time-axis e_( hat(0))=N^(-1)(del//del t)\boldsymbol{e}_{\hat{0}}=N^{-1}(\partial / \partial t). The two terms represent a nonrelativistic perfect fluid ( rho propV^(-1)\rho \propto V^{-1} ) and a radiation fluid ( rho propV^(-4//3)\rho \propto V^{-4 / 3} ), respectively, and lead to
This Hamiltonian, with its simple quadratic momentum dependence, differs in only two ways from the Hamiltonians of elementary mechanics, namely, (1) in the sign of the p_(alpha)^(2)p_{\alpha}{ }^{2} term and (2) in the detailed shape of the "potential" term as function of alpha\alpha and beta_(+-)\beta_{ \pm}, the study of which reduces to a study of the function V(beta)V(\beta). Hamilton's equations, from varying alpha,beta_(+-),p_(alpha)\alpha, \beta_{ \pm}, p_{\alpha}, and p_(+-)p_{ \pm}in equation (29), yield
Thus the sign of the p_(alpha)^(2)p_{\alpha}{ }^{2} term causes alpha\alpha to accelerate toward (rather than away from) higher values of the "potential" terms e^(4alpha)(V-1)+mue^(3alpha)+Gammae^(2alpha)e^{4 \alpha}(V-1)+\mu e^{3 \alpha}+\Gamma e^{2 \alpha}. When |V|≪1|V| \ll 1 (small anisotropy), equation (33) is identical to its form in the isotropic Friedmann model, and allows a deceleration only when alpha\alpha is large enough that the positive curvature term ( -e^(4alpha)-e^{4 \alpha} ) dominates over matter ( mue^(3alpha)\mu e^{3 \alpha} ) and radiation (Gammae^(2alpha))\left(\Gamma e^{2 \alpha}\right). Near the singularity (alpha longrightarrow-oo)(\alpha \longrightarrow-\infty), the positive curvature term is always negligible compared to radiation and matter.
For studies of the singularity behavior, it is sufficient to study the simplified super-Hamiltonian
since the other terms obviously vanish for alpha longrightarrow\alpha \longrightarrow-oo-\infty. This form retains only the VV term in ^(3)R_(IX)={ }^{3} R_{\mathrm{IX}}=(3)/(2)e^(-2alpha)(1-V)\frac{3}{2} e^{-2 \alpha}(1-V), which dominates when the curvature of this closed universe becomes negative, V≫1V \gg 1. If the term in V(beta)V(\beta) were also negligible, then K=-p_(alpha)^(2)+p_(+)^(2)+p_(-)^(2)\mathscr{K}=-p_{\alpha}{ }^{2}+p_{+}{ }^{2}+p_{-}{ }^{2} would make each p_(alpha),p_(+-)p_{\alpha}, p_{ \pm}constant, giving the Kasner behavior with
and |d beta//d alpha|^(2)=1|d \boldsymbol{\beta} / d \alpha|^{2}=1 as expected (since matter and curvature have been neglected). To proceed further, a study of V(beta)V(\beta) is required, based on equations (19) and (8), and their immediate consequence:
One finds that V(beta)V(\beta) is a positive definite "potential well" which has the same symmetries as an equi-
lateral triangle in the beta_(+)beta_(-)\beta_{+} \beta_{-}plane. Near the origin, beta_(+-)=0\beta_{ \pm}=0, the equipotentials are circles, since
These two asymptotic forms, together with the triangular symmetry, give a complete asymptotic description of V(beta)V(\beta), as sketched in the figure, where on successive levels separated by Delta beta=1\Delta \beta=1, the potential VV increases by a factor of e^(8)=e^{8}=3xx10^(3)3 \times 10^{3}.
"Bounce" Interrupts Kasner-like Steps Toward the Singularity
The dominant feature of the V(beta)V(\beta) potential is evidently its steep (exponential) triangular walls, with equation (37) representing the typical one for study. Under the influence of this potential wall, the evolution of this model universe is governed by the super-Hamiltonian
If alpha longrightarrow-oo\alpha \longrightarrow-\infty with dbeta_(+)//d alpha > 1//2d \beta_{+} / d \alpha>1 / 2 [recall dbeta_(+)//d alphad \beta_{+} / d \alpha== const., |d beta//d alpha|=1|d \boldsymbol{\beta} / d \alpha|=1, when the last term in (39) is small], then the potential term grows and will eventually become large enough to influence the motion. A simple "Lorentz" transformation, suggested by the superspace metric (coefficients of the
Some equipotentials, V(beta)=V(\beta)= constant, are shown for the function defined in equation (35). Equipotentials near the origin of the beta\beta-plane are closed curves for V < 1V<1 and are omitted here.
Box 30.1 (continued)
quadratic in the momenta) simplifies the computation further. Set
For this super-Hamiltonian both bar(p)_(alpha)\bar{p}_{\alpha} and p_(-)p_{-}are constants of motion, whereas the bar(beta)_(+)\bar{\beta}_{+}-Hamiltonian, bar(p)_(+)^(2)+(1)/(3)e^(-4sqrt3 bar(beta)_(+))\bar{p}_{+}{ }^{2}+\frac{1}{3} e^{-4 \sqrt{3} \bar{\beta}_{+}}, represents a simple bounce against a one-dimensional potential wall with the initial and final values of bar(p)_(+)\bar{p}_{+}different only in sign. The behavior of the anisotropy parameters beta_(+-)\beta_{ \pm}near the singularity thus consists of a simple Kasner step (where dbeta_(+-)//d alpha=d \beta_{ \pm} / d \alpha= const., with the dbeta_(+)//d alpha >= (1)/(2)d \beta_{+} / d \alpha \geq \frac{1}{2}, or conditions equivalent by symmetry, satisfied relative to one of the three walls), followed by a bounce against that wall, beginning a new Kasner step with other Kasner parameters. [The most detailed description of this behavior and its relation to more general cosmological models can be found in Belinsky, Khalatnikov, and Lifshitz (1970)-see also the briefer report, Khalatnikov and Lifshitz (1970)-using quite different methods. For detailed developments by Hamiltonian methods, which supercede the partial Lagrangian methods of Misner (1969b), see Misner (1970, 1972a), and Ryan (1972a,b).]
Steady-State, Quasiperiodic Infinity of "Bounces" Approaching the Singularity
Some comprehensive features of the singularity behavior, involving many Kasner-like steps, can be exhibited by another transformation of the parameter space (superspace) of the metric field.
The transformation introduces a "radial" tt coordinate out from the origin of alphabeta_(+-)\alpha \beta_{ \pm}space, while respecting the metric properties of this superspace implied by the form of the super-Hamiltonian. Thus one defines (for any constant alpha_(0)\alpha_{0} )
The advantage of this transformation is that in the limit t longrightarrow oo(alpha longrightarrow-oot \longrightarrow \infty(\alpha \longrightarrow-\infty, singularity) the potential terms become, in first approximation, independent of tt. Thus equation (37) gives, for one potential wall,
For t longrightarrow oot \longrightarrow \infty this expression evidently tends to either zero or infinity, depending on the sign of the expression in parentheses. Therefore define the asymptotic potential walls by
in the sector |phi-pi| < pi//3|\phi-\pi|<\pi / 3, and equivalent formulae in which phi\phi is replaced by phi+-(2pi//3)\phi \pm(2 \pi / 3) for the other sides of the triangle. Consequently, an asymptotic approximation to the super-Hamiltonian is
where V^(')(zeta,phi)V^{\prime}(\zeta, \phi) vanishes inside the asymptotic walls (44) and equals +oo+\infty outside. Because the remaining tt-dependence is a common factor in (45), a simple change of independent variable e^(-2t)d lambda=e^{-2 t} d \lambda=dlambda^(')d \lambda^{\prime} in equation (29)-equivalent to the choice
in place of equation (28)-gives a new superHamiltonian K^(')=e^(2t)K\mathscr{K}^{\prime}=e^{2 t} \mathscr{K} with the variational integral
{:(47)I=intp_(t)dt+p_(s)d xi+p_(phi)d phi-K^(')dlambda^(').:}\begin{equation*}
I=\int p_{t} d t+p_{s} d \xi+p_{\phi} d \phi-\mathscr{K}^{\prime} d \lambda^{\prime} . \tag{47}
\end{equation*}
one immediately sees that p_(t)p_{t} is a constant of motion, and that the "bouncing" of the zeta phi\zeta \phi values within the asymptotic potential walls is a stationary, quasi-periodic process in this time-coordinate lambda^(')\lambda^{\prime} (or tt, since dt//dlambda^(')=-p_(t)=d t / d \lambda^{\prime}=-p_{t}= const). [More detailed studies based on this asymptotic superHamiltonian show that the motion is even ergodic, with zeta phi\zeta \phi approaching arbitrarily close to any given value infinitely many times as t longrightarrow oot \longrightarrow \infty; see Chitre (1972a).]
Summary
One has found the singularity behavior in this Mixmaster example to be extraordinarily active. In the simple Kasner singularity, two axes collapse, but the third is stretched in a simple tidal deformation accompanied by volume compression. But in the Mixmaster example, every such collapse attempt is defeated by the high negative curvature it implies. Or rather it is diverted to another attempt as compression continues inexorably, but the tidal deformations attempt first one configuration, then another, in an infinitely recurring probing of all possible configurations.
Speculations on Time and the Singularity
The cosmological singularity (in all examples where its character is not known to be unstable) involves infinite curvature and infinite density. One's abhorence of such a theoretical prediction is particularly heightened by the correlative prediction that these infinities occurred at a finite proper time in the past, and would-if they
recur-occur again at some finite proper time in the future. The singularity prediction would be more tolerable if the infinite densities could be removed to the infinitely distant past. The universe could then, as now, find its natural state to be one of expansion, so every finite density will have been experienced at some suitably remote past time, but infinite density becomes a formal abstraction never realized in the course of evolution.
To push infinite curvature out of the finite past might be achieved in two ways. It is not known which, if either, works. One way is to change the physical laws which require the singularity, changing them perhaps only in obvious and desirable ways, such as stating the laws of gravity in a proper quantum language. Computations of quantum geometry are not yet definitive, however, and some (perhaps inadequate) approximations [Misner (1972a)] do not remove the singularity problem.
Another way to discard the singularity is to accept the mathematics of the classical Einstein equations, but reinterpret it in terms of an infinite past time. There are, of course, simple and utterly inadequate ways to do this by arbitrary coordinate transformations such as t=ln taut=\ln \tau which change a tau=0\tau=0 singularity into one at t=-oot=-\infty. But an arbitrary coordinate is without significance. The problem is that the singularity occurs at a finite proper time in the past, and proper time is the most physically significant, most physically real time we know. It corresponds to the ticking of physical clocks and measures the natural rhythms of actual events. To reinterpret finite past time as infinite, one must attack proper time on precisely these grounds, and claim it is inadequately physical. On a local basis, where special relativity is valid, no challenge to the physical significance of proper time can succeed. It is on a more global scale that the physical primacy of proper time needs to be reviewed.
"The cosmological singularity occurred ten thousand million years ago." In this statement, take time to mean the proper time along the world line of the solar system, ephemeris time. Then the statement would have a most direct physical sig-
Box 30.1 (continued)
nificance if it meant that the Earth had completed 10^(10)10^{10} orbits about the sun since the beginning of the universe. But proper time is not that closely tied to actual physical phenomena. The statement merely implies that those 5xx10^(9)5 \times 10^{9} orbits which the earth may have actually accomplished give a standard of time which is to be extrapolated in prescribed ways, thus giving theoretical meaning to the other 5xx10^(9)5 \times 10^{9} years which are asserted to have preceeded the formation of the solar system.
A hardier standard clock changes the details of the argument, but not its qualitative conclusion. To interpret 10^(10)10^{10} years in terms of SI (Systeme Internationale) seconds assigns a past history containing some 3xx10^(27)3 \times 10^{27} oscillations of a hyperfine transition in neutral Cesium. But again the critical early ticks of the clock (needed to locate the singularity in time by actual physical events) are missing. The time needed for stellar nucleosynthesis to produce the first Cesium disqualifies this clock on historical grounds, and the still earlier high temperatures nearer the singularity would have ionized all Cesium even if this element had predated stars.
Thus proper time near the singularity is not a direct counting of simple and actual physical phenomena, but an elaborate mathematical extrapolation. Each actual clock has its "ticks" discounted by a suitable factor- 3xx10^(7)3 \times 10^{7} seconds per orbit from the Earth-sun system, 1.1 xx10^(-10)1.1 \times 10^{-10} seconds per oscillation for the Cesium transition, etc. Since no single clock (because of its finite size and strength) is conceivable all the way back to the singularity, a statement about the proper time since the singularity involves the concept of an infinite sequence of successively smaller and sturdier clocks with their ticks then discounted and
added. "Finite proper time," then, need not imply that any finite sequence of events was possible. It may describe a necessarily infinite number of events ("ticks") in any physically conceivable history, converted by mathematics into a finite sum by the action of a non-local convergence factor, the "discount" applied to convert "ticks" into "proper time."
Here one has the conceptual inverse of Zeno's paradox. One rejects Zeno's suggestion that a single swing of a pendulum is infinitely complicatedbeing composed of a half period, plus a quarter period, plus 2^(-n)2^{-n} ad infinitum-because the terms in his infinite series are mathematical abstractions, not physically achieved discrete acts in a drama that must be played out. By a comparable standard, one should ignore as a mathematical abstraction the finite sum of the proper-time series for the age of the universe, if it can be proved that there must be an infinite number of discrete acts played out during its past history. In both cases, finiteness would be judged by counting the number of discrete ticks on realizable clocks, not by assessing the weight of unrealizable mathematical abstractions.
Whether the universe is infinitely old by this standard remains to be determined. The quantum influences, in particular, remain to be calculated. The decisive question is whether each presentepoch event is subject to the influence of infinitely many previous discrete events. In that case statistical assumptions (large numbers, random phases, etc.) could enter in stronger ways into theories of cosmology. The Mixmaster cosmological model does have an infinite past history in this sense, since each "bounce" from one Kasner-like motion to another is a recognizable cosmological event, of which infinitely many must be realized between any finite epoch and the singularity.
§30.8. HORIZONS AND THE ISOTROPY OF the microwave background
The fundamental cosmological question-"Must a universe that is born chaotic necessarily become as homogeneous and isotropic as our universe is, and do so before life evolves?"-entails one further issue. This issue is horizons. As was discussed in §27.10, at any given epoch in the expansion of a Friedmann universe (e.g., the present epoch), there may be significant portions of the universe from which no light signal or other causally propagating influence will have yet reached Earth in the time available since the initial singularity. "If we should live so long," the question would arise, "will the new portions of the universe which first come into view during the next ten billion years look statistically identical to the neighboring portions which are already being seen?"
Fortunately, this question need not be posed only for the future. It can be asked as of some past time, and the answer then is yes. Microwave background radiation arrives at the earth from all directions in the sky with very nearly the same temperature. [The data of Boughn, Fram, and Partridge (1971) and of Conklin (1969) show Delta T//T <= 0.004\Delta T / T \leqq 0.004.] The plasma that emitted the microwave radiation coming to earth from one direction in the sky had not been able, before the epoch of emission, to communicate causally with the plasma emitting the radiation that arrives from other directions. If one adopts a Friedmann model of the universe, then different sectors of the microwave sky are disjoint from each other in this sense if they are separated from each other by more than 30^(@)30^{\circ}, even if the microwaves were emitted as recently as z=7z=7. (The critical angle is much smaller if the microwaves were last scattered at z=1,000z=1,000.) From this, one concludes that the foundations for the homogeneity and isotropy of the universe were laid long before the universe became approximately Friedmann, for if statistical homogeneity and isotropy of the universe had not already been achieved at the longest wavelengths earlier, these horizon limitations would have prevented any further synchronization of conditions over large scales while the universe was in a nearly Friedmann state, and small amplitude ( 10%10 \% ) deviations from isotropy should be observed now.
The mixmaster universe received its name from the hope that it could contribute to the solution of this problem. The very large uu values that occur sporadically an infinite number of times near the singularity in a mixmaster universe give a geometry close to that of the Kasner model with p_(1)=1,p_(2)=p_(3)=0p_{1}=1, p_{2}=p_{3}=0. This model can be written in the form
{:(30.6)ds^(2)=e^(2eta)(-deta^(2)+dx^(2))+dy^(2)+dz^(2):}\begin{equation*}
d s^{2}=e^{2 \eta}\left(-d \eta^{2}+d x^{2}\right)+d y^{2}+d z^{2} \tag{30.6}
\end{equation*}
where eta=ln t\eta=\ln t. If this metric is converted into a closed-universe model by interpreting x,y,zx, y, z as angle coordinates each with period 4pi4 \pi, then one sees that a light ray can circumnavigate the universe in the xx-direction in a time interval Delta eta=4pi\Delta \eta=4 \pi, which corresponds to a volume expansion by a factor sqrt(-g_(1))//sqrt(-g_(2))=e^(4pi)\sqrt{-g_{1}} / \sqrt{-g_{2}}=e^{4 \pi}. Unfortunately, a quantitative analysis of the degree and frequency with which the mixmaster universe achieves this specific Kasner form suggests that the horizon breaking
Horizons in a Friedmann universe
Observed isotropy of microwave radiation proves foundations for homogeneity were laid before universe became Friedmann-like
What made the universe homogeneous and isotropic?
(1) Mixmaster oscillations? probably not
(2) particle creation near singularity?
is inadequate to explain the present state of the universe [Doroshkevich, Lukash, and Novikov (1971); Chitre (1972)]. It may turn out that particle creation near the singularity can solve this horizon question, as well as provide for the dissipation of anisotropy. Hope is provided by the fact that particle creation, when described in purely classical terms, has some acausal appearances, even though it is a strictly causal process at the quantum level [Zel'dovich (1972)].
GRAVITATIONAL COLLAPSE AND BLACK HOLES
Wherein the reader is transported to the land of black holes, and encounters colonies of static limits, ergospheres, and horizons-behind whose veils are hidden gaping, ferocious singularities.
синотен 31
SCHWARZSCHILD GEOMETRY
§31.1. INEVITABILITY OF COLLAPSE FOR MASSIVE STARS
There is no equilibrium state at the endpoint of thermonuclear evolution for a star containing more than about twice the number of baryons in the sun (A > A_(max)∼2A_(o.))\left(A>A_{\max } \sim 2 A_{\odot}\right). This is one of the most surprising-and disturbing-consequences of the discussion in Chapter 24. Stated differently: A star with A > A_(max)∼2A_(o.)A>A_{\max } \sim 2 A_{\odot} must eject all but A_("max ")A_{\text {max }} of its baryons-e.g., by nova or supernova explosions-before settling down into its final resting state; otherwise there will be no final resting state for it to settle down into.
What is the fate of a star that fails to eject its excess baryons before nearing the endpoint of thermonuclear evolution? For example, after a very massive supernova explosion, what will become of the collapsed degenerate-neutron core when it contains more than A_("max ")A_{\text {max }} baryons? Such a supercritical mass cannot explode, since it is gravitationally bound and it has no more thermonuclear energy to release. Nor can it reach a static equilibrium state, since there exists no such state for so large a mass. There remains only one alternative; the supercritical mass must collapse through its "gravitational radius," r=2Mr=2 M, leaving behind a gravitating "black hole" in space.
The phenomenon of collapse through the gravitational radius, as described by classical general relativity, will be the subject of the next chapter. However, before tackling it, one must understand more fully than heretofore the Schwarzschild spacetime geometry, which surrounds black holes and collapsing stars as well as static stars.
This chapter will concern itself with two topics that, at first sight, appear to be disconnected. One is the fall of a test particle in a preexisting Schwarzschild geometry, which is regarded as static, but can also be visualized as all that remains of a star that underwent collapse some time ago. The second topic is the physical
This chapter, on Schwarzschild geometry, is key preparation for understanding gravitational collapse (next chapter) and black holes (following chapter)
character of this geometry, regarded in and by itself. For the exploration of this geometry, the test particle serves as the best of all explorers. But the test particle may also be regarded in another light. It can be viewed as a rag-tag johnny-comelately piece of the matter of the falling star. Regarded in this way, it provides the simplest of all illustrations of an asymmetry in the distribution of mass of a collapsing star. That this asymmetry irons itself out will therefore give one some preliminary insight into how more complicated asymmetries also iron themselves out. In brief, the motion of the test particle and the dynamics of the Schwarzschild geometry (for this geometry will prove to be dynamic), two apparently different problems, have the happy ability to throw light on each other.
§31.2. THE NONSINGULARITY OF THE GRAVITATIONAL RADIUS
The Schwarzschild spacetime geometry
{:(31.1)ds^(2)=-(1-(2M)/(r))dt^(2)+(dr^(2))/(1-2M//r)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-\left(1-\frac{2 M}{r}\right) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{31.1}
\end{equation*}
appears to behave badly near r=2Mr=2 M; there g_(tt)g_{t t} becomes zero, and g_(rr)g_{r r} becomes infinite. However, one cannot be sure without careful study whether this pathology in the line element is due to a pathology in the spacetime geometry itself, or merely to a pathology of the (t,r,theta,phi)(t, r, \theta, \phi) coordinate system near r=2Mr=2 M. (As an example of a coordinate-induced pathology, consider the neighborhood of theta=0\theta=0 on one of the invariant spheres, t=t= const and r=r= const. There g_(phi phi)g_{\phi \phi} becomes zero because the coordinate system behaves badly; however, the intrinsic, coordinate-independent geometry of the sphere is well-behaved there. For another example, see Figure 1.4.
The worrisome region of the Schwarzschild geometry, r=2Mr=2 M, is called the "gravitational radius," or the "Schwarzschild radius," or the "Schwarzschild surface," or the "Schwarzschild horizon," or the "Schwarzschild sphere." It is also called the "Schwarzschild singularity" in some of the older literature; but that is a misnomer, since, as will be shown, the spacetime geometry is not singular there.
To determine whether the spacetime geometry is singular at the gravitational radius, send an explorer in from far away to chart it. For simplicity, let him fall freely and radially into the gravitational radius, carrying his orthonormal tetrad with him as he falls. His trajectory through spacetime ["parabolic orbit"; radial geodesic of metric (31.1)] is
[See $25.5\$ 25.5 and especially equation (25.38) for derivation and discussion.] One obtains the rr coordinate of the explorer in terms of the proper time measured on a clock
he carries, r(tau)r(\tau), by inverting the first equation; one finds his rr coordinate in terms of coordinate time, r(t)r(t), by inverting the second equation.
Of all the features of the traveler's trajectory, one stands out most clearly and disturbingly: to reach the gravitational radius, r=2Mr=2 M, requires a finite lapse of proper time, but an infinite lapse of coordinate time:
{:[r//2M=1-(tau+" constant ")//2M quad" when near "r=2M],[(31.3)r//2M=1+" constant "xx exp(-t//2M)quad" in limit as "t longrightarrow oo]:}\begin{gather*}
r / 2 M=1-(\tau+\text { constant }) / 2 M \quad \text { when near } r=2 M \\
r / 2 M=1+\text { constant } \times \exp (-t / 2 M) \quad \text { in limit as } t \longrightarrow \infty \tag{31.3}
\end{gather*}
(see Fig. 25.5.) Of course, proper time is the relevant quantity for the explorer's heart-beat and health. No coordinate system has the power to prevent him from reaching r=2Mr=2 M. Only the coordinate-independent geometry of spacetime could possibly do that; and equation (31.3) shows it does not!
Let the explorer approach and reach r=2Mr=2 M, then. What spacetime geometry does he measure there? Is it singular or nonsingular? Restated in terms of measurements, do infinite tidal gravitational forces tear the traveler apart and crush him as he approaches r=2Mr=2 M, or does he feel only finite tidal forces which in principle his body can withstand?
The tidal forces felt by the explorer as he passes a given radius rr are measured by the components of the Riemann curvature tensor with respect to his orthonormal frame there (equation of geodesic deviation). To calculate those curvature components at rr, proceed in two steps. (1) Calculate the components, not in the traveler's frame, but rather in the "static" orthonormal frame
located at the event through which he is passing; the result [obtainable from equations (14.50) and (14.51) by setting {:e^(2phi)=e^(-2A)=1-2M//r]\left.e^{2 \phi}=e^{-2 A}=1-2 M / r\right] is
all other R_( hat(alpha) hat(beta) hat(gamma) hat(delta))R_{\hat{\alpha} \hat{\beta} \hat{\gamma} \hat{\delta}} vanish except those obtainable from the above by symmetries of Riemann.
(2) Calculate the components in the explorer's frame by applying to the "staticframe" components (31.4b) the appropriate transformation-for r > 2Mr>2 M, a Lorentz boost in the e_( hat(r))\boldsymbol{e}_{\hat{r}} direction with ordinary velocity v^( hat(r))v^{\hat{r}}; for r < 2Mr<2 M, not a "boost," but a transformation given by the standard boost formula (Box 2.4) with v^( hat(r)) > 1v^{\hat{r}}>1. Here
{:(31.5)v^( hat(r))=((g_(rr))^(1//2)dr)/((-g_(tt))^(1//2)dt)=(dr//dt)/(1-2M//r)=-((2M)/(r))^(1//2):}\begin{equation*}
v^{\hat{r}}=\frac{\left(g_{r r}\right)^{1 / 2} d r}{\left(-g_{t t}\right)^{1 / 2} d t}=\frac{d r / d t}{1-2 M / r}=-\left(\frac{2 M}{r}\right)^{1 / 2} \tag{31.5}
\end{equation*}
The amazing result (a consequence of special algebraic properties of the Schwarzschild geometry, and somewhat analogous to what happens-or, rather, does not hap-
An infalling observer reaches r=2Mr=2 M in finite proper time but infinite coordinate time
pen-to the components of the electromagnetic field, E\boldsymbol{E} and B\boldsymbol{B}, when they are both parallel to a boost) is this: all the components of Riemann are left completely unaffected by the boost. If e_( hat(rho))\boldsymbol{e}_{\hat{\rho}} is the traveler's radial basis vector, and e_( hat(tau))=u\boldsymbol{e}_{\hat{\tau}}=\boldsymbol{u} is his time basis vector, then
(See exercise 31.1.)
The payoff of this calculation: according to equations (31.6), none of the components of Riemann in the explorer's orthonormal frame become infinite at the gravitational radius. The tidal forces the traveler feels as he approaches r=2Mr=2 M are finite; they do not tear him apart-at least not when the mass MM is sufficiently great, because at r=2Mr=2 M the typical non-zero component R_( hat(alpha) hat(beta) hat(gamma)delta)R_{\hat{\alpha} \hat{\beta} \hat{\gamma} \delta} of the curvature tensor is of the order 1//M^(2)1 / M^{2}. The gravitational radius is a perfectly well-behaved, nonsingular region of spacetime, and nothing there can prevent the explorer from falling on inward.
By contrast, deep inside the gravitational radius, at r=0r=0, the traveler must encounter infinite tidal forces, independently of the route he uses to reach there. One says that " r=0r=0 is a physical singularity of spacetime." To see this, one need only calculate from equation (31.4b) or (31.6) the "curvature invariant":
Box 31.1 THE "SCHWARZSCHILD SINGULARITY": HISTORICAL REMARKS
Although Eddington (1924) was the first to construct a coordinate system that is nonsingular at r=2Mr=2 M, he seems not to have recognized the significance of his result. Lemaître (1933c, especially p. 82) appears to have been the first to recognize that the so-called "Schwarzschild singularity" at r=2Mr=2 M is not a singularity. He wrote, "La singularité du champ de Schwarzschild est donc une singularité fictive, analogue à celle qui se présentait a l'horizon du centre dans la forme originale de l'univers de de Sitter". He also provided a coordinate system to go through r=2Mr=2 M. However, his coordinate system, like Eddington's, covered only half of the Schwarzschild geometry:
regions I and II of Figure 31.3. Synge (1950) was the first to discover the incompleteness in the Eddington and Lemaittre coordinate systems, and to provide coordinates that cover the entire geometry (regions I, II, III, IV of Figure 31.3). Fronsdal (1959), unaware of Synge's work, rediscovered the global structure of the Schwarzschild geometry by means of embedding diagrams and calculations. The coordinate system that provides maximum insight into the Schwarzschild geometry is the one generally known as the Kruskal-Szekeres coordinate system. It was constructed independently by Kruskal (1960) and by Szekeres (1960).
In every local Lorentz frame this will be a sum of products of curvature components, and it will have the same value 48M^(2)//r^(6)48 M^{2} / r^{6}. Thus, in every local Lorentz frame, including the traveler's, Riemann will have one or more infinite components as r longrightarrow0r \longrightarrow 0; i.e., tidal forces will become infinite.
Exercise 31.1. TIDAL FORCES ON INFALLING EXPLORER
(a) Carry out the details of the derivation of the Riemann tensor components (31.6).
(b) Calculate, roughly, the critical mass M_("crit ")M_{\text {crit }} such that, if M > M_("crit ")M>M_{\text {crit }} the explorer's body (a human body made of normal flesh and bones) can withstand the tidal forces at r=2Mr=2 M, but if M < M_("crit ")M<M_{\text {crit }} his body is mutilated by them. [Answer: M_("crit ")∼1000M_(o.)M_{\text {crit }} \sim 1000 M_{\odot}. Evidently, if M∼M_(o.)M \sim M_{\odot} the physicist should transform himself into an ant before taking the plunge! For details see §32.6\S 32.6§.]
At r=0r=0 the curvature is infinite
EXERCISE
§31.3. BEHAVIOR OF SCHWARZSCHILD COORDINATES AT r=2Mr=2 M
Since the spacetime geometry is well behaved at the gravitational radius, the singular behavior there of the Schwarzschild metric components, g_(tt)=-(1-2M//r)g_{t t}=-(1-2 M / r) and g_(rr)=(1-2M//r)^(-1)g_{r r}=(1-2 M / r)^{-1}, must be due to a pathology there of the Schwarzschild coordinates t,r,theta,phit, r, \theta, \phi. Somehow one must find a way to get rid of that pathology-i.e., one must construct a new coordinate system from which the pathology is absent. Before doing this, it is helpful to understand better the precise nature of the pathology.
The most obvious pathology at r=2Mr=2 M is the reversal there of the roles of tt and rr as timelike and spacelike coordinates. In the region r > 2Mr>2 M, the tt direction, del//del t\partial / \partial t, is timelike (g_(tt) < 0)\left(g_{t t}<0\right) and the rr direction, del//del r\partial / \partial r, is spacelike (g_(rr) > 0)\left(g_{r r}>0\right); but in the region r < 2M,del//del tr<2 M, \partial / \partial t is spacelike (g_(tt) > 0)\left(g_{t t}>0\right) and del//del r\partial / \partial r is timelike (g_(rr) < 0)\left(g_{r r}<0\right).
What does it mean for rr to "change in character from a spacelike coordinate to a timelike one"? The explorer in his jet-powered spaceship prior to arrival at r=2Mr=2 \mathrm{M} always has the option to turn on his jets and change his motion from decreasing rr (infall) to increasing rr (escape). Quite the contrary is the situation when he has once allowed himself to fall inside r=2Mr=2 M. Then the further decrease of rr represents the passage of time. No command that the traveler can give to his jet engine will turn back time. That unseen power of the world which drags everyone forward willy-nilly from age twenty to forty and from forty to eighty also drags the rocket in from time coordinate r=2Mr=2 M to the later value of the time coordinate r=0r=0. No human act of will, no engine, no rocket, no force (see exercise 31.3) can make time stand still. As surely as cells die, as surely as the traveler's watch ticks away "the unforgiving minutes," with equal certainty, and with never one halt along the way, rr drops from 2M2 M to 0 .
At r=2Mr=2 M, where rr and tt exchange roles as space and time coordinates, g_(tt)g_{t t} vanishes while g_(rr)g_{r r} is infinite. The vanishing of g_(tt)g_{t t} suggests that the surface r=2Mr=2 M, which
Nature of the coordinate pathology at r=2Mr=2 \mathrm{M} :
(1) tt and rr reverse roles as timelike and spacelike coordinates
(2) the region r=2Mr=2 M, -oo < t < +oo-\infty<t<+\infty is two-dimensional rather than three
(3) radial geodesics reveal that the regions r=2Mr=2 M, t=+-oot= \pm \infty are "finite" parts of spacetime
appears to be three-dimensional in the Schwarzschild coordinate system (-oo < t <(-\infty<t<+oo,0 < theta < pi,0 < phi < 2pi+\infty, 0<\theta<\pi, 0<\phi<2 \pi ) has zero volume and thus is actually only two-dimensional, or else is null; thus,
{:[(31.8)int_(r=2M)|g_(tt)g_(theta theta)g_(phi phi)|^(1//2)dtd theta d phi=0],[int_((r=2M,t=" const "))|g_(theta theta)g_(phi phi)|^(1//2)d theta d phi=4pi(2M)^(2)]:}\begin{align*}
& \int_{r=2 M}\left|g_{t t} g_{\theta \theta} g_{\phi \phi}\right|^{1 / 2} d t d \theta d \phi=0 \tag{31.8}\\
& \int_{(r=2 M, t=\text { const })}\left|g_{\theta \theta} g_{\phi \phi}\right|^{1 / 2} d \theta d \phi=4 \pi(2 M)^{2}
\end{align*}
The divergence of g_(rr)g_{r r} at r=2Mr=2 M does not mean that r=2Mr=2 M is infinitely far from all other regions of spacetime. On the contrary, the proper distance from r=2Mr=2 M to a point with arbitrary rr is
{:(31.9)int_(2M)^(r)|g_(rr)|^(1//2)dr={[[r(r-2M)]^(1//2)+2M ln|(r//2M-1)^(1//2)+(r//2M)^(1//2)|],[" when "r > 2M],[-2Mcot^(-1)[r^(1//2)//(2M-r)^(1//2)]-[r(2M-r)]^(1//2)],[" when "r < 2M]:}:}\int_{2 M}^{r}\left|g_{r r}\right|^{1 / 2} d r=\left\{\begin{array}{r}
{[r(r-2 M)]^{1 / 2}+2 M \ln \left|(r / 2 M-1)^{1 / 2}+(r / 2 M)^{1 / 2}\right|} \tag{31.9}\\
\text { when } r>2 M \\
-2 M \cot ^{-1}\left[r^{1 / 2} /(2 M-r)^{1 / 2}\right]-[r(2 M-r)]^{1 / 2} \\
\text { when } r<2 M
\end{array}\right.
which is finite for all 0 < r < oo0<r<\infty.
Just how the region r < 2Mr<2 M is physically connected to the region r > 2Mr>2 M can be discovered by examining the radial geodesics of the Schwarzschild metric. Focus attention, for concreteness, on the trajectory of a test particle that gets ejected from the singularity at r=0r=0, flies radially outward through r=2Mr=2 M, reaches a maximum radius r_(max)r_{\max } ("top of orbit") at proper time tau=0\tau=0 and coordinate time t=0t=0, and then falls back down through r=2Mr=2 M to r=0r=0. The solution of the geodesic equation for such an orbit was derived in §25.5\S 25.5§ and described in Figure 25.3. It has the "cycloid form" (with the parameter eta\eta running from -pi-\pi to +pi+\pi ),
Figure 31.1 plots this orbit in the r,tr, t-coordinate plane (curve F-F^(')-F^('')F-F^{\prime}-F^{\prime \prime} ), along with several other types of radial geodesics.
Every radial geodesic except a "set of geodesics of measure zero" crosses the gravitational radius at t=+oot=+\infty (or at t=-oot=-\infty, or both), according to Figure 31.1 and the calculations behind that figure (exercises for the student! See Chapter 25). One therefore suspects that all the physics at r=2Mr=2 M is consigned to t=+-oot= \pm \infty by reason of some unhappiness in the choice of the Schwarzschild coordinates. A better coordinate system, one begins to believe, will take these two "points at infinity" and
Figure 31.1.
Typical radial geodesics of the Schwarzschild geometry, as charted in Schwarzschild coordinates (schematic). FF^(')F^('')F F^{\prime} F^{\prime \prime} [see equations (31.10)] is the timelike geodesic of a test particle that starts at rest at r=5.2Mr=5.2 \mathrm{M} and falls straight in, arriving in a finite proper time at the singularity r=0r=0 (zig-zag marking). The unhappiness of the Schwarzschild coordinate system shows in two ways: (1) in the fact that tt goes to oo\infty partway through the motion; and (2) in the fact that tt thereafter decreases as tau\tau (not shown) continues to increase. The course of the same trajectory prior to t=0t=0 may be constructed by reflecting the diagram in the horizontal axis ("time inversion"). The time-reversed image of FF " marks the ejection of the test particle from the singularity. AA^(')A^('')A A^{\prime} A^{\prime \prime} is a timelike geodesic which comes in from r=+oo*BB^(')B^('')r=+\infty \cdot B B^{\prime} B^{\prime \prime} is the null geodesic travelled by a photon that falls straight in (no summit; never at rest!). DD^(')D^('')D D^{\prime} D^{\prime \prime} is a spacelike radial geodesic. So is CC^(')C C^{\prime}, but E^(')E^('')E^{\prime} E^{\prime \prime} is timelike. Neither of the latter two ever succeed in crossing r=2Mr=2 M. (Unanswered questions about these geodesics will answer themselves in Figure 31.4, where the same world lines are charted in a "Kruskal-Szekeres diagram").
Described mathematically via equation (31.10), the geodesic F_("inverse ")^('')F_("inverse ")^(')FF^(')F^('')F_{\text {inverse }}^{\prime \prime} F_{\text {inverse }}^{\prime} F F^{\prime} F^{\prime \prime} starts with ejection at
(event F^('')F^{\prime \prime} in diagram).
spread them out into a line in a new ( r_("new "),t_("new ")r_{\text {new }}, t_{\text {new }} )-plane; and will squeeze the "line" (r=2M,t(r=2 M, t from -oo-\infty to +oo)+\infty) into a single point in the ( {:r_("new "),t_("new "))\left.r_{\text {new }}, t_{\text {new }}\right)-plane. One is the more prepared to accept this tentative conclusion and act on it because one has already seen (equation 31.8) that the region covering the (theta,phi)2(\theta, \phi) 2-sphere at r=2Mr=2 M, and extending from t=-oot=-\infty to t=+oot=+\infty, has zero proper volume. What timelier indication could one want that the "line" r=2M,-oo < t < oor=2 M,-\infty<t<\infty, is actually a point?
§31.4. SEVERAL WELL-BEHAVED COORDINATE SYSTEMS
The well-behaved coordinate system that is easiest to visualize is one in which the radially moving test particles of equations (31.10) remain always at rest ("comoving
Novikov coordinates:
(1) how constructed
(2) line element
coordinates"). Such coordinates were first used by Novikov (1963). Novikov attaches a specific value of his radial coordinate, R^(**)R^{*}, to each test particle as it emerges from the singularity of infinite tidal forces at r=0r=0, and insists that the particle carry that value of R^(**)R^{*} throughout its "cycloidal life"-up through r=2Mr=2 M to r=r_(max)r=r_{\max }, then back down through r=2Mr=2 M to r=0r=0. For definiteness, Novikov expresses the R^(**)R^{*} value for each particle in terms of the peak point on its trajectory by
As a time coordinate, Novikov uses proper time tau\tau of the test particles, normalized so tau=0\tau=0 at the peak of the orbit. Every particle in the swarm is ejected in such a manner that it arrives at the summit of its trajectory (r=r_(max),tau=0)\left(r=r_{\max }, \tau=0\right) at one and the same value of the Schwarzschild coordinate time; namely, at t=0t=0.
Simple though they may be conceptually, the Novikov coordinates are related to the original Schwarzschild coordinates by a very complicated transformation: (1) combine equations (31.10b) and (31.11) to obtain eta(tau,R^(**))\eta\left(\tau, R^{*}\right); (2) combine eta(tau,R^(**))\eta\left(\tau, R^{*}\right) with (31.10a) and (31.11) to obtain r(tau,R^(**))r\left(\tau, R^{*}\right); (3) combine eta(tau,R^(**))\eta\left(\tau, R^{*}\right) with (31.10c) and (31.11) to obtain t(tau,R^(**))t\left(\tau, R^{*}\right). The resulting coordinate transformation, when applied to the Schwarzschild metric (31.1), yields the line element
{:(31.12a)ds^(2)=-dtau^(2)+((R^(**2)+1)/(R^(**2)))((del r)/(delR^(**)))^(2)dR^(**2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-d \tau^{2}+\left(\frac{R^{* 2}+1}{R^{* 2}}\right)\left(\frac{\partial r}{\partial R^{*}}\right)^{2} d R^{* 2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{31.12a}
\end{equation*}
("Schwarzschild geometry in Novikov coordinates".) Here rr is no longer a radial coordinate; it is now a metric function r(tau,R^(**))r\left(\tau, R^{*}\right) given implicitly by
Figure 31.2 shows the locations of several key regions of Schwarzschild spacetime in this coordinate system. The existence of two distinct regions with r=0r=0 (singularities) and two distinct regions with r longrightarrow oor \longrightarrow \infty (asymptotically flat regions; recall that 4pir^(2)=4 \pi r^{2}= surface area!) will be discussed in §31.5\S 31.5§.
Figure 31.2.
The Novikov coordinate system for Schwarzschild spacetime (schematic). The dashed curves are curves of constant rr (recall: 4pir^(2)=4 \pi r^{2}= surface area about center of symmetry). The region shaded gray is not part of spacetime; it corresponds to r < 0r<0, a region that cannot be reached because of the singularity of spacetime at r=0r=0. Notice that the "line" (r=2M,-oo <(r=2 M,-\infty<t < +oot<+\infty ) of the Schwarzschild coordinate diagram (Figure 31.1) has been compressed into a point here, in accordance with the discussion at the end of §31.3\S 31.3§.
Although Novikov's coordinate system is very simple conceptually, the mathematical expressions for the metric components in it are rather unwieldy. Simpler, more usable expressions have been obtained in a different coordinate system ("Kruskal-Szekeres coordinates") by Kruskal (1960), and independently by Szekeres (1960).
Kruskal and Szekeres use a dimensionless radial coordinate uu and a dimensionless Kruskal-Szekeres coordinates time coordinate vv related to the Schwarzschild rr and tt by
{:(31.13b){:[u=(r//2M-1)^(1//2)e^(r//4M)cosh(t//4M)],[v=(r//2M-1)^(1//2)e^(r//4M)sinh(t//4M)]}" when "r > 2M",":}\left.\begin{array}{l}
u=(r / 2 M-1)^{1 / 2} e^{r / 4 M} \cosh (t / 4 M) \\
v=(r / 2 M-1)^{1 / 2} e^{r / 4 M} \sinh (t / 4 M) \tag{31.13b}
\end{array}\right\} \text { when } r>2 M,
(Motivation for introducing such coordinates is given in Box 31.2.) By making this change of coordinates in the Schwarzschild metric (31.1), one obtains the following line element:
{:(31.14a)ds^(2)=(32M^(3)//r)e^(-r//2M)(-dv^(2)+du^(2))+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=\left(32 M^{3} / r\right) e^{-r / 2 M}\left(-d v^{2}+d u^{2}\right)+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{31.14a}
\end{equation*}
("Schwarzschild geometry in Kruskal-Szekeres coordinates"). Here rr is to be regarded as a function of uu and vv defined implicitly by
Box 31.2 MOTIVATION FOR KRUSKAL-SZEKERES COORDINATES*
A. EDDINGTON-FINKELSTEIN COORDINATES
The motivation for the Kruskal-Szekeres system begins by introducing a different coordinate system, first devised by Eddington (1924) and rediscovered by Finkelstein (1958). Eddington and Finkelstein use as the foundation of their coordinate system, not freely falling particles as did Novikov, but freely falling photons. More particularly, they introduce coordinates widetilde(U)\widetilde{U} and widetilde(V)\widetilde{V}, which are labels for outgoing and ingoing, radial, null geodesics. The geodesics are given by
ds^(2)=0=-(1-2M//r)dt^(2)+(1-2M//r)^(-1)dr^(2).d s^{2}=0=-(1-2 M / r) d t^{2}+(1-2 M / r)^{-1} d r^{2} .
Equivalently, outgoing geodesics are given by widetilde(U)=\widetilde{U}= const, where
Ingoing Eddington-Finkelstein Coordinates-Adopt rr and widetilde(V)\widetilde{V} as coordinates in place of rr and tt
The Schwarzschild metric becomes.
{:(2)ds^(2)=-(1-2M//r)d widetilde(V)^(2)+2d widetilde(V)dr+r^(2)dOmega^(2):}\begin{equation*}
d s^{2}=-(1-2 M / r) d \widetilde{V}^{2}+2 d \widetilde{V} d r+r^{2} d \Omega^{2} \tag{2}
\end{equation*}
The radial light cone, ds^(2)=0d s^{2}=0, has one leg
{:(3a)d widetilde(V)//dr=0",":}\begin{equation*}
d \widetilde{V} / d r=0, \tag{3a}
\end{equation*}
and the other leg
{:(3b)(d( widetilde(V)))/(dr)=(2)/(1-2M//r):}\begin{equation*}
\frac{d \widetilde{V}}{d r}=\frac{2}{1-2 M / r} \tag{3b}
\end{equation*}
From this, and this alone, one can infer all features of the drawing.
Ingoing Eddington-Finkelstein coordinates (one rotational degree of freedom is suppressed; i.e., theta\theta is set equal to pi//2\pi / 2 ). Surfaces of constant widetilde(V)\widetilde{V}, being ingoing null surfaces, are plotted on a 45 -degree slant, just as they would be in flat spacetime. Equivalently, surfaces of constant
Outgoing Eddington-Finkelstein Coordinates-Adopt rr and widetilde(U)\widetilde{U} as coordinates in place of rr and tt
The Schwarzschild metric becomes
{:(4)ds^(2)=-(1-2M//r)d widetilde(U)^(2)-2d widetilde(U)dr+r^(2)dOmega^(2):}\begin{equation*}
d s^{2}=-(1-2 M / r) d \widetilde{U}^{2}-2 d \widetilde{U} d r+r^{2} d \Omega^{2} \tag{4}
\end{equation*}
Box 31.2 (continued)
The radial light cone, ds^(2)=0d s^{2}=0, has one leg
{:(5a)d widetilde(U)//dr=0",":}\begin{equation*}
d \widetilde{U} / d r=0, \tag{5a}
\end{equation*}
From this, and this alone, one can infer all features of the drawing.
Outgoing Eddington-Finkelstein coordinates (one rotational degree of freedom is suppressed). (Surfaces of constant widehat(U)\widehat{U}, being outgoing null surfaces, are plotted on a 45 -degree slant, just as they would be in flat spacetime.)
Notice that both Eddington-Finkelstein coordinate systems are better behaved at the gravitational radius than is the Schwarzschild coordinate system; but they are not fully well-behaved. The outgoing coordinates ( widetilde(U),r,theta,phi\widetilde{U}, r, \theta, \phi ) describe in a nonpathological manner the ejection of particles outward from r=0r=0 through r=2Mr=2 M; but their description of infall through r=2Mr=2 M has the same pathology as the description given by Schwarzschild coordinates (Figure 31.1). Similarly, the ingoing coordinates ( widetilde(V),r,theta,phi\widetilde{V}, r, \theta, \phi ) describe well the infall of a particle through r=2Mr=2 M, but they give a pathological description of outgoing trajectories. Moreover, the contrast between the two diagrams seems paradoxical: in one the gravitational radius is made up of world lines of outgoing photons; in the other it is made up of world lines of ingoing photons! To resolve the paradox, one must seek another, better-behaved coordinate system. [But note: because the ingoing Eddington-Finklestein coordinates describe infall so well, they are used extensively in discussions of gravitational collapse (Chapter 32) and black holes (Chapters 33 and 34).]
B. TRANSITION FROM EDDINGTON-FINKELSTEIN TO KRUSKAL-SZEKERES
Perhaps one would obtain a fully well-behaved coordinate system by dropping rr from view and using widetilde(U), widetilde(nu)\widetilde{U}, \widetilde{\nu}, as the two coordinates in the radial-time plane. The resulting coordinate system is related to Schwarzschild coordinates by [see equations (1)]
and the line element in terms of the new coordinates reads
{:(7)ds^(2)=-(1-2M//r)d widetilde(U)d widetilde(V)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-(1-2 M / r) d \widetilde{U} d \widetilde{V}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{7}
\end{equation*}
Contrary to one's hopes, this coordinate system is pathological at r=2Mr=2 M.
Second thoughts about the construction reveal the trouble: the surfaces widetilde(U)=\widetilde{U}= constant (outgoing null surfaces) used in constructing it are geometrically welldefined, as are the surfaces widetilde(V)=\widetilde{V}= constant (ingoing null surfaces); but the way of labeling them is not. Any relabeling, widetilde(u)=F( widetilde(U))\widetilde{u}=F(\widetilde{U}) and widetilde(v)=G( widetilde(V))\widetilde{v}=G(\widetilde{V}), will leave the surfaces unchanged physically. What one needs is a relabeling that will get rid of the singular factor 1-2M//r1-2 M / r in the line element (7). A successful relabeling is suggested by the equation
will remove the offending 1-2M//r1-2 M / r from the metric coefficients. In terms of these new coordinates, the line element reads
{:(10a)ds^(2)=-(32M^(3)//r)e^(-r//2M)d widetilde(v)d widetilde(u)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-\left(32 M^{3} / r\right) e^{-r / 2 M} d \widetilde{v} d \widetilde{u}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{10a}
\end{equation*}
Here rr is still defined by 4pir^(2)=4 \pi r^{2}= surface area, but it must be regarded as a function of widetilde(v)\widetilde{v} and widetilde(u)\widetilde{u} :
One can readily verify that this equation determines rr uniquely (recall: r > 0r>0 !) in terms of the product widetilde(u) widetilde(v)\widetilde{u} \widetilde{v} [details in Misner (1969a)].
The coordinates, widetilde(u), widetilde(v)\widetilde{u}, \widetilde{v}, which label the ingoing and outgoing null surfaces, are null coordinates; i.e.,
[see equation (10a)]. If one is not accustomed to working with null coordinates, it is helpful to replace widetilde(u)\widetilde{u} and widetilde(v)\widetilde{v} by spacelike and timelike coordinates, uu and vv (KruskalSzekeres coordinates!) defined by
{:(12)dv^(2)-du^(2)=d widetilde(v)d widetilde(u):}\begin{equation*}
d v^{2}-d u^{2}=d \widetilde{v} d \widetilde{u} \tag{12}
\end{equation*}
In terms of these coordinates, the line element has the Kruskal form (31.14), which is fully well-behaved at the gravitational radius.
Although the Kruskal-Szekeres line element is well behaved at r=2Mr=2 M, the transformation (11) from Schwarzschild to Kruskal-Szekeres is not; it becomes meaningless ( uu and vv "imaginary") when one moves from r > 2Mr>2 M to r < 2Mr<2 M. Of course, this is a manifestation of the pathologies of Schwarzschild coordinates. By trial and error, one readily finds a new transformation, to replace (11) at r < 2Mr<2 M, leading from Schwarzschild to Kruskal-Szekeres coordinates:
§31.5. RELATIONSHIP BETWEEN KRUSKAL-SZEKERES COORDINATES AND SCHWARZSCHILD COORDINATES
In the Kruskal-Szekeres coordinate system, the singularity r=0r=0 is located at v^(2)-u^(2)=1v^{2}-u^{2}=1. Thus there are actually two singularities, not one; both
{:(31.15)v=+(1+u^(2))^(1//2)" and "v=-(1+u^(2))^(1//2)" correspond to "r=0!:}\begin{equation*}
v=+\left(1+u^{2}\right)^{1 / 2} \text { and } v=-\left(1+u^{2}\right)^{1 / 2} \text { correspond to } r=0! \tag{31.15}
\end{equation*}
This is not the only surprise that lies hidden in the Kruskal-Szekeres line element (31.14). Notice also that r≫2Mr \gg 2 M (the region of spacetime far outside the gravitational radius) is given by u^(2)≫v^(2)u^{2} \gg v^{2}. Thus there are actually two exterior regions*; both
{:(31.16)u≫+|v|" and "u≪-|v|" correspond to "r≫2M!:}\begin{equation*}
u \gg+|v| \text { and } u \ll-|v| \text { correspond to } r \gg 2 M! \tag{31.16}
\end{equation*}
How can this be? When the geometry is charted in Schwarzschild coordinates, it contains one singularity and one exterior region; but when expressed in KruskalSzekeres coordinates, it shows two of each. The answer must be that the Schwarzschild coordinates cover only part of the spacetime manifold; they must be only a local coordinate patch on the full manifold. Somehow, by means of the coordinate transformation that leads to Kruskal-Szekeres coordinates, one has analytically extended the limited Schwarzschild solution for the metric to cover all (or more nearly all) of the manifold.
To understand this covering more clearly, transform back from Kruskal-Szekeres coordinates to Schwarzschild coordinates (see Figure 31.3). The transformation equations, as written down in (31.13) were valid only for the quadrants u > |v|u>|v| [equation (31.13a)] and v > |u|v>|u| [equation (31.13b)] of Kruskal coordinates. Denote these quadrants by the numerals I and II; and denote the other quadrants by III and IV (see Figure 31.3). In the other quadrants, one can also transform the KruskalSzekeres line element (31.14) into the Schwarzschild line element (31.1); but slightly different transformation equations are needed. One easily verifies that the following sets of transformations work:
Kruskal-Szekeres coordinates reveal that Schwarzschild spacetime has two " r=0r=0 singularities" and two " r longrightarrow oor \longrightarrow \infty exterior regions"
Transformation between Schwarzschild coordinates and Kruskal-Szekeres coordinates
Figure 31.3.
The transformation of the Schwarzschild vacuum geometry between Schwarzschild and Kruskal-Szekeres coordinates. Two Schwarzschild coordinate patches I, II, and III, IV (illustrated in the upper and lower portions of Figure 31.5 ,a) are required to cover the complete Schwarzschild geometry, whereas a single Kruskal-Szekeres coordinate system suffices. The Schwarzschild geometry consists of four regions I, II, III, IV. Regions I and III represent two distinct, but identical, asymptotically flat universes in which r > 2Mr>2 M; while regions II and IV are two identical, but time-reversed, regions in which physical singularities (r=0)(r=0) evolve. The transformation laws that relate the Schwarzschild and Kruskal-Szekeres coordinate systems to each other are given by equations (31.17) and (31.18). In the Kruskal-Szekeres u,vu, v-plane, curves of constant rr are hyperbolae with asymptotes u=+-vu= \pm v, while curves of constant tt are straight lines through the origin.
The inverse transformations are
{:[(31.18a)(r//2M-1)e^(r//2M)=u^(2)-v^(2)" in I, II, III, IV; "],[(31.18b)t={[4Mtanh^(-1)(v//u)" in I and III, "],[4Mtanh^(-1)(u//v)" in II and IV. "]:}]:}\begin{gather*}
(r / 2 M-1) e^{r / 2 M}=u^{2}-v^{2} \text { in I, II, III, IV; } \tag{31.18a}\\
t=\left\{\begin{array}{l}
4 M \tanh ^{-1}(v / u) \text { in I and III, } \\
4 M \tanh ^{-1}(u / v) \text { in II and IV. }
\end{array}\right. \tag{31.18b}
\end{gather*}
These coordinate transformations are exhibited graphically in Figure 31.3. Notice that two Schwarzschild coordinate patches, I, II, and III, IV, are required to cover the entire Schwarzschild geometry; but a single Kruskal coordinate system suffices. Schwarzschild patch I, II, is divided into two regions-region I, which is outside the gravitational radius (r > 2M)(r>2 M), and region II, which is inside the gravitational radius ( r < 2Mr<2 M ). Similarly, Schwarzschild patch III, IV, consists of an exterior region (III) and an interior region (IV).
Two Schwarzschild coordinate patches are required to cover all of spacetime
Figure 31.4.
(a) Typical radial timelike (A,E,F)(A, E, F), lightlike (B)(B), and spacelike (C,D)(C, D) geodesics of the Schwarzschild geometry, as seen in the Schwarzschild coordinate system (schematic only). This is a reproduction of Figure 31.1.
(b) The same geodesics, as seen in the Kruskal-Szekeres coordinate system, and as extended either to infinite length or to the singularity of infinite curvature at r=0r=0 (schematic only).
Equations (31.18) reveal that the regions of constant rr (constant surface area) are hyperbolae with asymptotes u=+-vu= \pm v in the Kruskal-Szekeres diagram, and that regions of constant tt are straight lines through the origin.
Several radial geodesics of the complete Schwarzschild geometry are depicted in the Kruskal-Szekeres coordinate system in Figure 31.4. Notice how much more reasonable the geodesic curves look in Kruskal-Szekeres coordinates than in Schwarzschild coordinates. Notice also that radial, lightlike geodesics (paths of radial light rays) are 45-degree lines in the Kruskal-Szekeres coordinate system. This can be seen from the Kruskal-Szekeres line element (31.14), for which du=+-dvd u= \pm d v guarantees ds=0d s=0. Because of this 45 -degree property, the radial light cone in a Kruskal-Szekeres diagram has the same form as in the space-time diagram of special relativity. Any radial curve that points "generally upward" (i.e., makes an angle of less than 45 degrees with the vertical, vv, axis) is timelike; and curves that point "generally outward" are spacelike. This property enables a Kruskal-Szekeres diagram to exhibit easily the causality relation between one event in spacetime and another (see exercises 31.2 to 31.4 ).
Properties of the Kruskal-Szekeres coordinate system
EXERCISES
Exercise 31.2. NONRADIAL LIGHT CONES
Show that the world line of a photon traveling nonradially makes an angle less than 45 degrees with the vertical vv-axis of a Kruskal-Szekeres coordinate diagram. From this, infer that particles with finite rest mass, traveling nonradially or radially, must always move "generally upward" (angle less than 45 degrees with vertical vv-axis).
Exercise 31.3. THE CRACK OF DOOM
Use a Kruskal diagram to show the following.
(a) If a man allows himself to fall through the gravitational radius r=2Mr=2 M, there is no way whatsoever for him to avoid hitting (and being killed in) the singularity at r=0r=0.
(b) Once a man has fallen inward through r=2Mr=2 M, there is no way whatsoever that he can send messages out to his friends at r > 2Mr>2 M, but he can still receive messages from them (e.g., by radio waves, or laser beam, or infalling "CARE packages").
Exercise 31.4. HOW LONG TO LIVE?
Show that once a man falling inward reaches the gravitational radius, no matter what he does subsequently (no matter in what directions, how long, and how hard he blasts his rocket engines), he will be pulled into the singularity and killed in a proper time of
[Hint: The trajectory of longest proper time lapse must be a geodesic. Use the mathematical tools of Chapter 25 to show that the geodesic of longest proper time lapse between r=2Mr=2 \mathrm{M} and r=0r=0 is the radial geodesic (31.10a), with r_(max)=2Mr_{\max }=2 M, for which the time lapse is pi M\pi M.]
Exercise 31.5. EDDINGTON-FINKELSTEIN AND KRUSKAL-SZEKERES COMPARED
Use coordinate diagrams to compare the ingoing and outgoing Eddington-Finkelstein coordinates of Box 31.2 with the Kruskal-Szekeres coordinates. Pattern the comparison after that between Schwarzschild and Kruskal-Szekeres in Figures 31.3 and 31.4.
Exercise 31.6. ANOTHER COORDINATE SYSTEM
Construct a coordinate diagram for the widetilde(U), widetilde(V),theta,phi\widetilde{U}, \widetilde{V}, \theta, \phi coordinate system of Box 31.2 [equations (6) and (7)]. Show such features as (1) the relationship to Schwarzschild and to KruskalSzekeres coordinates; (2) the location of r=2Mr=2 M; and (3) radial geodesics.
§31.6. DYNAMICS OF THE SCHWARZSCHILD GEOMETRY
What does the Schwarzschild geometry look like? This question is most readily answered by means of embedding diagrams analogous to those for an equilibrium star ( $23.8\$ 23.8; Figure 23.1; and end of Box 23.2) and for Friedmann universes of positive and negative spatial curvature [equations (27.23) and (27.24) and Box 27.2].
Examine, first, the geometry of the spacelike hypersurface v=0v=0, which extends from u=+oo(r=oo)u=+\infty(r=\infty) into u=0(r=2M)u=0(r=2 M) and then out to u=-oo(r=oo)u=-\infty(r=\infty). In Schwarzschild coordinates this surface is a slice of constant time, t=0t=0 [see equation (31.18b)]; it is precisely the surface for which an embedding diagram was calculated in equation (23.34b). The embedded surface, with one degree of rotational freedom suppressed, is described by the paraboloid of revolution
{:(31.20) bar(r)=2M+ bar(z)^(2)//8M:}\begin{equation*}
\bar{r}=2 M+\bar{z}^{2} / 8 M \tag{31.20}
\end{equation*}
Figure 31.5.
(a) The Schwarzschild space geometry at the "moment of time" t=v=0t=v=0, with one degree of rotational freedom suppressed (theta=pi//2)(\theta=\pi / 2). To restore that rotational freedom and obtain the full Schwarzschild 3-geometry, one mentally replaces the circles of constant bar(r)=( bar(x)^(2)+ bar(y)^(2))^(1//2)\bar{r}=\left(\bar{x}^{2}+\bar{y}^{2}\right)^{1 / 2} with spherical surfaces of area 4pi bar(r)^(2)4 \pi \bar{r}^{2}. Note that the resultant 3 -geometry becomes flat (Euclidean) far from the throat of the bridge in both directions (both "universes").
(b) An embedding of the Schwarzschild space geometry at "time" t=v=0t=v=0, which is geometrically identical to the embedding (a), but which is topologically different. Einstein's field equations fix the local geometry of spacetime, but they do not fix its topology; see the discussion at end of Box 27.2 . Here the Schwarzschild "wormhole" connects two distant regions of a single, asymptotically flat universe. For a discussion of issues of causality associated with this choice of topology, see Fuller and Wheeler (1962).
in the flat Euclidean space with metric
{:(31.21)dsigma^(2)=d bar(r)^(2)+d bar(z)^(2)+ bar(r)^(2)d bar(phi)^(2):}\begin{equation*}
d \sigma^{2}=d \bar{r}^{2}+d \bar{z}^{2}+\bar{r}^{2} d \bar{\phi}^{2} \tag{31.21}
\end{equation*}
(See Figure 31.5.)
Notice from the embedding diagram of Figure 31.5,a, that the Schwarzschild
The 3-surface v=t=0v=t=0 is a "wormhole" connecting two asymptotically flat universes, or two different regions of one universe
Schwarzschild geometry is dynamic in regions r < 2Mr<2 M
Time evolution of the wormhole: creation; expansion; recontraction; and pinch-off
Communication through the wormhole is impossible: it pinches off too fast
geometry on the spacelike hypersurface t=t= const consists of a bridge or "wormhole" connecting two distinct, but identical, asymptotically flat universes. This bridge is sometimes called the "Einstein-Rosen bridge" and sometimes the "Schwarzschild throat" or the "Schwarzschild wormhole." If one so wishes, one can change the topology of the Schwarzschild geometry by connecting the two asymptotically flat universes together in a region distant from the Schwarzschild throat [Fuller and Wheeler (1962); Fig. 31.5b]. The single, unique universe then becomes multiply connected, with the Schwarzschild throat providing one spacelike path from point aa to point B\mathscr{B}, and the nearly flat universe providing another. For concreteness, focus attention on the interpretation of the Schwarzschild geometry, not in terms of Wheeler's multiply connected single universe, but rather in terms of the EinsteinRosen double universe of Figure 31.5,a.
One is usually accustomed to think of the Schwarzschild geometry as static. However, the static "time translations," t longrightarrow t+Delta tt \longrightarrow t+\Delta t, which leave the Schwarzschild geometry unchanged, are time translations in the strict sense of the words only in regions I and III of the Schwarzschild geometry. In regions II and IV, t longrightarrow t+Delta tt \longrightarrow t+\Delta t is a spacelike motion, not a timelike motion (see Fig. 31.3). Consequently, a spacelike hypersurface, such as the surface t=t= const of Figure 31.5,a, which extends from region I through u=v=0u=v=0 into region III, is not static. As this spacelike hypersurface is pushed forward in time (in the +v+v direction of the Kruskal diagram), it enters region II, and its geometry begins to change.
In order to examine the time-development of the Schwarzschild geometry, one needs a sequence of embedding diagrams, each corresponding to the geometry of a spacelike hypersurface to the future of the preceding one. But how are the hypersurfaces to be chosen? In Newtonian theory or special relativity, one chooses hypersurfaces of constant time. But in dynamic regions of curved spacetime, no naturally preferred time coordinate exists. This situation forces one to make a totally arbitrary choice of hypersurfaces to use in visualizing the time-development of geometry, and to keep in mind how very arbitrary that choice was.
Figure 31.6 uses two very different choices of hypersurfaces to depict the timedevelopment of the Schwarzschild geometry. (Still other choices are shown in Figure 21.4.) Notice that the precise geometry of the evolving bridge depends on the arbitrary choice of spacelike hypersurfaces, but that the qualitative nature of the evolution is independent of the choice of hypersurfaces. Qualitatively speaking, the two asymptotically flat universes begin disconnected, with each one containing a singularity of infinite curvature (r=0)(r=0). As the two universes evolve in time, their singularities join each other and form a nonsingular bridge. The bridge enlarges, until it reaches a maximum radius at the throat of r=2Mr=2 M (maximum circumference of 4pi M4 \pi M; maximum surface area of 16 piM^(2)16 \pi M^{2} ). It then contracts and pinches off, leaving the two universes disconnected and containing singularities (r=0)(r=0) once again. The formation, expansion, and collapse of the bridge occur so rapidly, that no particle or light ray can pass across the bridge from the faraway region of the one universe to the faraway region of the other without getting caught and crushed in the throat as it pinches off. (To verify this, examine the Kruskal-Szekeres diagram of Figure 31.3, where radial light rays move along 45 -degree lines.)
Spacelike slices
History A-B-C-D-E-F-GA-B-C-D-E-F-G
History A-W-X-D-Y-Z-GA-W-X-D-Y-Z-G
Figure 31.6.
Dynamical evolution of the Einstein-Rosen bridge of the vacuum Schwarzschild geometry (schematic). Shown here are two sequences of embedding diagrams corresponding to two different ways of viewing the evolution of the bridge-History A-B-C-D-E-F-GA-B-C-D-E-F-G, and History A-W-X-D-Y-Z-GA-W-X-D-Y-Z-G. The embedding diagrams are skeletonized in that each diagram must be rotated about the appropriate vertical axis in order to become two-dimensional surfaces analogous to Figure 31.5 , a. [Notice that the hypersurfaces of which embedding diagrams are given intersect the singularity only tangentially. Hypersurfaces that intersect the singularity at a finite angle in the u,vu, v-plane are not shown because they cannot be embedded in a Euclidean space. Instead, a Minkowski space (indefinite metric) must be used, at least near r=0r=0. For an example of an embedding in Minkowski space, see the discussion of a universe with constant negative spatial curvature in equations (27.23) and (27.24) and Box 27.2C.] Figure 21.4 exhibits embedding diagrams for other spacelike slices in the Schwarzschild geometry.
From the Kruskal-Szekeres diagram and the 45 -degree nature of its radial light rays, one sees that any particle that ever finds itself in region IV of spacetime must have been "created" in the earlier singularity; and any particle that ever falls into region II is doomed to be crushed in the later singularity. Only particles that stay forever in one of the asymptotically flat universes I or III, outside the gravitational radius (r > 2M)(r>2 M), are forever safe from the singularities.
Some investigators, disturbed by the singularities at r=0r=0 or by the "double-universe" nature of the Schwarzschild geometry, have proposed modifications of its topology. One proposal is that the earlier and later singularities be identified with each other, so that a particle which falls into the singularity of region II, instead of being destroyed, will suddenly reemerge, being ejected, from the singularity of region IV. One cannot overstate the objections to this viewpoint: the region r=0r=0 is a physical singularity of infinite tidal gravitation forces and infinite Riemann curvature. Any particle that falls into that singularity must be destroyed by those
Creation and destruction in the singularities
Nonviable proposals for modifying the topology of Schwarzschild spacetime
forces. Any attempt to extrapolate its fate through the singularity using Einstein's field equations must fail; the equations lose their predictive power in the face of infinite curvature. Consequently, to postulate that the particle reemerges from the earlier singularity is to make up an ad hoc mathematical rule, one unrelated to physics. It is conceivable, but few believe it true, that any object of finite mass will modify the geometry of the singularity as it approaches r=0r=0 to such an extent that it can pass through and reemerge. However, whether such a speculation is correct must be answered not by ad hoc rules, but by concrete, difficult computations within the framework of general relativity theory (see Chapter 34).
A second proposal for modifying the topology of the Schwarzschild geometry is this: one should avoid the existence of two different asymptotically flat universes by identifying each point (v,u,theta,phi)(v, u, \theta, \phi) with its opposite point (-v,-u,theta,phi)(-v,-u, \theta, \phi) in the Kruskal-Szekeres coordinate system. Two objections to this proposal are: (1) it produces a sort of "conical" singularity (absence of local Lorentz frames) at (v,u)=(v, u)=(0,0)(0,0), i.e., at the neck of the bridge at its moment of maximum expansion; and (2) it leads to causality violations in which a man can meet himself going backward in time.
One good way for the reader to become conversant with the basic features of the Schwarzschild geometry is to reread $$31.1-31.4\$ \$ 31.1-31.4 carefully, reinterpreting everything said there in terms of the Kruskal-Szekeres diagram.
EXERCISES
Exercise 31.7. SCHWARZSCHILD METRIC IN ISOTROPIC COORDINATES
(a) Show that, rewritten in the isotropic coordinates of Exercise 23.1, the Schwarzschild metric reads
{:(31.22)ds^(2)=-((1-M//2( bar(r)))/(1+M//2( bar(r))))^(2)dt^(2)+(1+(M)/(2( bar(r))))^(4)[d bar(r)^(2)+ bar(r)^(2)(dtheta^(2)+sin^(2)theta dphi^(2))];:}\begin{equation*}
d s^{2}=-\left(\frac{1-M / 2 \bar{r}}{1+M / 2 \bar{r}}\right)^{2} d t^{2}+\left(1+\frac{M}{2 \bar{r}}\right)^{4}\left[d \bar{r}^{2}+\bar{r}^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)\right] ; \tag{31.22}
\end{equation*}
between the two radial coordinates.
(b) Which regions of spacetime (I, II, III, IV; see Figure 31.3) are covered by the isotropic coordinate patch, and which are not?
(c) Calculate and construct an embedding diagram for the spacelike hypersurface t=0t=0, 0 < bar(r) < oo0<\bar{r}<\infty.
(d) Find a coordinate transformation that interchanges the region near bar(r)=0\bar{r}=0 with the region near bar(r)=oo\bar{r}=\infty, while leaving the metric coefficients in their original form.
Exercise 31.8. REISSNER-NORDSTR O/\varnothing M GEOMETRY
(a) Solve the Einstein field equations for a spherically symmetric, static gravitational field
ds^(2)=-e^(2phi(r))dt^(2)+e^(2A(r))dr^(2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)),d s^{2}=-e^{2 \phi(r)} d t^{2}+e^{2 A(r)} d r^{2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right),
with no matter present, but with a radial electric field B=0,E=f(r)e_( hat(r))\boldsymbol{B}=0, \boldsymbol{E}=f(r) \boldsymbol{e}_{\hat{r}} in the static orthonormal frame
Use as a source in the Einstein field equations the stress-energy of the electric field. [Answer:
{:[(31.24a)E=(Q//r^(2))e_( hat(r))],[(31.24b)ds^(2)=-(1-(2M)/(r)+(Q^(2))/(r^(2)))dt^(2)+(1-(2M)/(r)+(Q^(2))/(r^(2)))^(-1)dr^(2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))]:}\begin{gather*}
\boldsymbol{E}=\left(Q / r^{2}\right) \boldsymbol{e}_{\hat{r}} \tag{31.24a}\\
d s^{2}=-\left(1-\frac{2 M}{r}+\frac{Q^{2}}{r^{2}}\right) d t^{2}+\left(1-\frac{2 M}{r}+\frac{Q^{2}}{r^{2}}\right)^{-1} d r^{2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{31.24b}
\end{gather*}
This is called the "Reissner (1916)-Nordstrom (1918) metric".]
(b) Show that the constant QQ is the total charge as measured by a distant observer (r >= 2M(r \geqslant 2 M and r≫Qr \gg Q ), who uses a Gaussian flux integral, or who studies the coulomb-force-dominated orbits of test charges with charge-to-mass ratio e//mu >= M//Qe / \mu \geqslant M / Q. What is the charge-to-mass ratio, in dimensionless units, for an electron? Show that the constant MM is the total mass as measured by a distant observer using the Keplerian orbits of electrically neutral particles.
(c) Show that for Q > MQ>M, the Reissner-Nordstrom coordinate system is well-behaved from r=oor=\infty down to r=0r=0, where there is a physical singularity and infinite tidal forces.
(d) Explore the nature of the spacetime geometry for Q < MQ<M, using all the techniques of this chapter (coordinate transformations, Kruskal-like coordinates, studies of particle orbits, embedding diagrams, etc.).
[Solution: see Graves and Brill (1960); also Fig. 34.4 of this book.]
(e) Similarly explore the spacetime geometry for Q=MQ=M. [Solution: see Carter (1966b).]
(f) For the case of a large ratio of charge to mass [Q > M[Q>M as in part (c)], show that the region near r=0r=0 is unphysical. More precisely, show that any spherically symmetric distribution of charged stressed matter that gives rise to the fields (31.24) outside its boundary must modify these fields for r < r_(0)=Q^(2)//2Mr<r_{0}=Q^{2} / 2 M. [Hint: Study the quantity m(r)m(r) defined in equations (23.18) and (32.22h), noting its values deduced from equation (31.24), on the one hand, and from the appropriate Einstein equation within the matter distribution, on the other hand. See Figure 26 of Misner (1969a) for a similar argument.]
cmorest 32
GRAVITATIONAL COLLAPSE
Now, here, you see, it takes all the running you can do, to keep in the same place. If you want to get somewhere else, you must run twice as fast as that.
The roles and relevance of the Schwarzschild geometry
§32.1. RELEVANCE OF SCHWARZSCHILD GEOMETRY
The story that unfolded in the preceeding chapter was fantastic! One began with the innocuous looking Schwarzschild line element
{:(32.1)ds^(2)=-(1-(2M)/(r))dt^(2)+(dr^(2))/(1-2M//r)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-\left(1-\frac{2 M}{r}\right) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{32.1}
\end{equation*}
which was derived originally as the external field of a static star. One asked what happens if the star is absent; i.e., one probed the nature of the Schwarzschild geometry when no star is present to generate it. One might have expected the geometry to be that of a point mass sitting at r=0r=0. But it was not. It turned out to represent a "wormhole" connecting two asymptotically flat universes. Moreover, the wormhole was dynamic. It was created by the "joining together" of two " r=0r=0 " singularities, one in each universe; it expanded to a maximum circumference of 4pi M4 \pi M; it then recontracted and pinched off, leaving the two universes disconnected once again, each with its own " r=0r=0 " singularity.
As a solution to Einstein's field equations, this expanding and recontracting wormhole must be taken seriously. It is an exact solution; and it is one of the simplest of all exact solutions. But there is no reason whatsoever to believe that such wormholes exist in the real universe! They can exist only if the expanding universe, ∼10 xx10^(9)\sim 10 \times 10^{9} years ago, was "born" with the necessary initial conditions-with " r=0r=0 "
Schwarzschild singularities ready and waiting to blossom forth into wormholes. There is no reason at all to believe in such pathological initial conditions!
Why, then, was so much time and effort spent in Chapter 31 on understanding the Schwarzschild geometry? (1) Because it illustrates clearly the highly nonEuclidean character of spacetime geometry when gravity becomes strong; (2) because it illustrates many of the techniques one can use to analyze strong gravitational fields; and most importantly (3) because, when appropriately truncated, it is the spacetime geometry of a black hole and of a collapsing star-as well as of a wormhole.
This chapter explores the role of the Schwarzschild geometry in gravitational collapse; the next chapter explores its role in black-hole physics.
§32.2. BIRKHOFF'S THEOREM
That the Schwarzschild geometry is relevant to gravitational collapse follows from Birkhoff's (1923) theorem: Let the geometry of a given region of spacetime (1) be spherically symmetric, and (2) be a solution to the Einstein field equations in vacuum. Then that geometry is necessarily a piece of the Schwarzschild geometry. The external field of any electrically neutral, spherical star satisfies the conditions of Birkhoff's theorem, whether the star is static, vibrating, or collapsing. Therefore the external field must be a piece of the Schwarzschild geometry.
Birkhoff's theorem is easily understood on physical grounds. Consider an equilibrium configuration that is unstable against gravitational collapse and that, like all equilibrium configurations (see §23.6\S 23.6§ ), has the Schwarzschild geometry as its external gravitational field. Perturb this equilibrium configuration in a spherically symmetric way, so that it begins to collapse radially. The perturbation and subsequent collapse cannot affect the external gravitational field so long as exact spherical symmetry is maintained. Just as Maxwell's laws prohibit monopole electromagnetic waves, so Einstein's laws prohibit monopole gravitational waves. There is no possible way for any gravitational influence of the radial collapse to propagate outward.
Not only is Birkhoff's theorem easy to understand, but it is also fairly easy to prove. Consider a spherical region of spacetime. Spherical symmetry alone is sufficient to guarantee that conditions (i), (ii), and (iii) of Box 23.3 are satisfied, and thus to guarantee that one can introduce Schwarzschild coordinates
{:[ds^(2)=-e^(2phi)dt^(2)+e^(2Lambda)dr^(2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))","],[(32.2)Phi=Phi(t","r)","" and "Lambda=Lambda(t","r).]:}\begin{gather*}
d s^{2}=-e^{2 \phi} d t^{2}+e^{2 \Lambda} d r^{2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right), \\
\Phi=\Phi(t, r), \text { and } \Lambda=\Lambda(t, r) . \tag{32.2}
\end{gather*}
[See Box 23.3 for proof; and notice that: (1) for generality one must allow g_(tt)=-e^(2Phi)g_{t t}=-e^{2 \Phi} and g_(rr)=e^(2Lambda)g_{r r}=e^{2 \Lambda} to be positive or negative (no constraint on sign!); (2) at events where the gradient of the "circumference function" rr is zero or null, Schwarzschild coordinates cannot be introduced. The special case (grad r)^(2)=0(\boldsymbol{\nabla} r)^{2}=0 is treated in exercise 32.1.]
The uniqueness of the Schwarzschild geometry: Birkhoff's theorem
The physics underlying Birkhoff's theorem
Proof of Birkhoff's theorem
Impose Einstein's vacuum field equation on the metric (32.2), using the orthonormal components of the Einstein tensor as derived in exercise 14.16:
Equation (32.3b) guarantees that Lambda\Lambda is a function of rr only, and equation (32.3a) then guarantees that Lambda\Lambda has the same form as for the Schwarzschild metric:
{:(32.4a)Lambda=-(1)/(2)ln |1-2M//r|:}\begin{equation*}
\Lambda=-\frac{1}{2} \ln |1-2 M / r| \tag{32.4a}
\end{equation*}
Equations (32.3c,d) then become two equivalent equations for Phi(t,r)\Phi(t, r)-equivalent by virtue of the Bianchi identity, grad*G=0\boldsymbol{\nabla} \cdot \boldsymbol{G}=0-whose solution is
{:(32.4b)Phi=(1)/(2)ln |1-2M//r|+f(t):}\begin{equation*}
\Phi=\frac{1}{2} \ln |1-2 M / r|+f(t) \tag{32.4b}
\end{equation*}
Here ff is an arbitrary function. Put expressions (32.4) into the line element (32.2); thereby obtain
ds^(2)=-e^(2f(t))(1-(2M)/(r))dt^(2)+(dr^(2))/(1-2M//r)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))d s^{2}=-e^{2 f(t)}\left(1-\frac{2 M}{r}\right) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)
Then redefine the time coordinate
t_(new)=inte^(f(t))dtt_{\mathrm{new}}=\int e^{f(t)} d t
and thereby bring the line element into the Schwarzschild form
ds^(2)=-(1-(2M)/(r))dt^(2)+(dr^(2))/(1-2M//r)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))d s^{2}=-\left(1-\frac{2 M}{r}\right) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)
Conclusion: When the spacetime surrounding any object has spherical symmetry and is free of charge, mass, and all fields other than gravity, then one can introduce coordinates in which the metric is that of Schwarzschild. Conclusion restated in coordinate-free language: the geometry of any spherically symmetric vacuum region of spacetime is a piece of the Schwarzschild geometry (Birkhoff's theorem). Q.e.D.
EXERCISE
Exercise 32.1. UNIQUENESS OF REISSNER-NORDSTR phi\phi M GEOMETRY [Track 2]
Prove the following generalization of Birkhoff's theorem. Let the geometry of a given region of spacetime (1) be spherically symmetric, and (2) be a solution to the Einstein field equations
with an electromagnetic field as source. Then that geometry is necessarily a piece of the Reissner-Nordstrom geometry [equation (31.24b)] with electric and magnetic fields, as measured in the standard static orthonormal frames
[Hints: (1) First consider regions of spacetime in which (grad r)^(2)!=0(\boldsymbol{\nabla} r)^{2} \neq 0, using the same methods as the text uses for Birkhoff's theorem. The result is the Reissner-Nordstrom solution. (2) Any region of dimensionality less than four, in which (grad r)^(2)=0(\boldsymbol{\nabla} r)^{2}=0 (e.g., the Schwarzschild radius), can be treated as the join between four-dimensional regions with (grad r)^(2)!=0(\boldsymbol{\nabla} r)^{2} \neq 0. Moreover, the geometry of such a region is determined uniquely by the geometry of the adjoining four-dimensional regions ("junction conditions"; §21.13). Since the adjoining regions are necessarily Reissner-Nordstrøm (step 1), then so are such "sandwiched" regions. (3) Next consider four-dimensional regions in which grad r=dr\boldsymbol{\nabla} r=\boldsymbol{d} r is null and nonzero. Show that in such regions there exist coordinate systems with
ds^(2)=-2Psi drdt+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))d s^{2}=-2 \Psi d r d t+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)
where Psi=Psi(r,t)\Psi=\Psi(r, t). Show further that the Ricci tensor for this line element has an orthonormalized component
These quantities, R_( hat(theta) hat(theta))R_{\hat{\theta} \hat{\theta}} and 8piT_( hat(theta) hat(theta))8 \pi T_{\hat{\theta} \hat{\theta}}, must be equal (Einstein's field equation) but cannot be because of their different rr-dependence. Thus, an electromagnetic field cannot generate regions with dr!=0,dr*dr=0\boldsymbol{d} r \neq 0, \boldsymbol{d} r \cdot \boldsymbol{d} r=0. (4) Finally, consider four-dimensional regions in which dr=0\boldsymbol{d} r=0. Denote the constant value of rr by aa, and show that any event can be chosen as the origin of a locally well-behaved coordinate system with
{:[ds^(2)=a^(2)(-d tilde(tau)^(2)+e^(2lambda)dz^(2)+dtheta^(2)+sin^(2)theta dphi^(2))],[lambda=lambda( widetilde(tau)","z)","quad lambda( widetilde(tau)=0","z)=0","quadlambda^(˙)( tilde(tau)=0","z)=0.]:}\begin{gathered}
d s^{2}=a^{2}\left(-d \tilde{\tau}^{2}+e^{2 \lambda} d z^{2}+d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \\
\lambda=\lambda(\widetilde{\tau}, z), \quad \lambda(\widetilde{\tau}=0, z)=0, \quad \dot{\lambda}(\tilde{\tau}=0, z)=0 .
\end{gathered}
[Novikov-type coordinate system; see §31.4.] Show that, in the associated orthonormal frame, spherical symmetry demands
and that the Einstein field equations then require Q=aQ=a and e^(lambda)=cos tilde(tau)e^{\lambda}=\cos \tilde{\tau}, so that
ds^(2)=Q^(2)(-d tilde(tau)^(2)+cos^(2)( widetilde(tau))dz^(2)+dtheta^(2)+sin^(2)theta dphi^(2))d s^{2}=Q^{2}\left(-d \tilde{\tau}^{2}+\cos ^{2} \widetilde{\tau} d z^{2}+d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)
(5) This solution of the field equations [sometimes called the "Bertotti(1959)-Robinson (1959a) Electromagnetic Universe," and explored in this coordinate system by Lindquist (1960)] is actually the throat of the Reissner-Nordstrom solution for the special case Q=MQ=M. Verify this claim by performing the following coordinate transformation on the Reissner-Nordstrom throat region [equation (31.24b) with Q=MQ=M and |r-Q|≪Q|r-Q| \ll Q ]:
(6) Thus, each possible case leads either to no solution at all, or to a segment of the Reissner-Nordstrom geometry. Q.e.D.] Note: The missing case, (grad r)^(2)=0(\boldsymbol{\nabla} r)^{2}=0, in the text's proof of Birkhoff's theorem, is resolved by noting that, for Q=0Q=0, steps (3) and (4) above lead to no solutions at all. We thank G. F. R. Ellis for pointing out the omission of the case (grad r)^(2)=0(\nabla r)^{2}=0 in the preliminary version of this book.
Gravitational collapse analyzed by examining the star's exterior, Schwarzschild geometry
The gravitational radius as a point of no return, and the "crushing" at r=0r=0
§32.3. EXTERIOR GEOMETRY OF A COLLAPSING STAR
Consider a star that is momentarily static, but will subsequently begin to collapse. Its space geometry at the initial moment of Schwarzschild coordinate time, t=0t=0, has two parts: in the exterior, vacuum region (r > R > 2M)(r>R>2 M), it is the Schwarzschild geometry (Birkhoff's theorem!); but in the star's interior, it is some other, totally different geometry. Whatever the interior geometry may be, it has an embedding diagram at time t=0t=0 which is qualitatively like that of Figure 23.1. (For discussion and proof of this, see §23.8.) Notice that the star's space geometry is obtained by discarding the lower universe of the full Schwarzschild geometry (Figure 31.5,a), and replacing it with a smooth "bowl" on which the matter of the star is contained.
To follow the subsequent collapse of this star in the Schwarzschild coordinate system, or in the Kruskal-Szekeres coordinate system, or in an ingoing EddingtonFinkelstein coordinate system, one can similarly discard that part of the coordinate diagram which lies inside the star's surface, and keep only the exterior Schwarzschild region. (See Figure 32.1.) In place of the discarded interior Schwarzschild region, one must introduce some other coordinate system, line element, and diagram that correctly describe the interior of the collapsing star.
From truncated coordinate diagrams (such as Figures 32.1,a,b,c), one can readily discover and understand the various peculiar features of collapse through the gravitational radius.
(1) No matter how stiff may be the matter of which a (spherical) star is made, once its surface has collapsed within the gravitational radius, the star will continue to collapse until its surface gets crushed in the singularity at r=0r=0. This one discovers by recalling that the star's surface cannot move faster than the speed of light, so its world line must always make an angle of less than 45 degrees with the vv-axis of the Kruskal-Szekeres diagram.
(2) No signal (e.g., photon) emitted from the star's surface after it collapses inside the gravitational radius can ever escape to an external observer. Rather, all signals emitted from inside the gravitational radius get caught and destroyed by the collapse of the surrounding geometry into the singularity at r=0r=0 as space "pinches off" around the star.
(3) Consequently, an external observer can never see the star after it passes the gravitational radius; and he can never see the singularity that terminates its col-lapse-unless he chooses to fall through the gravitational radius himself and pay the price of death for the knowledge gained.
Does this mean that the collapsing star instantaneously and completely disappears from external view as it reaches the gravitational radius? No, not according to the analysis depicted in Figure 31.1,c: Place an astrophysicist on the surface of a collapsing star, and have him send a series of uniformly spaced signals to a distant astronomer, at rest at r≫2Mr \gg 2 M, to inform him of the progress of the collapse. These signals propagate along null lines in the spacetime diagram of Figure 31.1c. The signals originate on the world line of the stellar surface, and they are received by the distant astronomer when they intersect his world line, r=r= constant ≫M\gg M. As the star collapses closer and closer to its gravitational radius, R=2MR=2 M, the signals, which are sent at equally spaced intervals according to the astrophysicist's clock, are received by the astronomer at more and more widely spaced intervals. The astronomer does not receive a signal emitted just before the gravitational radius is reached until after an infinite amount of time has elapsed; and he never receives signals emitted after the gravitational radius has been passed. Those signals, like the astrophysicist who sends them, after brief runs get caught and destroyed by the collapsing geometry in the singularity, at r=0r=0. It is not only the star that collapses. The geometry around the star collapses.
Hence, to the distant astronomer, the collapsing star appears to slow down as it approaches its gravitational radius: light from the star becomes more and more red-shifted. Clocks on the star appear to run more and more slowly. It takes an infinite time for the star to reach its gravitational radius; and, as seen by the distant astronomer, the star never gets beyond there.
The optical appearance of a collapsing star was first analyzed mathematically, giving main attention to radially propagating photons, by J. R. Oppenheimer and H. Snyder (1939). More recently a number of workers have reexamined the problem [see, e.g., Podurets (1964), Ames and Thorne (1968) and Jaffe (1969)]. The most important quantitative results of these studies are as follows. In the late stages of collapse, when the distant astronomer sees the star to be very near its gravitational radius, he observes its total luminosity to decay exponentially in time
However, the light from the star is dominated in these late stages, not by photons flying along radial trajectories from near the gravitational radius, but by photons that were deposited by the star in unstable circular orbits as its surface passed through r=3Mr=3 M (see $25.6\$ 25.6 and Box 25.7 ). As time passes, these photons gradually leak out the diffuse spherical shell of trapped photons at r=3Mr=3 M and fly off to the distant observer, who measures them to have redshift z~~2z \approx 2. Consequently, in the late stages of collapse the star's spectral lines are broadened enormously, but they are brightest at redshift z~~2z \approx 2.
The redshift of signals emitted just before passage through the gravitational radius
Optical appearance of the collapsing star
(a) Schwarzschild
(b) Kruskal-Szekeres
Figure 32.1.
The free-fall collapse of a star of initial radius R_(i)=10MR_{i}=10 \mathrm{M} as depicted alternatively in (a) Schwarzschild coordinates, (b) Kruskal-Szekeres coordinates, and (c) ingoing Eddington-Finkelstein coordinates (see Box 31.2). The region of spacetime inside the collapsing star is grey, that outside it is white. Only the geometry of the exterior region is that of Schwarzschild. The curve separating the grey and white regions is the geodesic world line of the surface of the collapsing star (equations [31.10] or [32.10] with {:r_(max)=R_(i)=10M)\left.r_{\max }=R_{i}=10 \mathrm{M}\right). This world line is parameterized by proper time, tau\tau, as measured by an observer who sits on the surface of the star; the radial light cones, as calculated from ds^(2)=0d s^{2}=0, are attached to it.
Notice that, although the shapes of the light cones are not all the same relative to Schwarzschild coordinates or relative to Eddington-Finkelstein coordinates, they are all the same relative to KruskalSzekeres coordinates. This is because light rays travel along 45-degree lines in the u,vu, v-plane (dv=+-du)(d v= \pm d u), but they travel along curved paths in the r,tr, t-plane and r, widetilde(V)r, \widetilde{V}-plane.
The Kruskal-Szekeres spacetime diagram shown here is related to the Schwarzschild diagram by equations (31.13) plus a translation of Schwarzschild time: t longrightarrow t+42.8Mt \longrightarrow t+42.8 \mathrm{M}. The Eddington-Finklestein diagram is related to the Schwarzschild diagram by
widetilde(V)=t+r^(**)=t+r+2M ln |r//2M-1|\widetilde{V}=t+r^{*}=t+r+2 M \ln |r / 2 M-1|
(see Box 31.2).
It is evident from these diagrams that the free-fall collapse is characterized by a constantly diminishing radius, which drops from R=10 MR=10 M to R=0R=0 in a finite and short comoving proper time interval, Delta tau=35.1M\Delta \tau=35.1 \mathrm{M}. The point R=0R=0 and the entire region r=0r=0 outside the star make up a physical "singularity" at which infinite tidal gravitational forces-according to classical, unquantized general relativity-can and do crush matter to infinite density (see end of §31.2\S 31.2§; also §32.6\S 32.6§ ).
(c) Eddington-Finkelstein
The Eddington-Finkelstein diagram depicts a series of photons emitted radially from the surface of the collapsing star, and received by an observer at r=R_("initial ")=10Mr=R_{\text {initial }}=10 \mathrm{M}. The observer eventually receives all photons emitted radially from outside the gravitational radius; all photons emitted after the star passes through its gravitational radius eventually get pulled into the singularity at r=0r=0; and any photon emitted radially at the gravitational radius stays at the gravitational radius forever.
Non-free-fall collapse is similar to the collapse depicted here. When pressure gradients are present, only the detailed shape of the world line of the star's surface changes.
Notice how short is the characteristic ee-folding time for the decay of luminosity and for the radial redshift:
{:[tau_("char ")=2M~~1xx10^(-5)(M//M_(o.))" sec "],[(32.7)=((" light-travel time across a flat-space ")/(" distance equal to the gravitational radius ")).]:}\begin{align*}
\tau_{\text {char }} & =2 M \approx 1 \times 10^{-5}\left(M / M_{\odot}\right) \text { sec } \\
& =\binom{\text { light-travel time across a flat-space }}{\text { distance equal to the gravitational radius }} . \tag{32.7}
\end{align*}
Here M_(o.)M_{\odot} denotes one solar mass.
EXERCISE
The rest of this chapter is Track 2. No previous Track-2 material is needed as preparation for it, but it is needed as preparation for (1) the Track-2 part of Chapter 33 (black holes), and (2) Chapter 34 (singularities and global methods).
Exercise 32.2. REDSHIFTS DURING COLLAPSE
(a) Let a radio transmitter on the surface of a collapsing spherical star emit monochromatic waves of wavelength lambda_(e)\lambda_{e}; and let a distant observer, at the same theta,phi\theta, \phi, as the transmitter, receive the waves. Show that at late times the wavelength received varies as
[equation (32.6)], where tt is proper time as measured by the distant observer.
(b) [Track 2] Use kinetic theory for the outgoing photons (conservation of density in phase space: Liouville's theorem; $22.6\$ 22.6 ) to show that the energy flux of the radiation received (ergs //cm^(2)sec/ \mathrm{cm}^{2} \mathrm{sec} ) varies as
{:(32.8b)F prope^(-t//M):}\begin{equation*}
F \propto e^{-t / M} \tag{32.8b}
\end{equation*}
(c) Suppose that nuclear reactions at the center of the collapsing star generate neutrinos of energy E_(e)E_{e}, and that these neutrinos flow freely outward (negligible absorption in star). Show that the energy of the neutrinos received by a distant observer decreases at late times as
(d) Show that the flux of neutrino energy dies out at late times as
{:(32.9b)F prope^(-t//2M):}\begin{equation*}
F \propto e^{-t / 2 M} \tag{32.9b}
\end{equation*}
(e) Explain in elementary terms why the decay laws (32.8a) and (32.9a) for energy are the same, but the decay laws (32.8b) and (32.9b) for energy flux are different.
(f) Let a collapsing star emit photons from its surface at the black-body rate
(dN)/(d tau)=(1.5 xx10^(11)(" photons ")/(cm^(2)sec K^(3)))xx((" surface area ")/(" of star "))xx((" temperature ")/(" of surface "))^(3)\frac{d N}{d \tau}=\left(1.5 \times 10^{11} \frac{\text { photons }}{\mathrm{cm}^{2} \sec \mathrm{~K}^{3}}\right) \times\binom{\text { surface area }}{\text { of star }} \times\binom{\text { temperature }}{\text { of surface }}^{3}
Let a distant observer count the photons as they pass through his sphere of radius r≫Mr \gg M. Let him begin his count (time t=0t=0 ) when he sees (via photons traveling radially outward) the center of the star's surface pass through the radius r=3Mr=3 \mathrm{M}. Show that, in order of magnitude, the time he and his associates must wait, until the last photon that will ever get out has reached them, is
where T_(11)T_{11} is the star's surface temperature in units of 10^(11)K10^{11} \mathrm{~K}.
§32.4. COLLAPSE OF A STAR WITH UNIFORM DENSITY AND ZERO PRESSURE
When one turns attention to the interior of a collapsing star and to the precise world line that its surface follows in the Schwarzschild geometry, one encounters rather complicated mathematics. The simplest case to treat is that of a "star" with uniform density and zero pressure; and, indeed, until recently that was the only case which had been treated in detail. The original-and very complete-analysis of the collapse of such a uniform-density "ball of dust" was given in the classic paper of Oppenheimer and Snyder (1939). More recently, other workers have discussed it from slightly different points of view and using different coordinate systems. The approach taken here was devised by Beckedorff and Misner (1962).
Because no pressure gradients are present to deflect their motion, the particles on the surface of any ball of dust must move along radial geodesics in the exterior Schwarzschild geometry. For a ball that begins at rest with finite radius, R=R_(i)R=R_{i}, at time t=0t=0, the subsequent geodesic motion of its surface is given by equations (31.10):
Here RR is the Schwarzschild radial coordinate (i.e., 4piR^(2)4 \pi R^{2} is the star's surface area) at Schwarzschild time tt. This world line is plotted in Figure 32.1 for R_(i)=10MR_{i}=10 \mathrm{M}, in terms of Schwarzschild coordinates, Kruskal-Szekeres coordinates, and EddingtonFinkelstein coordinates. The proper time read by a clock on the surface of the collapsing star is given by equation (31.10b):
Note that the collapse begins when the parameter eta\eta is zero ( R=R_(i),t=tau=0R=R_{i}, t=\tau=0 ); and it terminates at the singularity ( R=0,eta=piR=0, \eta=\pi ) after a lapse of proper time, as measured on any test particle falling with the dust, equal to
It is interesting, though coincidental, that this is precisely the time-lapse required for free-fall collapse to infinite density in Newtonian theory [see equation (25.27'), Figure 25.3, and associated discussion].
What is the behavior of the interior of the ball of dust as it collapses? A variety of different interiors for pressureless dust can be conceived (exercise 32.8). But here attention focuses on the simplest of them: an interior that is homogeneous and isotropic everywhere, except at the surface-i.e., an interior locally identical to a dust-filled Friedmann cosmological model (Box 27.1). Is the Friedmann interior to be "open" (k=-1)(k=-1), "flat" (k=0)(k=0), or "closed" (k=+1)(k=+1) ? Only the closed case
The collapse, from rest, of a uniform-density ball of "dust":
(1) world line of ball's surface in exterior Schwarzschild coordinates
(2) interior of ball is identical to a portion of a closed Friedmann universe
is appropriate, since one has already demanded [equation (32.10)] that the star be at rest initially (initial rate of change of density equals zero; "moment of maximum expansion").
Using comoving hyperspherical coordinates, chi,theta,phi\chi, \theta, \phi, for the star's interior, and putting the origin of coordinates at the star's center, one can write the line element in the interior in the familiar Friedmann form
{:(32.11)ds^(2)=-dtau^(2)+a^(2)(tau)[dchi^(2)+sin^(2)chi(dtheta^(2)+sin^(2)theta dphi^(2))]:}\begin{equation*}
d s^{2}=-d \tau^{2}+a^{2}(\tau)\left[d \chi^{2}+\sin ^{2} \chi\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)\right] \tag{32.11}
\end{equation*}
Here a(tau)a(\tau) is given by the familiar cycloidal relation,
[see equations (1), (9), (4), and (5) of Box 27.1, with eta\eta replaced by eta+pi\eta+\pi ].
There is one possible difficulty with this interior solution. In the cosmological case, the solution was homogeneous and isotropic everywhere. Here homogeneity and isotropy are broken at the star's surface-which lies at some radius
for all tau\tau, as measured in terms of the hyperspherical polar angle chi\chi, a comoving coordinate (first picture in Box 27.2). At that surface (i.e., three-dimensional world tube enclosing the star's fluid) the interior Friedmann geometry must match smoothly onto the exterior Schwarzschild geometry. If the match cannot be achieved, then the Friedmann line element (32.11) cannot represent the interior of a collapsing star. An example of a case in which the matching could not be achieved is an interior of uniform and nonzero pressure, as well as uniform density. In that case there would be an infinite pressure gradient at the star's surface, which would blow off the outer layers of the star, and would send a rarefaction wave propagating inward toward its center. The uniform distribution of density and pressure would quickly be destroyed.
For the case of zero pressure, the match is possible. As a partial verification of the match, one can examine the separate and independent predictions made by the interior and exterior solutions for the star's circumference, C=2pi RC=2 \pi R, as a function of proper time tau\tau at the star's surface. The external Schwarzschild solution predicts the cycloidal relation,
A more complete verification of the match is given in exercise 32.4.
For further insight into this idealized model of gravitational collapse, see Box 32.1.
Exercise 32.3. EMBEDDING DIAGRAMS AND PHOTON PROPAGATION FOR COLLAPSING STAR
Verify in detail the features of homogeneous collapse described in Box 32.1.
Exercise 32.4. MATCH OF FRIEDMANN INTERIOR TO SCHWARZSCHILD EXTERIOR
The Einstein field equations are satisfied on a star's surface if and only if the intrinsic and extrinsic geometries of the surface's three-dimensional world tube are the same, whether measured on its interior or on its exterior (see §21.13\S 21.13§ for proof and discussion). Verify that for the collapsing star discussed above, the intrinsic and extrinsic geometries match at the join between the Friedmann interior and the Schwarzschild exterior. [Hints: (a) Use eta,theta,phi\eta, \theta, \phi, as coordinates on the world tube of the star's surface, and show that the intrinsic geometry has the same line element
{:(32.18a)ds^(2)=a^(2)(eta)[-deta^(2)+sin^(2)chi_(0)(dtheta^(2)+sin^(2)theta dphi^(2))]",":}\begin{equation*}
d s^{2}=a^{2}(\eta)\left[-d \eta^{2}+\sin ^{2} \chi_{0}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)\right], \tag{32.18a}
\end{equation*}
whether measured in the Schwarzschild exterior or in the Friedmann interior. (b) Show that the extrinsic curvature of the world tube has the same components
whether measured in the Schwarzschild exterior or in the Friedmann interior.]
Exercise 32.5. STARS THAT COLLAPSE FROM INFINITY
(a) Patch together a truncated Schwarzschild geometry and the geometry of a truncated "flat" ( k=0k=0 ), dust-filled Friedmann universe to obtain a model of a star that collapses from rest at an infinite initial radius. [Hint: The world line of the star's surface in the Schwarzschild geometry is given by equations (31.2).]
(b) Similarly patch together a truncated Schwarzschild geometry and the geometry of a truncated "open" (k=-1)(k=-1), dust-filled Friedmann universe to obtain a star which collapses from infinity with finite initial inward velocity.
Box 32.1 AN IDEALIZED COLLAPSING STAR WITH FRIEDMANN INTERIOR AND SCHWARZSCHILD EXTERIOR
(See §32.4 and exercises 32.3 and 32.4 for justification of the results described here.)
Initial State
(1) Take a Friedmann universe of radius a=a_(m)a=a_{m} at its moment of maximum expansion, eta=0\eta=0; and slice off and discard the region chi_(0) < chi <= pi\chi_{0}<\chi \leq \pi, where chi_(0)\chi_{0} is some angle less than pi//2\pi / 2. (2) Take a Schwarzschild geometry of mass M=(a_(m)//2)sin^(3)chi_(0)M=\left(a_{m} / 2\right) \sin ^{3} \chi_{0} at the moment t=0t=0; and slice off and discard the region r < R_(i)=a_(m)sin chi_(0)r<R_{i}=a_{m} \sin \chi_{0}. (3) Glue the remaining pieces of Friedmann and Schwarzschild geometry together smoothly along their cut surfaces. The resultant object will be a momentarily static star of uniform density rho_(i)=3//(8pia_(m)^(2))\rho_{i}=3 /\left(8 \pi a_{m}{ }^{2}\right), of mass M=(a_(m)//2)sin^(3)chi_(0)M=\left(a_{m} / 2\right) \sin ^{3} \chi_{0}, and of radius R_(i)=R_{i}=a_(m)sin chi_(0)a_{m} \sin \chi_{0}.
Subsequent Evolution
Release this star from its intial state, and let it collapse in accord with Einstein's field equations. The interior, truncated Friedmann universe and the exterior, truncated Schwarzschild geometry will evolve just as though they had never been cut up and patched together; and this evolution will preserve the smoothness of the match between interior and exterior!
Details of the Collapse
Probe the details of the collapse using sequences of embedding diagrams (histories ABCDA B C D and AWXYA W X Y ), and using photons that propagate radially outward (photons alpha,beta,gamma,delta,epsilon)\alpha, \beta, \gamma, \delta, \epsilon). The example shown here has chi_(0)=0.96\chi_{0}=0.96 and R_(i)//M=2//sin^(2)chi_(0)=3R_{i} / M=2 / \sin ^{2} \chi_{0}=3.
History of Collapse as Probed by Hypersurfaces ABCD:
(1) Initial configuration, A-A^(')-A^('')A-A^{\prime}-A^{\prime \prime}, is that constructed by cutting and sewing at times eta=t=0\eta=t=0.
(2) Each subsequent configuration has as its interior a slice of constant Friedmann time eta\eta.
(3) The interior remains always a spherical cup of half-angle chi_(0)\chi_{0}; but it contracts from radius R=R_(i)=a_(m)sin chi_(0)R=R_{i}=a_{m} \sin \chi_{0} to R=0R=0 as time increases.
(4) The matter in the star is all crushed simultaneously to infinite density when RR reaches zero, and the external Schwarzschild "funnel" develops a cusp-like singularity at that point.
(5) As time increases further, this cusp pulls the region r < 2Mr<2 M of the funnel into r=0r=0 so fast that the outward-traveling photon delta\delta is gobbled up and crushed.
These embedding diagrams must be rotated about the vertical axes in order to become 2-dimensional surfaces analogous to Figure 23.1.
Box 32.1 (continued)
History as Probed by Hypersurfaces AWXY
(1) Initial configuration, A-A^(')-A^('')A-A^{\prime}-A^{\prime \prime}, is again that constructed by cutting and sewing at eta=t=0\eta=t=0.
(2) Subsequent hypersurfaces are very different from eta=\eta= const.
(3) As time passes, a neck develops in the geometry just outside the surface of the star.
(4) This neck becomes tighter and tighter and then pinches off, leaving the star completely isolated from the rest of the universe, and leaving a deadly cusp-like singularity in the exterior geometry where the star used to be.
(5) The isolated star, in its own little closed universe, continues to contract until it is crushed to infinite density, while the external geometry begins to develop another neck and the cusp quickly gobbles up photon delta\delta.
X
The extreme difference between histories ABCDA B C D and AWXYA W X Y dramatizes the "many-fingered time" of general relativity. The hypersurface on which one explores the geometry can be pushed ahead faster in time in one region, at the option of the party of explorers. Thus whether one region of the star collapses first, or another, or the entire star collapses simultaneously, is a function both of the spacetime geometry and of the choice of slicing. The party of explorers has this choice of slicing at their own control, and thus they themselves to this extent govern what kind of spacelike slices they will see as their exploration moves forward in time. The spacetime geometry that they slice, however, is in no way theirs to control or to change. To the extent that their masses are negligible and they serve merely as test objects, they have no influence whatsoever on the spacetime. It was fixed completely by the specification of the initial conditions for the collapse. In brief, spacetime is four-dimensional and slices are only three-dimensional (and in the pictures here look only two-dimensional or one-dimensional). Any one set of slices captures only a one-sided view of the whole story. To see the entire picture one must either examine the dynamics of the geometry as it reveals itself in varied choices of the slicing or become accustomed to visualizing the spacetime geometry as a whole.
§32.5. SPHERICALLY SYMMETRIC COLLAPSE WITH INTERNAL PRESSURE FORCES
So far as the external gravitational field is concerned, the only difference between a freely collapsing star and a collapsing, spherically symmetric star with internal pressure is this: that the surfaces of the two stars move along different world lines in the exterior Schwarzschild geometry. Because the exterior geometry is the same in both cases, the qualitative aspects of free-fall collapse as described in the last section can be carried over directly to the case of nonnegligible internal pressure.
An important and fascinating question to ask is this: can large internal pressures in any way prevent a collapsing star from being crushed to infinite density by infinite tidal gravitational forces? From the Kruskal-Szekeres diagram of Figure 32.1,b, it is evident that, once a star has passed inside its gravitational radius ( R < 2MR<2 M ), no internal pressures, regardless of how large they may be, can prevent the star's surface from being crushed in a singularity. The surface must move along a time-like world line, and all such world lines inside r=2Mr=2 M hit r=0r=0. Although there is no such theorem now available, one can reasonably conjecture that, if the surface of a spherical configuration is crushed in the r=0r=0 singularity, the entire interior must also be crushed.
The details of the interior dynamics of a spherically symmetric collapsing star with pressure are not so well-understood as the exterior Schwarzschild dynamics. However, major advances in one's understanding of the interior dynamics are now being made by means of numerical computations and analytic analyses [see Misner (1969a) for a review]. In these computations and analyses, no new features (at least, no unexpected ones) have been encountered beyond those that occurred in the simple uniform-density, free-fall collapse of the last section.
Exercise 32.6. GENERAL SPHERICAL COLLAPSE: METRIC IN COMOVING COORDINATES
Consider an inhomogeneous star with pressure, undergoing spherical collapse. Spherical symmetry alone is enough to guarantee the existence of a Schwarzschild coordinate system (t,r,theta,phi)(t, r, \theta, \phi) throughout the interior and exterior of the star [see equation (32.2) and preceding discussion]. Label each spherical shell of the star by a parameter aa, which tells how many baryons are contained interior to that shell. Then r(a,t)r(a, t) is the world line of the shell with label aa. The expression for these world lines can be inverted to obtain a(t,r)a(t, r), the number of baryons interior to radius rr at time tt. Show that there exists a new time coordinate widetilde(t)(t,r)\widetilde{t}(t, r), such that the line element (32.2), rewritten in the coordinates ( widetilde(t),a,theta,phi)(\widetilde{t}, a, \theta, \phi), has the form
{:[(32.19a)ds^(2)=-e^(2 widetilde(Phi))d tilde(t)^(2)+[((del r//del a)_( tilde(t))da)/(Gamma)]^(2)+r^(2)dOmega^(2)","],[(32.19b) widetilde(Phi)= tilde(Phi)( tilde(t)","a)","quad r=r( tilde(t)","a)","quad Gamma=Gamma( tilde(t)","a).]:}\begin{align*}
& d s^{2}=-e^{2 \widetilde{\Phi}} d \tilde{t}^{2}+\left[\frac{(\partial r / \partial a)_{\tilde{t}} d a}{\Gamma}\right]^{2}+r^{2} d \Omega^{2}, \tag{32.19a}\\
& \widetilde{\Phi}=\tilde{\Phi}(\tilde{t}, a), \quad r=r(\tilde{t}, a), \quad \Gamma=\Gamma(\tilde{t}, a) . \tag{32.19b}
\end{align*}
Spherical collapse with pressure is qualitatively the same as without pressure
These are "comoving, synchronous coordinates" for the stellar interior.
Exercise 32.7. ADIABATIC SPHERICAL COLLAPSE: EQUATIONS OF EVOLUTION [Misner and Sharp (1964)]
Describe the interior of a collapsing star by the comoving, synchronous metric (32.19), by the number density of baryons nn, by the total density of mass-energy rho\rho, and by the pressure pp. The 4 -velocity of the star's fluid is
since the fluid is at rest in the coordinate system. Let a dot denote a proper time derivative as seen by the fluid-e.g.,
n^(˙)-=u[n]=e^(- widetilde(Phi))(del n//del widetilde(t))_(a)\dot{n} \equiv \boldsymbol{u}[n]=e^{-\widetilde{\Phi}}(\partial n / \partial \widetilde{t})_{a}
and let a prime denote a partial derivative with respect to baryon number,-e.g.
n^(')-=(del n//del a)_( tilde(t)).n^{\prime} \equiv(\partial n / \partial a)_{\tilde{t}} .
Denote by UU the rate of change of (1//2pi)xx(1 / 2 \pi) \times (circumference of shell), as measured by a man riding in a given shell:
{:(32.21a)U-=r^(˙):}\begin{equation*}
U \equiv \dot{r} \tag{32.21a}
\end{equation*}
and denote by m( tilde(t),a)m(\tilde{t}, a) the "total mass-energy interior to shell aa at time tilde(t)\tilde{t} :
{:(32.2lb)m( widetilde(t)","a)-=int_(0)^(a)4pir^(2)rho( widetilde(t)","a)r^(')da:}\begin{equation*}
m(\widetilde{t}, a) \equiv \int_{0}^{a} 4 \pi r^{2} \rho(\widetilde{t}, a) r^{\prime} d a \tag{32.2lb}
\end{equation*}
(See Box 23.1 for discussion of this method of localizing mass-energy.) Assume that the collapse is adiabatic (no energy flow between adjacent shells; stress-energy tensor entirely that of a perfect fluid).
(a) Show that the equations of collapse [baryon conservation, (22.3); local energy conservation, (22.11a); Euler equation, (22.13); and Einstein field equations (ex. 14.16)] can be reduced to the following eight equations for the eight functions widetilde(Phi),Gamma,r,n,rho,p,U,m\widetilde{\Phi}, \Gamma, r, n, \rho, p, U, m :
{:(32.22h){:[r^(˙),=U,,(" dynamic equation for "r);],[((nr^(2))^(*))/(nr^(2)),=-(U^('))/(r^(')),," (dynamic equation for "n);],[((rho^(˙)))/(rho+p),=((n^(˙)))/(n)","quad,,{:[" except at a shock front, where adiabaticity "],[" breaks down "]:}],[U^(˙),=-(Gamma^(2))/(rho+p)(p^('))/(r^('))-(m+4pir^(3)p)/(r^(2)),," (dynamic equation for "rho" ); "],[p,=p(n","rho),," (equation of state); "],[ widetilde(Phi)^('),=-p^(')//(rho+p)"," widetilde(Phi)=0" at star's surface ",," (source equation for " widetilde(Phi)" ); "],[m^('),=4pir^(2)rhor^(')",",m=0" at "a=0",",],[Gamma," (source equation for "m" ); "],[Gamma,sign(r^('))(1+U^(2)-2m//r)^(1//2),," (algebraic equation for "Gamma" ). "]:}:}\begin{array}{rlrl}
\dot{r} & =U & & (\text { dynamic equation for } r) ; \\
\frac{\left(n r^{2}\right)^{\cdot}}{n r^{2}} & =-\frac{U^{\prime}}{r^{\prime}} & & \text { (dynamic equation for } n) ; \\
\frac{\dot{\rho}}{\rho+p} & =\frac{\dot{n}}{n}, \quad & & \begin{array}{l}
\text { except at a shock front, where adiabaticity } \\
\text { breaks down }
\end{array} \\
\dot{U} & =-\frac{\Gamma^{2}}{\rho+p} \frac{p^{\prime}}{r^{\prime}}-\frac{m+4 \pi r^{3} p}{r^{2}} & & \text { (dynamic equation for } \rho \text { ); } \\
p & =p(n, \rho) & & \text { (equation of state); } \\
\widetilde{\Phi}^{\prime} & =-p^{\prime} /(\rho+p), \widetilde{\Phi}=0 \text { at star's surface } & & \text { (source equation for } \widetilde{\Phi} \text { ); } \\
m^{\prime} & =4 \pi r^{2} \rho r^{\prime}, & m=0 \text { at } a=0, & \\
\Gamma & \text { (source equation for } m \text { ); } \tag{32.22h}\\
\Gamma & \operatorname{sign}\left(r^{\prime}\right)\left(1+U^{2}-2 m / r\right)^{1 / 2} & & \text { (algebraic equation for } \Gamma \text { ). }
\end{array}
(b) The preceding equations are in a form useful for numerical calculations. [For particular numerical solutions and for the handling of shocks, see May and White (1966).] For analytic work it is often useful to replace (32.22b) by
{:('")"m^(˙)=-4pir^(2)pU:}\begin{equation*}
\dot{m}=-4 \pi r^{2} p U \tag{$\prime$}
\end{equation*}
Derive these equations.
(c) Explain why equations (32.22g)(32.22 \mathrm{~g}) and (32.22d^('))\left(32.22 \mathrm{~d}^{\prime}\right) justify the remarks made in Box 23.1 about localizability of energy.
Exercise 32.8. ANALYTIC SOLUTIONS FOR PRESSURE-FREE COLLAPSE [Tolman (1934b); Datt (1938)]
Show that the general solution to equations (32.22) in the case of zero pressure can be generated as follows.
(a) Specify the mass inside shell a,m(a)a, m(a); by equation ( 32.22d^(')32.22 \mathrm{~d}^{\prime} ), with p=0p=0, it will not change with time tilde(t)\tilde{t}.
(b) Assume that all the dust particles have rest masses mu\mu that depend upon radius, mu(a)\mu(a); so
{:(32.23a)rho=mu n:}\begin{equation*}
\rho=\mu n \tag{32.23a}
\end{equation*}
it will be independent of tilde(t)\tilde{t}.
(d) Specify an initial distribution of circumference 2pi r2 \pi r as function of aa, and solve the dynamic equation
to obtain the subsequent evolution of r( tilde(t),a)r(\tilde{t}, a). Notice that this equation has identically the same form as in Newtonian theory!
(e) Calculate the remaining quantities of interest from the algebraic equations
{:[(32.23~d)ds^(2)=- tilde(dt^(2))+(r^(')da//Gamma)^(2)+r^(2)dOmega^(2)],[(32.23e)rho=mu n=m^(')//(4pir^(2)r^('))],[(32.23f) widetilde(Phi)=0","quad U=del r//del widetilde(t)]:}\begin{gather*}
d s^{2}=-\tilde{d t^{2}}+\left(r^{\prime} d a / \Gamma\right)^{2}+r^{2} d \Omega^{2} \tag{32.23~d}\\
\rho=\mu n=m^{\prime} /\left(4 \pi r^{2} r^{\prime}\right) \tag{32.23e}\\
\widetilde{\Phi}=0, \quad U=\partial r / \partial \widetilde{t} \tag{32.23f}
\end{gather*}
[Note: In this solution, successive "shells" may pass through each other, producing a surface of infinite density as they do (r^(')longrightarrow0:}\left(r^{\prime} \longrightarrow 0\right. where {:m^(')!=0)\left.m^{\prime} \neq 0\right), since there is no pressure built up to stop shell crossing. When this happens, the coordinate system becomes pathological ( aa no longer increases monotonically outward), but spacetime remains well-behaved. The surface of infinite density (1) produces negligible tidal forces on neighboring dust particles; and (2) like the surface layers of $21.13\$ 21.13, it is an idealization that gets smeared down to finite density by finite pressure.]
Exercise 32.9. COLLAPSE WITH UNIFORM DENSITY
Recover the Friedmann-Schwarzschild solution for collapse with uniform density and zero pressure by specifying appropriate forms for m(a)m(a) and r(a)r(a) in the prescription of exercise 32.8. In the interior of the star, give the dust particles nonzero rest masses, mu=\mu= constant !=0\neq 0; in the exterior give them zero rest masses, mu=0\mu=0 ("imaginary dust particles" in vacuum). Reduce the resulting metric ( 32.23 d ) to that of Friedmann inside the star, and to that of Novikov for the Schwarzschild geometry outside the star [equations (31.12)].
§32.6. THE FATE OF A MAN WHO FALLS INTO THE SINGULARITY AT r=0r=0
The effect of tidal forces on the body of a man falling into the r=0r=0 singularity:
Stage 1: body resists deformation; stresses build up
Consider the plight of an experimental astrophysicist who stands on the surface of a freely falling star as it collapses to R=0R=0.
As the collapse proceeds toward R=0R=0, the various parts of the astrophysicist's body experience different gravitational forces. His feet, which are on the surface of the star, are attracted toward the star's center by an infinitely mounting gravitational force; while his head, which is farther away, is accelerated downward by a somewhat smaller, though ever rising force. The difference between the two accelerations (tidal force) mounts higher and higher as the collapse proceeds, finally becoming infinite as RR reaches zero. The astrophysicist's body, which cannot withstand such extreme forces, suffers unlimited stretching between head and foot as RR drops to zero.
But this is not all. Simultaneous with this head-to-foot stretching, the astrophysicist is pulled by the gravitational field into regions of spacetime with ever-decreasing circumferential area, 4pir^(2)4 \pi r^{2}. In order to accomplish this, tidal gravitational forces must compress the astrophysicist on all sides as they stretch him from head to foot. The circumferential compression is actually more extreme than the longitudinal stretching; so the astrophysicist, in the limit R longrightarrow0R \longrightarrow 0, is crushed to zero volume and indefinitely extended length.
The above discussion can be put on a mathematical footing as follows.
There are three stages in the killing of the astrophysicist: (1) the early stage, when his body successfully resists the tidal forces; (2) the intermediate stage, when it is gradually succumbing; and (3) the final stage, when it has been completely overwhelmed.
During the early stage, one can analyze the tidal forces by means of the equation of geodesic deviation, evaluated in the astrophysicist's orthonormal frame omega^( hat(gamma)),omega^( hat(rho))\boldsymbol{\omega}^{\hat{\gamma}}, \boldsymbol{\omega}^{\hat{\rho}}, omega^( hat(theta)),omega^( hat(phi))\boldsymbol{\omega}^{\hat{\theta}}, \boldsymbol{\omega}^{\hat{\phi}} (see §31.2\S 31.2§ ). In this frame, the nonvanishing components of the Riemann curvature tensor are given by equations (31.6):
The equation of geodesic deviation says that two freely moving particles, momentarily at rest in the astrophysicist's local inertial frame, and separated by the 3-vector
Using the components (32.24a) of the curvature tensor, one sees that
{:[D^(2)xi^( hat(rho))//dtau^(2)=+(2M//r^(3))xi^( hat(rho))],[(32.24b)D^(2)xi^( hat(theta))//dtau^(2)=-(M//r^(3))xi^( hat(theta))],[D^(2)xi^( hat(rho))//dtau^(2)=-(M//r^(3))xi^( hat(phi))]:}\begin{align*}
& D^{2} \xi^{\hat{\rho}} / d \tau^{2}=+\left(2 M / r^{3}\right) \xi^{\hat{\rho}} \\
& D^{2} \xi^{\hat{\theta}} / d \tau^{2}=-\left(M / r^{3}\right) \xi^{\hat{\theta}} \tag{32.24b}\\
& D^{2} \xi^{\hat{\rho}} / d \tau^{2}=-\left(M / r^{3}\right) \xi^{\hat{\phi}}
\end{align*}
To apply these equations to the astrophysicist's body, idealize it (for simplicity) as a homogeneous rectangular box of mass mu~~165\mu \approx 165 pounds ~~75kg\approx 75 \mathrm{~kg}, of length ℓ~~70\ell \approx 70 inches ~~1.8m\approx 1.8 \mathrm{~m} in the e_( hat(rho))\boldsymbol{e}_{\hat{\rho}} direction, and of width and depth w~~10w \approx 10 inches ~~0.2m\approx 0.2 \mathrm{~m} in the e_( hat(theta))\boldsymbol{e}_{\hat{\theta}} and e_( hat(phi))\boldsymbol{e}_{\hat{\phi}} directions. Then calculate the stresses that must be set up in this idealized body to prevent its particles from moving along diverging (and converging) geodesics.
From the form of equations (32.24), it is evident that the principal directions of the stress will be e_( hat(rho)),e_( hat(theta))\boldsymbol{e}_{\hat{\rho}}, \boldsymbol{e}_{\hat{\theta}}, and e_( hat(phi))\boldsymbol{e}_{\hat{\phi}} (i.e., in the e_( hat(rho)),e_( hat(theta)),e_( hat(phi))\boldsymbol{e}_{\hat{\rho}}, \boldsymbol{e}_{\hat{\theta}}, \boldsymbol{e}_{\hat{\phi}} basis, the stress tensor will be diagonal). The longitudinal component of the stress, at the astrophysicist's center of mass, can be evaluated as follows. A volume element of his body with mass d mud \mu, located at a height hh above the center of mass (distance hh measured along e_( hat(rho))\boldsymbol{e}_{\hat{\rho}} direction) would accelerate with a=(2M//r^(3))ha=\left(2 M / r^{3}\right) h away from the center of mass, if it were allowed to move freely. To prevent this acceleration, the astrophysicist's muscles must exert a force
dF=ad mu=(2M//r^(3))hd mud F=a d \mu=\left(2 M / r^{3}\right) h d \mu
This force contributes to the stress across the horizontal plane ( e_( hat(theta))^^e_( hat(phi))\boldsymbol{e}_{\hat{\theta}} \wedge \boldsymbol{e}_{\hat{\phi}} plane) through the center of mass. The total force across that plane is the sum of the forces on all mass elements above it (which is also equal to the sum of the forces on the mass elements below it):
{:[F=int_((region above plane) )ad mu=int_(0)^(1//2)((2Mh)/(r^(3)))((mu)/(ℓw^(2)))(w^(2)dh)],[=(1)/(4)(mu Mℓ)/(r^(3)).]:}\begin{aligned}
F & =\int_{\text {(region above plane) }} a d \mu=\int_{0}^{1 / 2}\left(\frac{2 M h}{r^{3}}\right)\left(\frac{\mu}{\ell w^{2}}\right)\left(w^{2} d h\right) \\
& =\frac{1}{4} \frac{\mu M \ell}{r^{3}} .
\end{aligned}
The stress is this force divided by the cross-sectional area w^(2)w^{2}, with a minus sign because it is a tension rather than a pressure:
The components of the stress in the e_( hat(theta))\boldsymbol{e}_{\hat{\theta}} and e_( hat(phi))\boldsymbol{e}_{\hat{\phi}} directions at the center of mass are, similarly,
(Recall that one atmosphere of pressure is 1.01 xx10^(6)1.01 \times 10^{6} dynes //cm^(2)/ \mathrm{cm}^{2}.)
Stage 2: body gives way; man dies
Stage 3: body gets crushed and distended
The human body cannot withstand a tension or pressure of >= 100\geq 100 atmospheres ~~10^(8)\approx 10^{8} dynes //cm^(2)/ \mathrm{cm}^{2} without breaking. Consequently, an astrophysicist on a freely collapsing star of one solar mass will be killed by tidal forces when the star's radius is R∼200km≫2M~~3kmR \sim 200 \mathrm{~km} \gg 2 M \approx 3 \mathrm{~km}.
By the time the star is much smaller than its gravitational radius, the baryons of the astrophysicist's body are moving along geodesics; his muscles and bones have completely given way. In this final stage of collapse, the timelike geodesics are curves along which the Schwarzschild "time"-coordinate, tt, is almost constant [ cfc f. the narrowing down of the light cones near r=0r=0 in Figure 32.1, a; also equation (31.2) in the limit r≪2M]r \ll 2 M]. The astrophysicist's feet touch the star's surface at one value of tt-say t=t_(f)t=t_{f}-while his head moves along the curve t=t_(h) > t_(f)t=t_{h}>t_{f}. Consequently, the length of the astrophysicist's body increases according to the formula
Here tau=[-int^(R)|g_(rr)|^(1//2)dr+:}\tau=\left[-\int^{R}\left|g_{r r}\right|^{1 / 2} d r+\right. constant ]] is proper time as it would be measured by the astrophysicist if he were still alive, and tau_("collapse ")\tau_{\text {collapse }} is the time at which he hits r=0r=0. The gravitational field also constrains the baryons of the astrophysicist's body to fall along world lines of constant theta\theta and phi\phi during the final stages of collapse. Consequently, his cross-sectional area decreases according to the law
By combining equations ( 32.26a,b32.26 \mathrm{a}, \mathrm{b} ), one sees that the volume of the astrophysicist's body decreases, during the last few moments of collapse, according to the law
This crushing of matter to infinite density by infinitely large tidal gravitational forces can occur not only on the surface of the collapsing star, but also at any other point along the r=0r=0 singularity outside the surface of the star. Hence, any foolish rocketeer who ventures below the radius r=2Mr=2 M of the external gravitational field is doomed to destruction.
For further discussion of spacetime singularities, and of the possibility that quantum gravitational effects might force a reconsideration of the singularities predicted by classical gravitation theory, see Chapter 30, §34.6\S 34.6§, and Chapter 44.
Instability, implosion, horizon, and singularity; these are the key stages in the spherical collapse of any star. Instability: The star, having exhausted its nuclear fuel, and having contracted slowly inward, begins to squeeze its pressure-sustaining electrons or photons onto its atomic nuclei; this softens the equation of state, which induces an instability [see, e.g., §§10.15\S \S 10.15§§ and 11.4 of Zel'dovich and Novikov (1971)
for details]. Implosion: Within a fraction of a second the instability develops into a full-scale implosion; for realistic density distributions, the stellar core falls rapidly inward on itself, and the outer envelopes trail along behind [see, e.g., the numerical calculations of Colgate and White (1966), Arnett (1966, 1967), May and White (1966), and Ivanova, Imshennik, and Nadezhin (1969)]. Horizon: In the idealized spherical case, the star's surface falls through its gravitational radius ("horizon"; end of communication with the exterior; point of no return). From the star's vantage point this happens after a finite, short lapse of proper time. But from an external vantage point the star requires infinite time to reach the horizon, though it becomes black exponentially rapidly in the process [e-folding time ∼M∼10^(-5)(M//M_(o.))\sim M \sim 10^{-5}\left(M / M_{\odot}\right) sec]. The result is a "black hole", whose boundary is the horizon (gravitational radius), and whose interior can never communicate with the exterior. Singularity: From the star's interior vantage point, within a short proper time interval Delta tau∼M∼10^(-5)(M//M_(o.))\Delta \tau \sim M \sim 10^{-5}\left(M / M_{\odot}\right) sec after passing through the horizon, a singularity is reached (zero radius, infinite density, infinite tidal gravitational forces).
Does this basic picture-instability, implosion, horizon, singularity-have any relevance for real stars? Might complications such as rotation, nonsphericity, magnetic fields, and neutrino fluxes alter the qualitative picture? No, not for small initial perturbations from sphericity. Perturbation theory analyses described in Box 32.2 and exercise 32.10 show that realistic, almost-spherically symmetric collapse, like idealized collapse, is characterized by instability, implosion, horizon; and Penrose (1965b; see $34.6\$ 34.6 ) proves that some type of singularity then follows.
Highly nonspherical collapse is more poorly understood, of course. Nevertheless, a number of detailed calculations and precise theorems point with some confidence to two conclusions: (1) horizons (probably) form when and only when a mass MM gets compacted into a region whose circumference in EVERY direction is C <= 4pi M\mathcal{C} \leq 4 \pi M (Box 32.3)32.3); (2) the external gravitational field of a horizon (black hole), after all the "dust" and gravitational waves have cleared away, is almost certainly the Kerr-Newman generalization of the Schwarzschild geometry (Chapter 33). If so, then the external field is determined uniquely by the mass, charge, and angular momentum that went "down the hole." (This nearly proved theorem carries the colloquial title "A black hole has no hair.")
The interior of the horizon, and the endpoint (if any) of the collapse are very poorly understood today. The various possibilities will be reviewed in Chapter 34. That a singularity occurs one can state with much certainty, thanks to theorems of Penrose, Hawking, and Geroch. But whether all, only some, or none of the collapsing matter and fields ultimately encounter the singularity one does not know.
EXERCISES
Summary of 1972
knowledge about realistic, nonspherical collapse:
(1) horizon
(2) black hole
(3) singularity
Exercise 32.10. PRICE'S THEOREM FOR A SCALAR FIELD
Exercise 32.10. PRICE'S THEOREM FOR A SCALAR FIELD [See Price (1971, 1972a), also Thorne (1972), for more details than are presented here.]
A collapsing spherical star, with an arbitrary nonspherical "scalar charge distribution," generates an external scalar field Phi\Phi. The vacuum field equation for Phi\Phi is ◻Phi=Phi_(;alpha)^(alpha)=0\square \Phi=\Phi_{; \alpha}{ }^{\alpha}=0. Ignore the back-reaction of the field's stress-energy on the geometry of spacetime.
Box 32.2 COLLAPSE WITH SMALL NON-SPHERICAL PERTURBATIONS [based on detailed calculations by Richard H. Price (1971, 1972a,b)].
A. Density Perturbations
When star begins to collapse, it possesses a small nonspherical "lump" in its density distribution.
As collapse proceeds, lump grows larger and larger [instability of collapse against small per-turbations-a phenomenon well known in Newtonian theory; see, e.g., Hunter (1967); Lin, Mestel, and Shu (1965)].
The growing lump radiates gravitational waves.
Waves of short wavelength (lambda≪M)(\lambda \ll M), emitted from near horizon (r-2M <= M)(r-2 M \leqq M), partly propagate to infinity and partly get backscattered by the "background" Schwarzschild curvature of spacetime. Backscattered waves propagate into horizon (surface of black hole; gravitational radius) formed by collapsing star.
Waves of long wavelength (lambda≫M)(\lambda \gg M), emitted from near horizon (r-2M <= M)(r-2 M \leqq M), get fully backscattered by spacetime curvature; they never reach out beyond r∼3Mr \sim 3 M; they end up propagating "down the hole."
Is lump on star still there as star plunges through horizon, and does star thereby create a deformed (lumpy) horizon? Yes, according to calculations.
But external observers can only learn about existence of "final lump" by examining deformation (quadrupole moment) in final gravitational field. That final deformation in field does not and cannot propagate outward with infinite speed (no instantaneous "action at a distance"). It propagates with speed of light, in form of gravitational waves with near-infinite wavelength (infinite redshift from edge of horizon to any external radius). Deformation in final field, like any other wave of long wavelength, gets fully backscattered by curvature of spacetime at r <= 3Mr \leqq 3 M; it cannot reach external observers. External observers can never learn of existence
of final lump. Final external field is perfectly spherical, lump-free, Schwarzschild geometry!
Even in region of backscatter (2M < r <= 3M)(2 M<r \leqq 3 M), final external field is lump-free. Backscattered waves, carrying information about existence of final lump, interfere destructively with outgoing waves carrying same information. Result is destruction of all deformation in external field and in horizon!
Final black hole is a Schwarzschild black hole!
B. Perturbations in Angular Momentum
When star begins to collapse, it possesses a small, nonzero intrinsic angular momentum ("spin") SS.
As collapse proceeds, S\boldsymbol{S} is conserved (except for a tiny, negligible change due to angular momentum carried off by waves; that change is proportional to square of amplitude of waves, i.e., to square of amplitude of perturbations in star, i.e., to S^(2)\boldsymbol{S}^{2} ).
Consequently, external field always and everywhere carries imprint of angular momentum SS (on imprints, see Chapter 19). There is no need for that imprint to propagate outward from near horizon. Moreover, it could not so propagate even if it tried, because of the conservation law for S\boldsymbol{S} (absence of dipole gravitational waves; see $§36.1\$ \S 36.1§ and 36.10 ).
Hence, the final external field is that of an undeformed, slowly rotating black hole:
{:[ds^(2)=-ubrace((1-(2M)/(r))dt^(2)+(dr^(2))/(1-2M//r)+r^(2)dOmega^(2)ubrace)_("Schwarzschild geometry ")],[-ubrace((4S sin theta)/(r^(2)))(r sin theta d phi)dt.ubrace)]:}\begin{aligned}
d s^{2}=- & \underbrace{\left(1-\frac{2 M}{r}\right) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2} d \Omega^{2}}_{\text {Schwarzschild geometry }} \\
& -\underbrace{\left.\frac{4 S \sin \theta}{r^{2}}\right)(r \sin \theta d \phi) d t .}
\end{aligned}
rotational imprint, see exercise 26.1; also Chapter 19.
Here the polar axis has been oriented along S\boldsymbol{S}.
C. Perturbations in Electromagnetic Field
Star possesses a magnetic field generated by currents in its interior, and an electric field due to an arbitrary internal charge distribution; and electromagnetic radiation is emitted by its hot matter. For simplicity, S\boldsymbol{S} is assumed zero.
Evolution of external electromagnetic field is similar to evolution of perturbations in external gravitational field. Distant observer can never learn "final" values of changeable quantities (magnetic dipole moment, electric dipole moment, quadrupole moments, . . .). Final values try to propagate out from horizon, carried by electromagnetic waves of near-infinite wavelength. But they cannot get out: spacetime curvature reflects them back down the hole; and they superpose destructively with their outgoing counterparts, to produce zero net field.
By contrast with all other quantities, which are changeable, the electric monopole moment (total flux of electric field; equal to 4pi4 \pi times total electric charge) is conserved. It never changes from before star collapses, through the collapse stage, into the quiescent black-hole stage.
Hence, the final external electromagnetic field is a spherically symmetric coulomb field {:E=(Q//r^(2))e_(r))\left.\boldsymbol{E}=\left(Q / r^{2}\right) \boldsymbol{e}_{r}\right) as measured by static B=0quad\boldsymbol{B}=0 \quad observer (r,theta,phi(r, \theta, \phi, constant ));
and the final spacetime geometry is that of Reissner and Nordstrom (charged black hole; see exercises 31.8 and 32.1 ):
{:[ds^(2)=-(1{:-(2M)/(r)+(Q^(2))/(r^(2)))dt^(2)],[+(dr^(2))/((1-2M//r+Q^(2)//r^(2)))+r^(2)dOmega^(2).]:}\begin{aligned}
d s^{2}=-(1 & \left.-\frac{2 M}{r}+\frac{Q^{2}}{r^{2}}\right) d t^{2} \\
& +\frac{d r^{2}}{\left(1-2 M / r+Q^{2} / r^{2}\right)}+r^{2} d \Omega^{2} .
\end{aligned}
D. Generalization; Price's Theorem
Let the star generate a "zero-rest-mass, integerspin field." ["Zero rest mass" refers to the quantized particles associated with the classical field. Classically it means the field has a Cou-lomb-law ( 1//r1 / r ) fall off at large distances. The spin also is a property of the quantized particles; classically it is most easily visualized as describing the symmetries of a monochromatic plane wave under rotations about the direction of propagation; see §35.6\S 35.6§. A scalar field has spin zero; an electromagnetic field has spin one; Einstein's gravitational field has spin two; . . . . Of such fields, only gravitational (s=2)(s=2) and electromagnetic ( s=1s=1 ) are known to exist in the real universe. See, e.g., Dirac (1936), Gårding (1945), Bargmann and Wigner (1948), Penrose (1965a), for further discussion.]
Let the spin-s field be sufficiently weak that its stress-energy perturbs the star's external, Schwarzschild geometry only very slightly.
Resolve the external field into spherical harmonics (scalar spherical harmonics for s=0s=0; vector spherical harmonics for s=1s=1; tensor spherical harmonics for s >= 2s \geq 2 ); and label the spherical harmonics by the usual integer ℓ(ℓ=\ell(\ell= 0 for monopole; ℓ=1\ell=1 for dipole; ℓ=2\ell=2 for quadrupole; etc.).
All multipole fields with ℓ < s\ell<s are conserved during the collapse (theorem from classical radiation theory). A scalar field (s=0)(s=0) conserves
nothing. The electromagnetic field (s=1)(s=1) conserves only its monopole parts (electric Coulomb field, and vanishing magnetic Coulomb field). The gravitational field (s=2)(s=2) conserves its monopole part (with imprint equal to mass), and its dipole parts (with imprints measuring the angular momentum, and the standard gravitational dipole moment-which vanishes if coordinate system is centered on star).
For ℓ >= s\ell \geq s, and only for ℓ >= s\ell \geq s, radiation is possible (scalar waves can have any multipolarity; electromagnetic waves must be dipole and higher; gravitational waves must be quadrupole and higher; see §36.1).
Price's theorem states that, as the nearly spherical star collapses to form a black hole, all things that can be radiated (all multipoles ℓ >= s\ell \geq s ) get radiated completely away-in part "off to infinity"; in part "down the hole" ("what is permitted is compulsory"). The final field is characterized completely by its conserved quantities (multipole moments with ℓ < s\ell<s ).
For proof of Price's theorem in the case of a scalar field, see exercise 32.10.
E. Generalization to Nonclassical Fields
See Hartle (1971,1972)(1971,1972) and Teitelboim (1972b,c) for neutrino fields; Bekenstein (1972a,b) and Teitelboim (1972a) for pion fields.
Box 32.3 COLLAPSE IN ONE AND TWO DIMENSIONS
A. The Question
To produce a black hole (horizon from which nothing can emerge), must one compact matter strongly in all three spatial dimensions, to circumferences C <= 4pi M\mathcal{C} \leqq 4 \pi M (quasispherical compaction); or is it sufficient to compact only in one or two dimensions?
B. The Answer for One Dimension
Consider, as an example readily generalized, the gravitational collapse of a spheroid of dust (zero pressure). Let the spheroid be highly Newtonian ( r⋙2Mr \ggg 2 M ) in its initial, momentary state of rest; and let it be slightly flattened (oblate). In Newtonian theory, any homogeneous, nonrotating spheroid of dust remains homogeneous as it collapses; but its deformations grow [see, e.g., Lin, Mestel, and Shu (1965) for details]. Hence, the spheroid of interest implodes to form a pancake of infinite density but finite mass per unit surface area. The final kinetic energy of the dust particles is roughly equal to their final potential energy:
{:[(1)/(2)v^(2)∼(M)/((C//2pi))],[M=" mass of spheroid, "],[C=" circumference of final pancake. "]:}\begin{aligned}
& \frac{1}{2} v^{2} \sim \frac{M}{(\mathcal{C} / 2 \pi)} \\
& M=\text { mass of spheroid, } \\
& \mathcal{C}=\text { circumference of final pancake. }
\end{aligned}
Consequently, so long as C//2pi≫2M\mathcal{C} / 2 \pi \gg 2 M, the collapse velocities remain much smaller than light, and the gravitational energy remains much smaller than the rest mass-energy. This means that for C//2pi≫\mathcal{C} / 2 \pi \gg2M2 M, the Newtonian analysis is an excellent approximation to general relativity all the way down to the pancake endpoint. Hence, no horizon can form, hardly any gravitational waves are emitted, and the whole story is exceedingly simple and fully Newtonian. However, since the pancake endpoint is not a singularity of spacetime (see the remarks at end of exercise 32.8), the evolution can proceed beyond it; and as E\mathcal{E} contracts to <= 4pi M\leqq 4 \pi M, the evolu-
tion will become very complicated and highly relativistic (see the "collapse, pursuit, and plunge scenario" of Figure 24.3).
C. The Answer for Two Dimensions
Consider, as an example not so readily generalized, the gravitational collapse of a homogeneous prolate spheroid of dust, initially highly Newtonian. Such a spheroid collapses to form a thin "thread" or "spindle" [see Lin, Mestel, and Shu (1965)]. Assume that the spheroid is still Newtonian when its threadlike state is reached. It then has a length ℓ\ell, a mass per unit length lambda=M//ℓ≪1\lambda=M / \ell \ll 1, and a rapidly contracting equatorial radius R≪ℓR \ll \ell. Subsequently, each segment of the thread collapses radially as though it were part of an infinite cylinder. [Ignore the instability of breakup into "beads"; see, e.g., Hunter (1967), Chandrasekhar (1968).] The radial collapse velocity approaches the speed of light and the gravitational energy approaches the rest mass-energy only when the thread has become exceedingly thin, R≲R_("crit ")∼ℓexp(-1//4lambda)R \lesssim R_{\text {crit }} \sim \ell \exp (-1 / 4 \lambda). At this stage, relativistic deviations from Newtonian collapse come into play. Thorne (1972) and Morgan and Thorne (1973) have analyzed the relativistic effects using an idealized infinite-cylinder model. The results are very different from either the spherical case or the pancake case. The collapsing cylinder emits a large flux of gravitational waves; but they are powerless to halt the collapse. The collapse proceeds inward to a thread-like singularity, without the creation of any horizon (no black hole!).
D. Objection to the Answer, a Reply, and a Conjecture
One can object that the collapses of both pancake and cylinder can be halted short of their endpoints, especially that of the pancake. As the thickness of
Box 32.3 (continued)
the pancake approaches zero, the vertical pull of gravity remains finite, but the pressure gradient caused by any finite pressure goes to infinity. Hence, pressure halts the collapse. Subsequently the rim of the pancake contracts toward the relativistic regime C//2pi <= 2M\mathcal{C} / 2 \pi \leq 2 M. In the collapse of a cylinder according to Newtonian theory, with a pressure-density relation p proprho^(gamma)p \propto \rho^{\gamma}, the gravitational acceleration a_(g)a_{g} and pressure-buoyancy acceleration a_(p)a_{p} vary as
Hence, for gamma > 1\gamma>1 (the most common realistic case) pressure halts the collapse, but for gamma < 1\gamma<1 it does
not. Whether this is true also after the relativistic domain is reached, one does not yet know.
Actually, the ability of pressure to halt the collapse is of no importance to the issue of black holes and horizons. The important thing is that in oblate collapse with final circumference C > 4pi M\mathcal{C}>4 \pi M, and also in prolate collapse with final thread length ℓ≫2M\ell \gg 2 M, no horizons are created. This, coupled with the omnipresent horizons in nearly spherical collapse (Box 32.2) suggests the following conjecture [Thorne (1972)]: Black holes with horizons form when and only when a mass MM gets compacted into a region whose circumference in EVER Y direction is C <= 4pi M\mathcal{C} \leq 4 \pi M. (Like most conjectures, this one is sufficiently vague to leave room for many different mathematical formulations!)
(a) Resolve the external field into scalar spherical harmonics, using Schwarzschild coordinates for the external Schwarzschild geometry:
Part of this effective potential [ℓ(ℓ+1)//r^(2)]\left[\ell(\ell+1) / r^{2}\right] is the "centrifugal barrier," and part [2M//r][2 M / r] is due to the curvature of spacetime. Notice the similarity of this effective potential for scalar waves, to the effective potentials for particles and photons moving in the Schwarzschild geometry,
(Boxes 25.6 and 25.7). The scalar-wave potential, like the photon potential, is positive for all r > 2Mr>2 M. It rises, from 0 at r=2Mr=2 M, to a barrier summit; then falls back to 0 at r=oor=\infty.
(b) Show that there exist no physically acceptable, static scalar-wave perturbations of a Schwarzschild black hole. [More precisely, show that all static solutions to equation (32.27b) become infinite at either the horizon (r=2M,r^(**)=-oo)\left(r=2 M, r^{*}=-\infty\right) or at radial infinity.] This suggests that somehow the black hole formed by collapse must divest itself of the star's external scalar field before it can settle down into a quiescent state.
(c) The general solution to the wave equation (32.27b) can be written in terms of a Fourier transform. For waves that begin near the horizon, propagate outward, and are partially transmitted and partially reflected ("rightward-propagating waves"), show that the general solution is
{:(32.28a)Psi_(ℓ)(t,r^(**))=int_(-oo)^(oo)A(k)R_(k)^(ℓ)(r^(**))e^(-ikt)dk:}\begin{equation*}
\Psi_{\ell}\left(t, r^{*}\right)=\int_{-\infty}^{\infty} A(k) R_{k}^{\ell}\left(r^{*}\right) e^{-i k t} d k \tag{32.28a}
\end{equation*}
where
{:[(32.28b)d^(2)R_(k)^(t)//dr^(**2)=[-k^(2)+V_(efr)(r)]R_(k)^(t)","],[R_(k)^(t)=e^(ikr^(**))+Gamma_(k)^((R))e^(-ikr^(**))" as "r^(**)longrightarrow-oo","],[(32.28c)R_(k)^(ℓ)=T_(k)^((R))e^(ikr^(**))" as "r^(**)longrightarrow oo]:}\begin{gather*}
d^{2} R_{k}^{t} / d r^{* 2}=\left[-k^{2}+V_{\mathrm{efr}}(r)\right] R_{k}^{t}, \tag{32.28b}\\
R_{k}^{t}=e^{i k r^{*}}+\Gamma_{k}^{(R)} e^{-i k r^{*}} \text { as } r^{*} \longrightarrow-\infty, \\
R_{k}^{\ell}=T_{k}^{(R)} e^{i k r^{*}} \text { as } r^{*} \longrightarrow \infty \tag{32.28c}
\end{gather*}
Show that the "reflection and transmission coefficients for rightward-propagating waves," Gamma_(k)^((R))\Gamma_{k}^{(R)} and T_(k)^((R))T_{k}^{(R)}, have the following asymptotic forms for |k|≪1//M|k| \ll 1 / M (short wave number; long wavelength):
where alpha\alpha and beta\beta are constants of order unity. Give a similar analysis for waves that impinge on a Schwarzschild black hole from outside ("leftward-propagating waves").
(d) Show that, as the star collapses into the horizon, the world line of its surface in (t,r^(**))\left(t, r^{*}\right) coordinates is
where R_(0)^(**)R_{0}{ }^{*} is related to the magnitude aa of the surface's 4 -acceleration ( a > 0a>0 for outward 4 -acceleration) by
{:(32.29b)R_(0)^(**)=(8M//e){1+16 Ma[Ma+(M^(2)a^(2)+(1)/(8))^(1//2)]}.:}\begin{equation*}
R_{0}^{*}=(8 M / e)\left\{1+16 M a\left[M a+\left(M^{2} a^{2}+\frac{1}{8}\right)^{1 / 2}\right]\right\} . \tag{32.29b}
\end{equation*}
Thus, the world line of the surface appears to become null near the horizon (t+r^(**)-=( widetilde(V))=:}\left(t+r^{*} \equiv \widetilde{V}=\right. constant); of course, this is due to pathology of the coordinate system there. Show, further, that the scalar field on the star's surface ( vec(V)=(\vec{V}= constant )) must vary as
when the star is approaching the horizon (t longrightarrow oo,r^(**)longrightarrow-oo,( widetilde(U))longrightarrow oo)\left(t \longrightarrow \infty, r^{*} \longrightarrow-\infty, \widetilde{U} \longrightarrow \infty\right), in order that the rate of change of Psi_(ℓ)\Psi_{\ell} be finite as measured on the star's surface. Notice that Q_(10)Q_{10} is the "final value" of the scalar field on the star's surface. It can be regarded as an outgoing wave with zero wave number (infinite wavelength); and, consequently, it gets completely and
destructively reflected by the effective potential [see equation (32.28d); also Box 32.2]. Conclusion: All multipoles of the scalar field die out at finite r^(**)r^{*} as t longrightarrow oot \longrightarrow \infty. (Price's theorem for a scalar field.) For a more detailed analysis, including the rates at which the multipoles die out, see Price (1971, 1972a) or Thorne (1972).
Exercise 32.11. NEWMAN-PENROSE "CONSTANTS"
[See Press and Bardeen (1971), Bardeen and Press (1972), and Piir (1971) for more details than are presented here.]
Wheeler (1955) showed that Maxwell's equations for an ℓ\ell-pole electromagnetic field residing in the Schwarzschild geometry can be reduced to the wave equation
[electromagnetic analogue of ( 32.27 b)) ]. After this equation has been solved, the components of the electromagnetic field can be obtained by applying certain differential operators to Psi_(l)(t,r^(**))Y_(ℓm)(theta,phi)\Psi_{l}\left(t, r^{*}\right) Y_{\ell m}(\theta, \phi).
(a) Show that the general solution to the electromagnetic wave equation (32.30) for dipole (ℓ=1)(\ell=1) fields, with outgoing-wave boundary conditions at r^(**)longrightarrow+oor^{*} \longrightarrow+\infty, has the form
{:[ widetilde(U)=t-r^(**)" is "retarded time", and "],[(32.31b)f_(1)^(')=f_(0)","quadf_(2)^(')=0","quad dots","quadf_(n)^(')=-((n+1)(n-2))/(2n)f_(n-1)+(n-2)Mf_(n-2)]:}\begin{gather*}
\widetilde{U}=t-r^{*} \text { is "retarded time", and } \\
f_{1}^{\prime}=f_{0}, \quad f_{2}^{\prime}=0, \quad \ldots, \quad f_{n}^{\prime}=-\frac{(n+1)(n-2)}{2 n} f_{n-1}+(n-2) M f_{n-2} \tag{32.31b}
\end{gather*}
When spacetime is flat (M=0)(M=0), this solution becomes
{:(32.31~N)Psi_(1)=f_(1)^(')( widetilde(U))+f_(1)( widetilde(U))//r:}\begin{equation*}
\Psi_{1}=f_{1}^{\prime}(\widetilde{U})+f_{1}(\widetilde{U}) / r \tag{32.31~N}
\end{equation*}
[The 1//r1 / r fall-off of the radiation field f_(1)^(')( widetilde(U))f_{1}^{\prime}(\widetilde{U}) has been factored out of Psi_(1)\Psi_{1}; see the scalar-wave function (32.27a).] The terms f_(2)( widehat(U))//r^(2)+dotsf_{2}(\widehat{U}) / r^{2}+\ldots, which are absent in flat spacetime, are attributable to backscatter of the outgoing waves by the curvature of spacetime. They are sometimes called the "tail" of the waves.
(b) Show that the general static dipole field has the form (32.31a) with
(c) Consider a star (not a black hole!) with a dipole field that is initially static. At time t=0t=0, let the star suddenly change its dipole moment to a new static value D^(')D^{\prime}. Equations (32.31b) demand that f_(2)f_{2} be conserved ["Newman-Penrose (1965) constant"]. Hence, f_(2)f_{2} will always exhibit a value, _(2)^(3)MD{ }_{2}^{3} M D, corresponding to the old dipole moment; it can never change to _(2)^(3)MD^('){ }_{2}^{3} M D^{\prime}. This is a manifestation of the tail of the waves that are generated by the sudden change in dipole moment. To understand this tail effect more clearly, and to discover an important flaw in the above result, evaluate the solution (32.31) for retarded time widetilde(U) > 0\widetilde{U}>0, using the assumptions
(1) field has static form (32.32) for widetilde(U) < 0\widetilde{U}<0,
(2) f_(1)=D^(')f_{1}=D^{\prime} for U^(⏜) > 0\overparen{U}>0.
Put the answer in the form Psi_(1)=(D^('))/(r)+((3)/(2)MD)/(r^(2))+sum_(n=3)^(oo)(2M(D^(')-D)(-1)^(n+1)(n+1) widetilde(U)^(n-2))/((2r)^(n))+O((M^(2))/(r^(3)),(M^(2)( widetilde(U)))/(r^(4)))\Psi_{1}=\frac{D^{\prime}}{r}+\frac{\frac{3}{2} M D}{r^{2}}+\sum_{n=3}^{\infty} \frac{2 M\left(D^{\prime}-D\right)(-1)^{n+1}(n+1) \widetilde{U}^{n-2}}{(2 r)^{n}}+O\left(\frac{M^{2}}{r^{3}}, \frac{M^{2} \widetilde{U}}{r^{4}}\right).
(The terms neglected are small compared to those kept for all widetilde(U)//r\widetilde{U} / r, so long as r≫r \gg.) Evidently, so long as the series converges the Newman-Penrose "constant" (coefficient of 1//r^(2)1 / r^{2} ) remembers the old DD value and is conserved, as claimed above. Show, however, that the series diverges for widetilde(U) > 2r\widetilde{U}>2 r-i.e., it diverges inside a sphere that moves outward with asymptotically (1)/(3)\frac{1}{3} the speed of light. Thus, the Newman-Penrose "constant" is well-defined and conserved only outside the " (1)/(3)\frac{1}{3}-speed-of-light cone."
(d) Sum the series in (32.34) to obtain a solution valid even for widehat(U) > 2r\widehat{U}>2 r :
=" the series "(32.34)" for " widetilde(U) < 2r" (domain of convergence of that series) "=\text { the series }(32.34) \text { for } \widetilde{U}<2 r \text { (domain of convergence of that series) }
=(D^('))/(r)+(3)/(2)(MD^('))/(r^(2))+O((M)/(( widetilde(U))r),(M^(2))/(r^(3)))quad" for " widetilde(U)≫r≫M=\frac{D^{\prime}}{r}+\frac{3}{2} \frac{M D^{\prime}}{r^{2}}+O\left(\frac{M}{\widetilde{U} r}, \frac{M^{2}}{r^{3}}\right) \quad \text { for } \widetilde{U} \gg r \gg M
From this result conclude that at fixed rr and late times the electromagnetic field becomes asymptotically static, and the Newman-Penrose "constant" assumes the new value _(2)^(3)MD{ }_{2}^{3} M D ' appropriate to the new dipole moment.
сниетег 33
BLACK HOLES
A luminous star, of the same density as the Earth, and whose diameter should be two hundred and fifty times larger than that of the Sun, would not, in consequence of its attraction, allow any of its rays to arrive at us; it is therefore possible that the largest luminous bodies in the universe may, through this cause, be invisible.
P. S. LAPLACE (1798)
A dialog explaining why black holes deserve their name
§33.1. WHY "BLACK HOLE"?
Sagredus. What is all this talk about "black holes"? When an external observer watches a star collapse, he sees it implode with ever-increasing speed, until the relativistic stage is reached. Then it appears to slow down and become "frozen," just outside its horizon (gravitational radius). However long the observer waits, he never sees the star proceed further. How can one reasonably give the name "black hole" to such a frozen object, which never disappears from sight?
Salvatius. Let us take the name "black hole" apart. Consider first the blackness. Surely nothing can be blacker than a black hole. The very redshift that makes the collapsing star appear to freeze also makes it darken and become black. In the continuum approximation, where one ignores the discreteness of photons, the intensity of the radiation received by distant observers decreases exponentially as time passes, L prop exp(-t//3sqrt3M)L \propto \exp (-t / 3 \sqrt{3} M), with an exceedingly short ee-folding time
Within a fraction of a second, the star is essentially black. Discreteness of photons makes it even blacker. The number of photons emitted before the star crosses its horizon is finite, so the exponential decay cannot continue
For a more detailed exposition of the foundations of "black-hole physics," see DeWitt and DeWitt (1973).
forever. Eventually-only 10^(-3)(M//M_(o.))10^{-3}\left(M / M_{\odot}\right) seconds after the star begins to dim\operatorname{dim} (see exercise 32.2)-the last photon that will ever get out reaches the distant observers. Thereafter nothing emerges. The star is not merely "essentially black"; it is "absolutely black."
Sagredus. Agreed. But it is the word "hole" that concerns me, not "black." How can one possibly regard the name "hole" as appropriate for an object that hovers forever just outside its horizon. True, absence of light makes the object invisible. But couldn't one always see it by shining a flashlight onto its surface? And couldn't one always fly down to its surface in a rocket ship and scoop up a few of the star's baryons? After all, as seen from outside the baryons at its surface will never, never, never manage to fall into the horizon!
Salvatius. Your argument sounds persuasive. To test its validity, examine the collapse of a spherically symmetric system, using the ingoing Eddington-Finkelstein diagram of Figure 33.1. Let a family of external observers shine their flashlights onto the star's surface, as you have suggested. Let the surface of the star be carefully silvered so it reflects back all light that reaches it. Initially (low down in the spacetime diagram of Figure 33.1) the ingoing light beams
Figure 33.1.
Spherical gravitational collapse of a star to form a black hole, as viewed in ingoing Eddington-Finkelstein coordinates. The "surface of last influence," R\mathscr{R}, is an ingoing null surface that intersects the horizon in coincidence with the surface of the collapsing star. After an external observer, moving forward in time, has passed through the surface of last influence, he cannot interact with and influence the star before it plunges through the horizon. Thus, one can think of the surface of last influence as the "birthpoint" of the black hole. Before passing through this surface, the external observer can say his flashlight is probing the shape of a collapsing star; afterwards, he can regard his signals as probes of a black hole. For further discussion, see text.
reach the star's surface and are reflected back to the flashlights with no difficulty. But there is a critical point-an ingoing radial null surface X\mathscr{X} beyond which reflection is impossible. Photons emitted inward along X\mathscr{X} reach the star just as it is passing through its horizon. After reflection these photons fly "outward" along the horizon, remaining forever at r=2Mr=2 M. Other photons, emitted inward after the flashlight has passed through R\mathscr{R}, reach the surface of the star and are reflected only after the star is inside its horizon. Such photons can never return to the shining flashlights. Once inside the horizon, they can never escape. Thus, the total number of photons returned is finite and is subject to the same blackness decay law as is the intrinsic luminosity of the star. Moreover, if the observers do not turn on their flashlights until after they pass through the null surface R\mathscr{R}, they can never receive back any reflected photons! Evidently, flashlights are of no help in seeing the "frozen star."
Sagredus. I cannot escape the logic of your argument. Nevertheless, seeing is not the only means of interacting with the frozen star. I have already suggested swooping down in a rocket ship and stealing a few baryons from its surface. Similarly, one might let matter fall radially inward onto the frozen star. When the matter hits the star's surface, its great kinetic energy of infall will be converted into heat and into outpouring radiation.
Salvatius. Thus it might seem at first sight. But examine again Figure 33.1. No swooping rocket ship and no infalling matter can move inward faster than a light ray. Thus, if the decision to swoop is made after the ship passes through the surface R\mathscr{R}, the rocket ship has no possibility of reaching the star before it plunges through the horizon; the rocket and pilot cannot touch the star, sweep up baryons, and return to tell their tale. Similarly, infalling matter to the future of ℜ\mathscr{\Re} can never hit the star's surface before passing through the horizon. The surface R\mathscr{R} is, in effect, a "surface of last influence." Once anybody or anything has passed through R\mathscr{R}, he or it has no possibility whatever of influencing or interacting with the star in any way before it plunges through the horizon. Thus, from a "causal" or "interaction" standpoint, the collapsing star becomes a hole in space at the surface ℜ\mathscr{\Re}. This hole is not black at first. Radiation from the collapsing star still emerges after R\mathscr{R} because of finite light-propagation times, just as radiation still reaches us today from the hot big-bang explosion of the universe. But if an observer at radius r≫2Mr \gg 2 M waits for rays emitted at R∼3MR \sim 3 M to get back to observer), then he will see the newly formed hole begin to turn black; and within a time Delta t∼(10^(-3):}\Delta t \sim\left(10^{-3}\right. seconds )(M//M_(o.)))\left(M / M_{\odot}\right) thereafter, it will be completely black.
Sagredus. You have convinced me. For all practical purposes the phrase "black hole" is an excellent description. The alternative phrases "frozen star" and "collapsed star," which I find in the pre-1969 physics literature, emphasize an optical-illusion aspect of the phenomenon. Attention must be directed away from the star that created the black hole, because beyond the surface of last influence one has no means to interact with that star. The star is irrelevant
to the subsequent physics and astrophysics. Only the horizon and its external spacetime geometry are relevant for the future. Let us agree to call that horizon the "surface of a black hole," and its external geometry the "gravitational field of the black hole."
Salvatius. Agreed.
§33.2. THE GRAVITATIONAL AND ELECTROMAGNETIC FIELDS OF A BLACK HOLE
The collapse of an electrically neutral star endowed with spherical symmetry produces a spherical black hole with external gravitational field described by the Schwarzschild line element
{:(33.1)ds^(2)=-(1-2M//r)dt^(2)+(dr^(2))/(1-2M//r)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2)):}\begin{equation*}
d s^{2}=-(1-2 M / r) d t^{2}+\frac{d r^{2}}{1-2 M / r}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right) \tag{33.1}
\end{equation*}
The surface of the black hole, i.e., the horizon, is located at r=2M=r=2 M= (gravitational radius). Only the region on and outside the black hole's surface, r >= 2Mr \geq 2 M, is relevant to external observers. Events inside the horizon can never influence the exterior.
The gravitational collapse of a realistic star (nonspherical, collapse with small but nonzero net charge of one sign or the other) produces a black hole somewhat different from the simple Schwarzschild hole. For collapse with small charge and small asymmetries, perturbation-theory calculations (Box 32.2) predict a final black hole with external field determined entirely by the mass MM, charge QQ, and intrinsic angular momentum SS of the collapsing star. For fully relativistic collapse, with large asymmetries and possibly a large charge, the final black hole (if one forms) is also characterized uniquely by M,QM, Q, and SS. This is the conclusion that strongly suggests itself in 1972 from a set of powerful theorems described in Box 33.1.
Why M,QM, Q, and SS should be the complete governors of the final external field of the black hole, one can understand heuristically as follows. Of all quantities intrinsic to any isolated source of gravity and electromagnetism, only M,QM, Q, and SS possess (and are defined in terms of) unique, conserved imprints in the distant external fields of the source (conserved Gaussian flux integrals; see Box 19.1 and §20.2). When a star collapses to form a black hole, its distant external fields are forced to maintain unchanged the imprints of M,QM, Q, and SS. In effect, M,QM, Q, and SS provide anchors or constraints on the forms of the fields. Initially other constraints are produced by the distributions of mass, momentum, stress, charge, and current inside the star. But ultimately the star plunges through a horizon, cutting itself off causally from the external universe. (The nonpropagation of long-wavelength waves through curved spacetime plays a key role in this cutoff; see Box 32.2.) Subsequently, the only anchors remaining for the external fields are the conserved imprints of MM, QQ, and SS. Consequently, the external fields quickly settle down into unique shapes corresponding to the given M,QM, Q, and SS. Of course, the settling down involves dynamic changes of the fields and an associated outflow of gravitational and electro-
The structure of a black hole is determined uniquely by its mass MM, charge QQ, and intrinsic angular momentum, SS
Heuristic explanation of the M-QM-Q-S uniqueness
Box 33.1 A BLACK HOLE HAS NO "HAIR"
The following theorems come close to proving that the external gravitational and electromagnetic fields of a stationary black hole (a black hole that has settled down into its "final" state) are determined uniquely by the hole's mass MM, charge QQ, and intrinsic angular momentum SS-i.e., the black hole can have no "hair" (no other independent characteristics). For a detailed review, see Carter (1973).
I. Stephen Hawking (1971b, 1972a): A stationary black hole must have a horizon with spherical topology; and it must either be static (zero angular momentum), or axially symmetric, or both.
II. Werner Israel (1967a, 1968): Any static black hole with event horizon of spherical topology has external fields determined uniquely by its mass MM and charge QQ; moreover, those external fields are the Schwarzschild solution if Q=0Q=0, and the Reissner-Nordstrøm solution (exercises 31.8 and 32.1) if Q!=0Q \neq 0 (both special cases of Kerr-Newman; see §33.2).
III. Brandon Carter (1970): "All uncharged, stationary, axially symmetric black holes with event horizons of spherical topology fall into disjoint families not deformable into each other. The black holes in each family have external gravitational fields determined uniquely by two parameters: the mass MM and the angular momentum SS." (Note: the "Kerr solutions"-i.e., "Kerr-Newman" with Q=0Q=0 -form one such family; it is very likely that there are no others, but this has not been proved as of December 1972. It is also likely that Carter's theorem can be extended to the case with charge; but this has also not yet been done.)
IV. Conclusions made by combining all three theorems:
(a) All stationary black holes are axially symmetric.
(b) All static (nonrotating) black holes are characterized uniquely by MM and QQ, and have the Reissner-Nordstrom form.
(c) All uncharged, rotating black holes fall into distinct and disjoint families, with each black hole in a given family characterized uniquely by MM and SS. The Kerr solutions form one such family. There may well be no other family.
V. Remarks and Caveats:
(a) The above statements of the theorems are all somewhat heuristic. Each theorem makes several highly technical assumptions, not stated here, about the global properties of spacetime. These assumptions seem physically reasonable and innocuous, but they might not be.
(b) Progress in black-hole physics is so rapid that, by the time this book is published, there may well exist theorems more powerful than the above, which really prove that "a black hole has no hair."
(c) For insight into the techniques of "global geometry" used in proving the above theorems and others like them, see Chapter 34 ; for greater detail see the forthcoming book by Hawking and Ellis (1973).
(d) For analyses which show that a black hole cannot exert any weak-interaction forces caused by the leptons which have gone down it, see Hartle (1971,1972)(1971,1972) and Teitelboim (1972b,c). For similar analyses which show absence of strong-interaction forces from baryons that have gone down the hole, see Bekenstein (1972a,b) and Teitelboim (1972a).
magnetic waves. And, of course, the outflowing waves carry off mass and angular momentum (but not charge), thereby leaving MM and SS changed. And, of course, the external fields must then readjust themselves to the new MM and SS. But the process will quickly converge, producing a black hole with specific final values of M,QM, Q, and SS and with external fields determined uniquely by those values.
The problem of calculating the external fields for given M,QM, Q, and SS and their given imprints, is analogous to the problem of Plateau-to calculate the shape of a soap film anchored to a wire of given shape.* One calculates the shape of the soap film by seeking a surface of minimum area spanning the bent wire. The condition of minimum area leads to a differential equation describing the soap film, which must be solved subject to the constraint imposed by the shape of the wire.
To calculate the external fields of a black hole, one can extremize the "action integral" int(R+E)sqrt(-g)d^(4)x\int(\mathscr{R}+\mathcal{E}) \sqrt{-g} d^{4} x for interacting gravitational and electromagnetic fields (see Chapter 21) subject to the anchored-down imprints of M,QM, Q, and SS at radial infinity, and subject to the existence of a physically nonsingular horizon (no infinite curvature at horizon!). Extremizing the action is equivalent to solving the coupled Einstein-Maxwell field equations subject to the constraints imprinted by M,QM, Q, and SS, and the existence of the horizon. The derivation of the solution and the proof of its uniqueness are much too complex to be given here. (See references cited in Box 33.1.) However, the solution turns out to be the "Kerr-Newman geometry" and its associated electromagnetic field. †\dagger
Written in the t,r,theta,phit, r, \theta, \phi coordinates of Boyer and Lindquist (1967) (generalization of Schwarzschild coordinates), the Kerr-Newman geometry has the form
{:[ds^(2)=-(Delta)/(rho^(2))[dt-asin^(2)theta d phi]^(2)+(sin^(2)theta)/(rho^(2))[(r^(2)+a^(2))d phi-adt]^(2)],[(33.2)+(rho^(2))/(Delta)dr^(2)+rho^(2)dtheta^(2)","]:}\begin{align*}
d s^{2}= & -\frac{\Delta}{\rho^{2}}\left[d t-a \sin ^{2} \theta d \phi\right]^{2}+\frac{\sin ^{2} \theta}{\rho^{2}}\left[\left(r^{2}+a^{2}\right) d \phi-a d t\right]^{2} \\
& +\frac{\rho^{2}}{\Delta} d r^{2}+\rho^{2} d \theta^{2}, \tag{33.2}
\end{align*}
{:(33.4)a-=S//M-=" angular momentum per unit mass. ":}\begin{equation*}
a \equiv S / M \equiv \text { angular momentum per unit mass. } \tag{33.4}
\end{equation*}
The corresponding electromagnetic field tensor, written as a 2 -form (recall: dx^(alpha)^^\boldsymbol{d} x^{\alpha} \wedge{:dx^(beta)-=dx^(alpha)ox dx^(beta)-dx^(beta)ox dx^(alpha))\left.\boldsymbol{d} x^{\beta} \equiv \boldsymbol{d} x^{\alpha} \otimes \boldsymbol{d} x^{\beta}-\boldsymbol{d} x^{\beta} \otimes \boldsymbol{d} x^{\alpha}\right) is
{:[(33.5)F=Qrho^(-4)(r^(2)-a^(2)cos^(2)theta)dr^^[dt-asin^(2)theta d phi]],[+2Qrho^(-4)ar cos theta sin theta d theta^^[(r^(2)+a^(2))d phi-adt].]:}\begin{align*}
\boldsymbol{F}= & Q \rho^{-4}\left(r^{2}-a^{2} \cos ^{2} \theta\right) \boldsymbol{d} r \wedge\left[\boldsymbol{d} t-a \sin ^{2} \theta \boldsymbol{d} \phi\right] \tag{33.5}\\
& +2 Q \rho^{-4} a r \cos \theta \sin \theta \boldsymbol{d} \theta \wedge\left[\left(r^{2}+a^{2}\right) \boldsymbol{d} \phi-a \boldsymbol{d} t\right] .
\end{align*}
Variational principle for black-hole structure
Details of black-hole structure:
(1) metric ("Kerr-Newman geometry")
(2) electromagnetic field
Expressions (33.2) for the metric and (33.5) for the electromagnetic field are sufficiently long to be somewhat frightening. Therefore, it is helpful to develop some qualitative insight into them and into their implications before attempting detailed computations with them. Boxes 33.2, 33.3, and 33.4 develop qualitative insight by presenting, without derivation, a summary of the key features of the Kerr-Newman geometry and a summary of the physics and astrophysics of black holes. The remainder of this chapter is a Track-2 justification and derivation of some, but not all, of the results cited in Boxes 33.2-33.4.
Box 33.2 KERR-NEWMAN GEOMETRY AND ELECTROMAGNETIC FIELD
I. Equations for metric and electromagnetic field
A. Parameters appearing in equations: M=M= mass, Q=Q= charge, a-=S//M=a \equiv S / M= angular momentum per unit mass, all as measured by their standard imprints on the distant fields.
B. Constraint on parameters:
The Kerr-Newman geometry has a horizon, and therefore describes a black hole, if and only if M^(2) >= Q^(2)+a^(2)M^{2} \geq Q^{2}+a^{2}. It seems likely that in any collapsing body which violates this constraint, centrifugal forces and/or electrostatic repulsion will halt the collapse before a size ∼M\sim M is reached; see equation (33.56).
C. Limiting cases:
Q=0Q=0,
Kerr (1963) geometry;
S=0S=0,
Reissner-Nordstrom geometry and electromagnetic field
(exercises 31.8 and 32.1);
Reissner-Nordstrom geometry and electromagnetic field
(exercises 31.8 and 32.1);| Reissner-Nordstrom geometry and electromagnetic field |
| :--- |
| (exercises 31.8 and 32.1); |
Q=S=0Q=S=0,
Schwarzschild geometry;
M^(2)=Q^(2)+a^(2)M^{2}=Q^{2}+a^{2}
"Extreme Kerr-Newman geometry."
Q=0, Kerr (1963) geometry;
S=0, "Reissner-Nordstrom geometry and electromagnetic field
(exercises 31.8 and 32.1);"
Q=S=0, Schwarzschild geometry;
M^(2)=Q^(2)+a^(2) "Extreme Kerr-Newman geometry."| $Q=0$, | Kerr (1963) geometry; |
| :--- | :--- |
| $S=0$, | Reissner-Nordstrom geometry and electromagnetic field <br> (exercises 31.8 and 32.1); |
| $Q=S=0$, | Schwarzschild geometry; |
| $M^{2}=Q^{2}+a^{2}$ | "Extreme Kerr-Newman geometry." |
D. Boyer-Lindquist (1967) coordinates (t,r,theta,phi-(t, r, \theta, \phi- generalization of Schwarzschild coordinates; black hole rotates in phi\phi direction):
{:[(1)ds^(2)=-(Delta//rho^(2))[dt-asin^(2)theta d phi]^(2)+(sin^(2)theta//rho^(2))[(r^(2)+a^(2))d phi-adt]^(2)],[+(rho^(2)//Delta)dr^(2)+rho^(2)dtheta^(2);],[(2)Delta-=r^(2)-2Mr+a^(2)+Q^(2)","quadrho^(2)-=r^(2)+a^(2)cos^(2)theta.],[F=],[(3)quad Qrho^(-4)(r^(2)-a^(2)cos^(2)theta)dr^^[dt-asin^(2)theta d phi]],[+2Qrho^(-4)ar cos theta sin theta d theta^^[(r^(2)+a^(2))d phi-adt].]:}\begin{align*}
& d s^{2}=-\left(\Delta / \rho^{2}\right)\left[d t-a \sin ^{2} \theta d \phi\right]^{2}+\left(\sin ^{2} \theta / \rho^{2}\right)\left[\left(r^{2}+a^{2}\right) d \phi-a d t\right]^{2} \tag{1}\\
&+\left(\rho^{2} / \Delta\right) d r^{2}+\rho^{2} d \theta^{2} ; \\
& \Delta \equiv r^{2}-2 M r+a^{2}+Q^{2}, \quad \rho^{2} \equiv r^{2}+a^{2} \cos ^{2} \theta . \tag{2}\\
& \boldsymbol{F}= \\
& \quad Q \rho^{-4}\left(r^{2}-a^{2} \cos ^{2} \theta\right) \boldsymbol{d} r \wedge\left[\boldsymbol{d} t-a \sin ^{2} \theta \boldsymbol{d} \phi\right] \tag{3}\\
&+2 Q \rho^{-4} a r \cos \theta \sin \theta \boldsymbol{d} \theta \wedge\left[\left(r^{2}+a^{2}\right) \boldsymbol{d} \phi-a \boldsymbol{d} t\right] .
\end{align*}
E. Kerr coordinates [ widetilde(V),r,theta, widetilde(phi)[\widetilde{V}, r, \theta, \widetilde{\phi}-generalization of ingoing Eddington-Finkelstein coordinates; ( widetilde(V),theta, widetilde(varphi))=(\widetilde{V}, \theta, \widetilde{\varphi})= constant is an ingoing, "radial," null geodesic; black hole rotates in widetilde(phi)\widetilde{\phi} direction]:
Relationship to Boyer-Lindquist:
{:[d widetilde(V)=dt+(r^(2)+a^(2))(dr//Delta)","],[(4)d widetilde(phi)=d phi+a(dr//Delta)],[ds^(2)=-[1-rho^(-2)(2Mr-Q^(2))]d widetilde(V)^(2)+2drd widetilde(V)+rho^(2)dtheta^(2)],[(5)+rho^(-2)[(r^(2)+a^(2))^(2)-Deltaa^(2)sin^(2)theta]sin^(2)theta d widetilde(phi)^(2)-2asin^(2)theta d widetilde(phi)dr],[-2arho^(-2)(2Mr-Q^(2))sin^(2)theta d widetilde(phi)d widetilde(V).],[(6)F=Qrho^(-4)[(r^(2)-a^(2)cos^(2)theta)dr^^d( widetilde(V))-2a^(2)r cos theta sin theta d theta^^d( widetilde(V)):}],[{:-asin^(2)theta(r^(2)-a^(2)cos^(2)theta)dr^^d( widetilde(phi))+2ar(r^(2)+a^(2))cos theta sin theta d theta^^d( widetilde(phi))].]:}\begin{gather*}
\boldsymbol{d} \widetilde{V}=\boldsymbol{d} t+\left(r^{2}+a^{2}\right)(\boldsymbol{d} r / \Delta), \\
\boldsymbol{d} \widetilde{\phi}=\boldsymbol{d} \phi+a(\boldsymbol{d} r / \Delta) \tag{4}\\
d s^{2}=-\left[1-\rho^{-2}\left(2 M r-Q^{2}\right)\right] d \widetilde{V}^{2}+2 d r d \widetilde{V}+\rho^{2} d \theta^{2} \\
+\rho^{-2}\left[\left(r^{2}+a^{2}\right)^{2}-\Delta a^{2} \sin ^{2} \theta\right] \sin ^{2} \theta d \widetilde{\phi}^{2}-2 a \sin ^{2} \theta d \widetilde{\phi} d r \tag{5}\\
-2 a \rho^{-2}\left(2 M r-Q^{2}\right) \sin ^{2} \theta d \widetilde{\phi} d \widetilde{V} . \\
\boldsymbol{F}=Q \rho^{-4}\left[\left(r^{2}-a^{2} \cos ^{2} \theta\right) \boldsymbol{d} r \wedge \boldsymbol{d} \widetilde{V}-2 a^{2} r \cos \theta \sin \theta \boldsymbol{d} \theta \wedge \boldsymbol{d} \widetilde{V}\right. \tag{6}\\
\left.-a \sin ^{2} \theta\left(r^{2}-a^{2} \cos ^{2} \theta\right) \boldsymbol{d} r \wedge \boldsymbol{d} \widetilde{\phi}+2 a r\left(r^{2}+a^{2}\right) \cos \theta \sin \theta \boldsymbol{d} \theta \wedge \boldsymbol{d} \widetilde{\phi}\right] .
\end{gather*}
II. Properties of spacetime geometry
A. Symmetries (§33.4):
The metric coefficients in Boyer-Lindquist coordinates are independent of tt and phi\phi, and in Kerr coordinates are independent of widetilde(V)\widetilde{V} and widetilde(phi)\widetilde{\phi}. Thus the spacetime geometry is "time-independent" (stationary) and axially symmetric. The "Killing vectors" ( $25.2\$ 25.2 ) associated with these two symmetries are (del//del t)_(r,theta,phi)=(del//del widetilde(V))_(r,theta, tilde(phi))(\partial / \partial t)_{r, \theta, \phi}=(\partial / \partial \widetilde{V})_{r, \theta, \tilde{\phi}} and (del//del phi)_(t,r,theta)=(del//del widetilde(varphi))_( tilde(V),r,theta)(\partial / \partial \phi)_{t, r, \theta}=(\partial / \partial \widetilde{\varphi})_{\tilde{V}, r, \theta}.
B. Dragging of inertial frames and static limit (§33.4):
The "dragging of inertial frames" by the black hole's angular momentum produces a precession of gyroscopes relative to distant stars. By this precession one defines and measures the angular momentum of the black hole (see §$19.2\S \$ 19.2§ and 19.3).
The dragging becomes more and more extreme the nearer one approaches the horizon of the black hole. Before the horizon is reached, at a surface described by
the dragging becomes so extreme that no observer can possibly remain at rest there (i.e., be "static") relative to the distant stars. At and inside this surface
(called the "static limit"), all observers with fixed rr and theta\theta must orbit the black hole in the same direction in which the hole rotates:
{:[Omega-=d phi//dt],[ > (a sin theta-sqrtDelta)/((r^(2)+a^(2))sin theta-sqrtDeltaasin^(2)theta)],[( >= 0" for "a=S//M > 0" and "r <= r_(0)).]:}\begin{aligned}
& \Omega \equiv d \phi / d t \\
&>\frac{a \sin \theta-\sqrt{\Delta}}{\left(r^{2}+a^{2}\right) \sin \theta-\sqrt{\Delta} a \sin ^{2} \theta} \\
&\left(\geq 0 \text { for } a=S / M>0 \text { and } r \leq r_{0}\right) .
\end{aligned}
No matter how hard an observer, at fixed (r,theta)(r, \theta) inside the static limit, blasts his rocket engines, he can never halt his angular motion relative to the distant stars.
3. The mathematical foundation for the above statement is this: world lines of the form (r,theta,phi)=(r, \theta, \phi)= constant [tangent vector prop del//del t=\propto \partial / \partial t= "Killing vector in time direction"] change from being timelike outside the static limit to being spacelike inside it. Therefore, on and inside the static limit, no observer can remain at rest.
C. Horizon (§33.4):
As with the Schwarzschild horizon of a nonrotating black hole, so also here, particles and photons can fall inward through the horizon; but no particle or
Box 33.2 (continued)
photon can emerge outward through it.
3. The horizon is "generated" by outgoing null geodesics (outgoing photon world lines).
D. Ergosphere (§33.4):
The "ergosphere" is the region of spacetime between the horizon and the static limit. It plays a fundamental role in the physics of black holes (Box 33.3; §33.7).
The static limit and the horizon touch at the point where they are cut by the axis of rotation of the black hole (theta=0,pi)(\theta=0, \pi); they are well-separated elsewhere with the static limit outside the horizon, unless a=0a=0 (no rotation). When a=0a=0, the static limit and horizon coincide; there is no dragging of inertial frames; there is no ergosphere.
Qualitative representation of horizon, ergosphere, and static limit [adapted from Ruffini and Wheeler (1971b)].
E. Singularity in Boyer-Lindquist coordinates:
For a nonrotating black hole, the Schwarzschild coordinates become singular at the horizon. One manifestation
of the singularity is the infinite amount of coordinate time required for any particle or photon to fall inward through the horizon, t longrightarrow oot \longrightarrow \infty as r longrightarrow2Mr \longrightarrow 2 M. One way to remove the singularity (Eddington-Finkelstein way) is to replace tt by a null coordinate
widetilde(V)=t+r+2M ln |r//2M-1|\widetilde{V}=t+r+2 M \ln |r / 2 M-1|
attached to infalling photons [so (del//del r)_( tilde(V),theta,phi)(\partial / \partial r)_{\tilde{V}, \theta, \phi} is vector tangent to photon world lines].
2. For a rotating black hole, the BoyerLindquist coordinates, being generalizations of the Schwarzschild coordinates, are also singular at the horizon. It requires an infinite coordinate time for any particle or photon to fall inward through the horizon, t longrightarrow oot \longrightarrow \infty as r longrightarrowr_(+)r \longrightarrow r_{+}. But that is not all. The dragging of inertial frames forces particles and photons near the horizon to orbit the black hole with Omega-=d phi//dt > 0\Omega \equiv d \phi / d t>0. Consequently, for a particle falling through the horizon ( r longrightarrowr_(+)r \longrightarrow r_{+}), just as t longrightarrow oot \longrightarrow \infty, so also phi longrightarrow oo\phi \longrightarrow \infty (infinite twisting of world lines around horizon).
3. To remove the coordinate singularity, one must perform an infinite compression of coordinate time, and an infinite untwisting in the neighborhood of the horizon. Kerr coordinates achieve this by replacing tt with a null coordinate widetilde(V)\widetilde{V}, and phi\phi with an untwisted angular coordinate widetilde(phi)\widetilde{\phi} :
{:[d widetilde(V)=dt+(r^(2)+a^(2))(dr//Delta)","],[d widetilde(phi)=d phi+a(dr//Delta)]:}\begin{aligned}
& \boldsymbol{d} \widetilde{V}=\boldsymbol{d} t+\left(r^{2}+a^{2}\right)(\boldsymbol{d} r / \Delta), \\
& d \widetilde{\phi}=\boldsymbol{d} \phi+a(\boldsymbol{d} r / \Delta)
\end{aligned}
Both of the new coordinates are attached to the world lines of a particular family of infalling photons; (del//del r)_( tilde(v),theta, tilde(phi))(\partial / \partial r)_{\tilde{v}, \theta, \tilde{\phi}} is the field of vectors tangent to the world lines of this family of photons (ingoing principal null congruence; §33.6).
F. Spacetime diagram:
A spacetime diagram in Kerr coordinates looks much like an EddingtonFinkelstein diagram for the Schwarzschild geometry. In both cases, one plots the surfaces of constant widetilde(V)\widetilde{V} not as horizontal planes, but as "backward light cones" (" 45 -degree surfaces"), because they are generated by the world lines of ingoing photons. Equivalently, one plots surfaces of constant widetilde(t)-= widehat(V)-r\widetilde{t} \equiv \widehat{V}-r as horizontal planes.
The key differences between a Kerr diagram and an Eddington-Finkelstein diagram are: (a) Because the KerrNewman geometry is not spherical, a Kerr diagram with one rotational degree of freedom suppressed loses information about the geometry. Kerr diagrams are usually made for the equatorial "plane," theta=pi//2\theta=\pi / 2. (b) Just as the horizon pulls the light cones inward, so the dragging of inertial frames tilts the light cones in the direction of increasing widetilde(phi)\widetilde{\phi}, for a > 0a>0 and r=r= constant.
(c) The ingoing edge of a light cone (dr//d widetilde(V)=-oo)(d r / d \widetilde{V}=-\infty) does not tilt toward increasing widetilde(phi)\widetilde{\phi}; the transformation from Boyer-Lindquist coordinates to Kerr coordinates untwists the tilt with decreasing rr, which would otherwise be produced by "frame dragging."
The shapes of the light cones reveal the special features of the static limit and horizon. At the static limit, a vertical world line [r,theta,( widetilde(phi)):}\left[r, \theta, \widetilde{\phi}\right. constant; (del//del widetilde(V))_(r,theta, widetilde(Phi))(\partial / \partial \widetilde{V})_{r, \theta, \widetilde{\Phi}}=(del//del t)_(r,theta,phi)==(\partial / \partial t)_{r, \theta, \phi}= tangent vector] lies on the light cone. At the horizon the light cones tilts fully inward, except for a single line of tangency to the horizon. Notice that the line of tangency has d widehat(phi)//d widetilde(V)=a//(r_(+)^(2)+a^(2))!=0d \widehat{\phi} / d \widetilde{V}=a /\left(r_{+}^{2}+a^{2}\right) \neq 0. Equivalently, the outgoing null geodesics, which generate the horizon, twist about it ("barber-pole-twist")-yet another manifestation of the dragging of inertial frames.
Kerr diagram for equatorial slice ( theta=pi//2\theta=\pi / 2 ) through the spacetime of an "extreme Kerr" black hole (Q=0,a=M)(Q=0, a=M).
View from above showing the shapes of the light cones as a function of radius
Box 33.2 (continued)
The Kerr diagram, like the EddingtonFinkelstein diagram, describes infall through the horizon in a faithful, nonsingular way.
[The term "Kerr diagram" is a misnomer. Kerr has not published such diagrams himself, though nowadays others construct such diagrams using his coordinate system. Penrose is the originator and greatest exploiter of such diagrams (see, e.g., Penrose, 1969). But several other types of diagrams bear Penrose's name, so it would be confusing to name them all after him.]
G. Maximal analytic extension of Kerr-Newman geometry:
When one abstracts the Schwarzschild geometry away from all sources (Chapter 31), one discovers that it describes an expanding and recontracting bridge, connecting two different universes. But in the context of black holes, only half of the Schwarzschild geometry (regions I and II) is relevant. The other half (regions III and IV) gets fully replaced by the interior of the star that collapsed to form the black hole. Because only a
part of the Schwarzschild geometry comes into play, ingoing EddingtonFinkelstein coordinates-which describe I and II well, but III and IV badly-are well-suited to black-hole physics.
Similarly, when one abstracts the KerrNewman geometry away from all sources, one discovers that it describes a much larger, and more complex spacetime manifold than one might ever have suspected. This "maximum analytic extension" of the Kerr-Newman geometry has been analyzed in detail by Boyer and Lindquist (1967) and by Carter (1966a, 1968a). But it is totally irrelevant to the subject of black holes, for two reasons. First, as with Schwarzschild, the star that collapsed to form the black hole replaces most of the inward extension of the Kerr-Newman manifold. Second, even outside the star, the Kerr-Newman geometry does not properly represent the true geometry at early times. At early times the star has not got far down the road to collapse. Gravitational moments of the star arise from mountains or prominences or turbulence or other particularities that have not yet gone into the meat grinder. The geometry departs from flatness (1) by a term that varies for large distances as mass divided by distance, and (2) by another term that varies as angular momentum divided by the square of the distance and multiplied by a spherical harmonic of order one, but also (3) by higher-order terms proportional to higher-order mass moments multiplied by higher spherical harmonics. These higher-order terms normally will deviate at early times from the corresponding terms in the mathematical analysis of the Kerr-Newman geometrythough the deviations will die out as time passes. For a system endowed with spherical symmetry, no such higherorder terms do occur or can occur. Therefore the geometry outside is
Schwarzschild in character at all stages of the collapse. However, when the system lacks spherical symmetry, the geometry outside initially departs from Kerr-Newman character. Only well after the collapse occurs (asymptotic future), and in the region at and outside the horizon, is the Kerr-Newman geometry a faithful descriptor of a black hole. This region is described in a nonsingular manner by Kerr coordinates and Kerr diagrams; and it is the only region that this book will explore.
H. Test-particle orbits
See §§33.5-33.8\S \S 33.5-33.8§§ and Box 33.5 .
III. Properties of electromagnetic field (§33.3):
A. Far from the black hole, where spacetime is nearly flat, in the usual spherical orthonormal frame (omega^( hat(t))=dt,omega^( hat(ı))=dr:}\left(\boldsymbol{\omega}^{\hat{t}}=\boldsymbol{d} t, \boldsymbol{\omega}^{\hat{\imath}}=\boldsymbol{d} r\right., {:omega^( hat(theta))=rd theta,omega^( hat(phi))=r sin theta d phi)\left.\boldsymbol{\omega}^{\hat{\theta}}=r \boldsymbol{d} \theta, \boldsymbol{\omega}^{\hat{\phi}}=r \sin \theta \boldsymbol{d} \phi\right), the electric and magnetic fields have dominant components
These reveal that Q=Q= charge of black hole, R-=Qa=\mathscr{R} \equiv Q a= magnetic dipole moment of black hole.
B. Notice that the gyromagnetic ratio, gamma-=\gamma \equiv (magnetic moment)/(angular momentum), is equal to Q//M=Q / M= (charge/mass), just as for an electron!
C. Notice that the value of the magnetic moment, like all other features of the black hole, is determined uniquely by the hole's mass, charge, and angular momentum: pi=QS//M\mathscr{\pi}=Q S / M. This illustrates the theorem (Box 33.1) that a black hole has no "hair."
D. Other electric and magnetic moments are nonzero, but are determined uniquely by M,SM, S, and QQ.
E. Near the black hole, the curvature of spacetime deforms the electric and magnetic fields produced by the charged, rotating black hole. For a mathematical description of this deformed field, see Cohen and Wald (1971); for a diagrammatic representation, Hanni and Ruffini (1973).
Box 33.3 THE ASTROPHYSICS OF BLACK HOLES
Black holes in nature should participate in astrophysical processes that are as varied as those for stars. By searching for observable phenomena associated with these processes, astronomers have a good chance of discovering the first black hole sometime during the 1970's. This box lists some possible astrophysical processes, and a few relevant references.
I. Mechanisms of Formation
A. "Direct, in isolation": A massive star (M >=(M \geq3M_(o.)3 M_{\odot} ) collapses, almost spherically, producing a collapsed neutron-star core that is too massive to support itself against gravity. Gravity pulls the core on inward, producing a horizon and black hole. [May
Box 33.3 (continued)
and White (1966, 1967); Chapter 32 of this book.]
B. "Indirect, in isolation": "Collapse, pursuit, and plunge scenario" depicted in Figure 24.3 [Ruffini and Wheeler (1971b).]
C. "In the thick of things": Stars collected into a dense cluster (e.g., the nucleus of a galaxy) exchange energy. Some acquire energy and move out into a halo. Others lose energy and make a more compact cluster. This process of segregation continues. The cluster becomes so compact that collisions ensue and gas is driven off. The gas moves toward the center of the gravitational potential well. Out of it new stars form. The process continues. Eventually star-star collisions may become sufficiently energetic and inelastic that the centers of the colliding stars coalesce. In this way supermassive objects may be built up and may evolve. Ultimately (1) many "small" stars may collapse to form "small" black holes (M∼M_(o.));(2)\left(M \sim M_{\odot}\right) ;(2) one or more supermassive stars may collapse to form huge black holes ( M∼10^(4)M_(o.)M \sim 10^{4} M_{\odot} to 10^(9)M_(o.)10^{9} M_{\odot} ); (3) the entire conglomerate of stars and gas and holes may become so dense that it collapses to form a single gigantic hole. [Sanders (1970), Spitzer (1971), LyndenBell (1967, 1969), Colgate (1967), $$24.5\$ \$ 24.5, 24.6, 25.7 of this book.]
D. "Primordially": Perturbations in the initial density distribution of the expanding universe may produce collapse, resulting in "primordial black holes." Those holes would subsequently grow by accretion of radiation and matter. By today all such holes might have grown into enormous objects [M∼10^(17)M_(o.);:}\left[M \sim 10^{17} M_{\odot} ;\right. Zel'dovich and Novikov (1966)]; but some of them might have avoided such growth and might be as small as 10^(-5)10^{-5} grams [Hawking (1971a)].
II. How many black holes are there in our galaxy today?
Peebles (1972) has given an excellent review of this issue and of prospects for finding black holes in the near future. He says "a good fraction of the mass of the disc of our galaxy was deposited [long ago] in stars capable of collapsing to black holes. ... The indication is that the galaxy's disk may contain on the order of 10^(9)10^{9} black holes."
III. "Live" black holes versus "dead" black holes
A. A Schwarzschild black hole is "dead" in the sense that one can never extract from it any of its mass-energy. One aspect of this "deadness"-the fact that a Schwarzschild black hole is stable against small perturbations-is essential (1) to the identification of a black hole with the ultimate "ground state" of a large mass, and (2) to any assertion that general relativity theory predicts the possible existence of black holes. [For a proof of stability see Vishveshwara (1970). The problem was formulated, and most of the necessary techniques developed, by Regge and Wheeler (1957), with essential contributions also by Zerilli (1970a).] Thus a small pulse of gravitational (or other) radiation impinging on a Schwarzschild black hole does not initiate a transition of the black hole into a very different object or state.
B. A Kerr-Newman black hole-which is rotating or charged or both-is not dead. The rotational and electromagnetic contributions to the mass-energy can be extracted. (See §§33.7\S \S 33.7§§ and 33.8 for mathematical details.) Thus, such black holes are "live"; they can inject energy into their surroundings. By a suitable arrangement of external apparatus, one can trigger an exponentially growing energy release [Press and Teukolsky (1972).] But for a perturbed
black hole in isolation, the release is always "controlled" and damped; i.e., Kerr black holes are stable in any classical context [Press and Teukolsky (1973)].
C. Most objects (massive stars; galactic nuclei; ...) that can collapse to form black holes have so much angular momentum that the holes they produce should be "very live" ( aa nearly equal to M;SM ; S nearly equal to M^(2)M^{2} ). [Bardeen (1970a).]
D. By contrast, it is quite probable (but far from certain) that no black hole in the universe has substantial charge-i.e., that all black holes have Q≪MQ \ll M. A black hole with Q∼MQ \sim M (say, Q > 0Q>0 for concreteness) would exert attractive electrostatic forces on electrons, and repulsive electrostatic forces on protons, that are larger than the hole"s gravitational pull by the factor
Here ee is the electron charge and mu\mu is the electron (or proton) mass. Such huge differential forces are likely to pull in enough charge from outside the hole to neutralize it.
E. But one has learned from the "unipolar induction process" for neutron stars [Goldreich and Julian (1968)] that charge neutralization can sometimes be circumvented. Whether any black-hole process can possibly prevent neutralization one does not know in 1972.
IV. Interaction of a black hole with its environment
A. Gravitational pull: A black hole exerts a gravitational pull on surrounding matter and stars. The pull is indistinguishable, at radii r≫Mr \gg M, from the pull of a star with the same mass.
B. Accretion and emission of xx-rays and gamma\gamma rays: Gas surrounding a black hole gets
pulled inward and is heated by adiabatic compression, by shock waves, by turbulence, by viscosity, etc. Before it reaches the horizon, the gas may become so hot that it emits a large flux of xx-rays and perhaps even gamma\gamma-rays. Thus, accreting matter can convert a black hole into a glowing "white" body [for a review of the literature, see Novikov and Thorne (1973)]. Accretion from a nonrotating gas cloud tends to decrease the angular momentum of a black hole [preferential accretion of particles with "negative" angular momentum; Doroshkevich (1966), Godfrey (1970a)]. But the gas surrounding a hole is likely to be rotating in the same direction as the hole itself, and to maintain S∼M^(2)S \sim M^{2} [more precisely, S~~0.998M^(2);S \approx 0.998 M^{2} ; Thorne (1973b)].
C. A lump of matter (an "asteroid" or a "planet" or a star) falling into a black hole should emit a burst of gravitational waves as it falls. The total energy radiated is E∼0.01 mu(mu//M)E \sim 0.01 \mu(\mu / M), where mu\mu is the mass of the object. [Zerilli (1970b); Davis, Ruffini, Press, and Price (1971); Figure 36.2 of this book.]
D. An object in a stable orbit around a black hole should spiral slowly inward because of loss of energy through gravitational radiation, until it reaches the most tightly bound, stable circular orbit. It should then fall quickly into the hole, emitting a "lastgasp burst" of waves. The total energy radiated during the slow inward spiral is equal to the binding energy of the last stable circular orbit:
={[0.0572 mu" for Schwarzschild hole, "],[0.4235 mu" for Kerr hole with "]}quad S=M^(2),Q=0.$=\left\{\begin{array}{l}
0.0572 \mu \text { for Schwarzschild hole, } \\
0.4235 \mu \text { for Kerr hole with }
\end{array}\right\} \quad S=M^{2}, Q=0 . ~ \$
Here mu\mu is the rest mass of the captured object. [Box 33.5.] The total energy in the last-gasp burst is E∼0.01 mu(mu//M)E \sim 0.01 \mu(\mu / M) if mu≪M\mu \ll M. [Fig. 36.2.]
Box 33.3 (continued)
E. When matter falls down a black hole, it can excite the hole's external spacetime geometry into vibration. The vibrations are gradually converted into gravitational waves, some of which escape, others go down the hole. [Press (1971, Goebel (1972).] These vibrations are analogous to an "incipient gravitational geon" [Wheeler (1962); Christodoulou (1971)]-except that for a vibrating black hole the background Kerr geometry holds the vibration energy together (prevents it from propagating away immediately), whereas in a geon it is curvature produced by the "vibration energy" itself that prevents disruption.
F. By a non-Newtonian, induction-zone (i.e., nonradiative) gravitational interaction, a black hole gradually transfers its angular momentum to any non-axially-symmetric, nearby distribution of matter or fields. [Hawking (1972a); Ipser (1971), Press (1972), Hawking and Hartle (1972).]
G. A star or planet falling into a large black hole will get torn apart by tidal gravitational forces. If the tearing occurs near but outside the horizon, it may eject a blob of stellar matter that goes out with relativistic velocity ("tube-of-toothpaste effect"). Moreover, the outgoing jet may extract a substantial amount of rotational energy from the hole's ergosphere-i.e., the hole might throw it off with a rest mass plus kinetic energy in excess of the rest mass of the original infalling object. [Wheeler (1971d); §§33.7 and 33.8.]
H. The magnetic field lines of a charged black hole may be anchored to surrounding plasma, may get wound up as the hole rotates, and may shake, twitch, and excite the plasma.
V. Collisions between black holes
A. Two black holes can collide and coalesce; but there is no way to blast a black hole apart into several black holes [Hawking (1972a); exercise 34.4].
B. When two black holes collide and coalesce, the surface area of the final black hole must exceed the sum of the surface areas of the two initial black holes ("second law of black-hole dynamics"; Hawking (1971a,b); Box 33.4; §34.5). This constraint places an upper limit on the amount of gravitational radiation emitted in the collision. For example, if all three holes are of the Schwarzschild variety and the two initial holes have equal masses M//2M / 2, then
VI. Where and how to search for a black hole [For a detailed review, see Peebles (1971)]:
A. When it forms, by the burst or bursts of gravitational radiation given off during formation [Figure 24.3].
B. In a binary star system: black-hole component optically invisible, but may emit xx-rays and gamma\gamma-rays due to accretion; visible component shows telltale Doppler shifts [Hoyle, Fowler, Burbidge, and Burbidge (1964); Zel'dovich and Guseynov (1965); Trimble and Thorne (1969); Pringle and Rees (1972); Shakura and Sunyaev (1973)]. The velocity of the visible component and the period give information on the mass of the invisible component. If
mass of this invisible component is four solar masses or more, it cannot be an ordinary star, because an ordinary star of that mass would have (4)^(3)=64(4)^{3}=64 times the luminosity of the sun. Neither can it be a white dwarf or a neutron star because either object, so heavy, would instantly collapse to a black hole. Therefore, it is attrac-tive-though not necessarily compelling [see Trimble and Thorne (1969)]-to identify the invisible object as a black hole.
C. [But one must not expect to see any noticeable gravitational lens action from a black hole in a binary system: if it taxed the abilities of astronomers for decades to see the black disc of Mercury, 4,800km4,800 \mathrm{~km} in diameter, swim across the great face of the sun, little hope there is to see a black hole with an effective radius of only ∼3\sim 3 km , enormously more remote, occult a companion star. Significant lens action requires that the lens (black hole) be separated by a normal interstellar distance from the star it focuses; whence the impact parameter of the focused rays is more than a stellar radius, so the lens action is not more than that of a normal star. Moreover,
even with 10^(9)10^{9} black holes in the galaxy, only one per year would pass directly between the Earth and a more distant star, and produce significant lens action (Refsdal, 1964). Chance of watching the right spot on the sky at the right time with a sufficiently strong telescope: nil!]
D. At the center of a globular cluster, where a black hole may settle down, attract normal stars to its vicinity, and thereby produce a cusp in the distribution of light from the cluster. [Cameron and Truran (1971), Peebles (1971).]
E. In the nucleus of a galaxy, including even the Milky Way, where a single huge black hole ( M∼10^(4)M \sim 10^{4} to 10^(8)M_(o.)10^{8} M_{\odot} ) might sit as an end-product of earlier activity of the galactic nucleus. Such a hole will emit gravitational waves, light, and radio waves as it accretes matter. Much of the light may be converted into infrared radiation by surrounding dust. The black hole may also produce jets and other nuclear activity. [Lynden-Bell (1969), Lynden-Bell and Rees (1971), Wheeler (1971d), Peebles (1971).]
Box 33.4 THE LAWS OF BLACK-HOLE DYNAMICS
The black-hole processes described in Box 33.3 are governed by the standard laws of physics: general relativity, plus Maxwell electrodynamics, plus hydrodynamic, quantum mechanical, and other laws for the physics of matter and radiation. From these standard laws of physics, one can derive certain "rules" or "constraints," which all black-hole processes must satisfy. Those rules have a power, elegance, and simplicity that rival and resemble the power, elegance, and simplicity of the laws of thermodynamics. Therefore, they have been given the analogous name "the laws of black-hole dynamics" (Israel 1971). This box states two of the laws of black-hole
Box 33.4 (continued)
dynamics and some of their ramifications. Two additional laws, not discussed here, have been formulated by Bardeen, Carter, and Hawking (1973).
I. The First and Second Laws of Black-Hole Dynamics.
A. The first law.
Like the first law of thermodynamics, the first law of black-hole dynamics is the standard law of conservation of total energy, supplemented by the laws of conservation of total momentum, angular momentum, and charge. For detailed discussions of these conservation laws, see Box 19.1 and Chapter 20.
Specialized to the case where matter falls down a black hole and gravitational waves pour out, the first law takes the form depicted and discussed near the end of Box 19.1.
Specialized to the case of infalling electric charge, the first law says that the total charge QQ of a black hole, as measured by the electric flux emerging from it, changes by an amount equal to the total charge that falls down the hole,
Delta Q=q_("that falls in ")\Delta Q=q_{\text {that falls in }}
Specialized to the case where two black holes collide and coalesce (example given in Box 33.3), the first law says: (a) Let P_(1)\boldsymbol{P}_{1} and P_(2)\boldsymbol{P}_{2} be the 4-momenta of the two black holes as measured gravitationally, when they are so well-separated that they have negligible influence on each other. ( P_(1)\boldsymbol{P}_{1} and P_(2)\boldsymbol{P}_{2} are 4 -vectors in the surrounding asymptotically flat spacetime.) Similarly, let J_(1)\boldsymbol{J}_{1} and J_(2)\boldsymbol{J}_{2} be their total angular-momentum tensors (not intrinsic angular-momentum vectors!) relative to some arbitrarily chosen origin of coordinates, P_(0)\mathscr{P}_{0}, in the surrounding asymptotically flat spacetime ( J_(1)\boldsymbol{J}_{1} and J_(2)\boldsymbol{J}_{2} contain orbital angular momentum, as well as intrinsic angular momentum; see Box 5.6.). (b) Let P_(3)\boldsymbol{P}_{3} and J_(3)\boldsymbol{J}_{3} be the similar total 4-momentum and angular momentum of the final black hole. (c) Let P_(r)\boldsymbol{P}_{r} and J_(r)\boldsymbol{J}_{r} be the total 4-momentum and angular momentum radiated as gravitational waves during the collision and coalescence. Then
[Note: to calculate the mass and intrinsic angular momentum of the final black hole from a knowledge of P_(3)\boldsymbol{P}_{3} and J_(3)\boldsymbol{J}_{3}, follow the prescription of Box 5.6. In that prescription, the world line of the final black hole is that world line, in the distant asymptotically Lorentz coordinates, on which the hole's distant spherical field is centered.
B. The second law [expounded and applied by Hawking (1971b, 1972a)]. When anything falls down a black hole, or when several black holes collide and coalesce or collide and scatter, or in any other process whatsoever involving black holes, the sum of the surface areas (or squares of "irreducible masses"-see equation 3 below) of all black holes involved can never decrease. (See §34.5\S 34.5§ for proof.) This is the second law of black-hole dynamics.
II. Reversible and Irreversible Transformations; Irreducible Mass [Christodoulou (1970); Christodoulou and Ruffini (1971)—results derived independently of and simultaneously with Hawking's discovery of the second law.]
A. Consider a single Kerr-Newman black hole interacting with surrounding matter and fields. Its surface area, at any moment of time, is given in terms of its momentary mass MM, charge QQ, and intrinsic angular momentum per unit mass a-=S//Ma \equiv S / M by
(exercise 33.12). Interaction with matter and fields may change M,QM, Q, and aa in various ways; MM can even be decreased-i.e., energy can be extracted from the black hole! [Penrose (1969); §33.7.] But whatever may be the changes, they can never reduce the surface area AA. Moreover, if any change in M,QM, Q, and aa ever increases the surface area, no future process can ever reduce it back to its initial value.
B. Thus, one can classify black-hole processes into two groups.
Reversible transformations change M,QM, Q, or aa or any set thereof, while leaving the surface area fixed. They can be reversed, bringing the black hole back to its original state.
Irreversible transformations change M,QM, Q, or aa or any set thereof, and increase the surface area in the process. Such a transformation can never be reversed. The black hole can never be brought back to its original state after an irreversible transformation.
C. Examples of reversible transformations and of irreversible transformations induced by infalling particles are presented in §§33.7\S \S 33.7§§ and 33.8 .
D. The reversible extraction of charge and angular momentum from a black hole (decrease in QQ and aa holding AA fixed) necessarily reduces the black hole's mass (energy extraction!). By the time all charge and angular momentum have been removed, the mass has dropped to a final "irreducible value" of
{:(2)M_(ir)=(A//16 pi)^(1//2)=((" mass of Schwarzschild ")/(" black hole of surface area "A)):}\begin{equation*}
M_{\mathrm{ir}}=(A / 16 \pi)^{1 / 2}=\binom{\text { mass of Schwarzschild }}{\text { black hole of surface area } A} \tag{2}
\end{equation*}
Box 33.4 (continued)
E. Expressed in terms of this final, irreducible mass, the initial mass-energy of the black hole (with charge QQ and intrinsic angular momentum SS ) is
[This formula, derived by Christodoulou and Ruffini, may be obtained by combining equations (1), (2), and S=Ma]S=M a].
F. Thus, one can regard the total mass-energy of a black hole as made up of an irreducible mass, an electromagnetic mass-energy, and a rotational energy. But one must resist the temptation to think of these contributions as adding linearly. On the contrary, they combine in a way [equation (3)] analogous to the way rest mass and linear momentum combine to give energy, E^(2)=E^{2}=m^(2)+p^(2)m^{2}+\boldsymbol{p}^{2}.
G. Contours of constant M//M_(ir)M / M_{i r} are depicted below in the "charge-angular momentum plane." Black holes can exist only in the interior of the region depicted ( Q^(2)+a^(2) <= M^(2)Q^{2}+a^{2} \leq M^{2} ). [Diagram adapted from Christodoulou (1971).]
H. Since a black hole's irreducible mass is proportional to the square root of its surface area, one can restate the second law of black-hole dynamics as follows:
In black-hole processes the sum of the squares of the irreducible masses of all black holes involved can never decrease.
§33.3. MASS, ANGULAR MOMENTUM, CHARGE, AND MAGNETIC MOMENT
It is instructive to verify that the constants M,QM, Q, and aa, which appear in equations (33.2)-(33.5) for the Kerr-Newman geometry and electromagnetic field, are actually the black hole's mass, charge, and angular momentum per unit mass, as claimed above.
Mass and angular momentum are defined by their imprints on the spacetime geometry far from the black hole. Therefore, to calculate the mass and angular momentum, one can expand the line element (33.2) in powers of 1//r1 / r and examine the leading terms:
{:[ds^(2)=-[1-(2M)/(r)+O((1)/(r^(2)))]dt^(2)-[(4aM)/(r)sin^(2)theta+O((1)/(r^(2)))]dtd phi],[(33.6)+[1+O((1)/(r))][dr^(2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(˙)^(2))]]:}\begin{align*}
d s^{2}= & -\left[1-\frac{2 M}{r}+O\left(\frac{1}{r^{2}}\right)\right] d t^{2}-\left[\frac{4 a M}{r} \sin ^{2} \theta+O\left(\frac{1}{r^{2}}\right)\right] d t d \phi \\
& +\left[1+O\left(\frac{1}{r}\right)\right]\left[d r^{2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \dot{\phi}^{2}\right)\right] \tag{33.6}
\end{align*}
The examination is facilitated by transforming to asymptotically Lorentz coordi-nates- x=r sin theta cos phi,y=r sin theta sin phi,z=r cos thetax=r \sin \theta \cos \phi, y=r \sin \theta \sin \phi, z=r \cos \theta :
{:[ds^(2)=-[1-(2M)/(r)+O((1)/(r^(2)))]dt^(2)-[(4aM)/(r^(3))+O((1)/(r^(4)))][xdy-ydx]],[(33.6')+[1+O((1)/(r))][dx^(2)+dy^(2)+dz^(2)].]:}\begin{align*}
d s^{2}= & -\left[1-\frac{2 M}{r}+O\left(\frac{1}{r^{2}}\right)\right] d t^{2}-\left[\frac{4 a M}{r^{3}}+O\left(\frac{1}{r^{4}}\right)\right][x d y-y d x] \\
& +\left[1+O\left(\frac{1}{r}\right)\right]\left[d x^{2}+d y^{2}+d z^{2}\right] . \tag{33.6'}
\end{align*}
Direct comparison with the "standard form" [equation (19.13)] of the metric far from a stationary rotating source reveals that (1) the parameter MM is, indeed, the mass of the black hole; and (2) the intrinsic angular momentum vector of the black hole is
{:(33.7)S=(aM)del//del z=(aM)*((" unit vector pointing along polar axis ")/(" of Boyer-Lindquist coordinates ")).:}\begin{equation*}
\boldsymbol{S}=(a M) \partial / \partial z=(a M) \cdot\binom{\text { unit vector pointing along polar axis }}{\text { of Boyer-Lindquist coordinates }} . \tag{33.7}
\end{equation*}
The charge is defined for the black hole, as for any source, by a Gaussian flux integral of its electric field over a closed surface surrounding the hole. The electric
The rest of this chapter is Track 2. To be prepared for it, one needs to have covered the Track-2 part of Chapter 32 (gravitational collapse). In reading it, one will be helped greatly by Chapter 25 (orbits in Schwarzschild geometry). The rest of this chapter is needed as preparation for Chapter 34 (singularities and global methods).
The metric far outside a black hole: imprints of mass and angular momentum
The electromagnetic field far outside a black hole:
(1) electric field
(2) magnetic field
(3) magnetic dipole moment
Nonspherical shape of hole's geometry
field in the asymptotic rest frame of the black hole has as its orthonormal components
{:[E_( hat(r))=E_(r)=F_(rt)=Q//r^(2)+O(1//r^(3))],[(33.8)E_( hat(theta))=E_(theta)//r=F_(theta t)//r=O(1//r^(4))],[E_( hat(phi))=E_(phi)//r sin theta=F_(theta t)//r sin theta=0]:}\begin{gather*}
E_{\hat{r}}=E_{r}=F_{r t}=Q / r^{2}+O\left(1 / r^{3}\right) \\
E_{\hat{\theta}}=E_{\theta} / r=F_{\theta t} / r=O\left(1 / r^{4}\right) \tag{33.8}\\
\mathrm{E}_{\hat{\phi}}=E_{\phi} / r \sin \theta=F_{\theta t} / r \sin \theta=0
\end{gather*}
Hence, the electric field is purely radial with a Gaussian flux integral of 4pi Q4 \pi Q, which reveals QQ to be the black hole's charge.
A similar calculation of the dominant components of the magnetic field reveals
for the magnetic moment of the black hole.
Just as the rotation of the black hole produces a magnetic field, so it also produces nonspherical deformations in the gravitational field of the black hole [see Hernandez (1967) for quantitative discussion]. But those deformations, like the magnetic moment, are not freely specifiable. They are determined uniquely by the mass, charge, and angular momentum of the black hole.
§33.4. SYMMETRIES AND FRAME DRAGGING
The metric components (33.2) of a Kerr-Newman black hole are independent of the Boyer-Lindquist time coordinate tt and angular coordinate phi\phi. This means (see §25.2) that
are Killing vectors associated with the stationarity (time-translation invariance) and axial symmetry of the black hole. The scalar products of these Killing vectors with themselves and each other are
Since Killing vectors are geometric properties of spacetime, with existence independent of any and all coordinate systems, their scalar products also have coor-dinate-free meaning. It so happens (not by chance, but by careful choice of coordinates!) that the Boyer-Lindquist metric components g_(tt),g_(t phi)g_{t t}, g_{t \phi}, and g_(phi phi)g_{\phi \phi} are equal to these coordinate-independent scalar products. Thus g_(tt),g_(t phi)g_{t t}, g_{t \phi}, and g_(phi phi)g_{\phi \phi} can be thought of as three scalar fields which embody information about the symmetries of spacetime. By contrast, the metric coefficients g_(rr)=rho^(2)//Deltag_{r r}=\rho^{2} / \Delta and g_(theta theta)=rho^(2)g_{\theta \theta}=\rho^{2} carry no information at all about the symmetries.* They depend, for their existence and values, on the specific Boyer-Lindquist choice of coordinates.
Any observer who moves along a world line of constant (r,theta)(r, \theta) with uniform angular velocity sees an unchanging spacetime geometry in his neighborhood. Hence, such an observer can be thought of as "stationary" relative to the local geometry. If and only if his angular velocity is zero, that is, if and only if he moves along a world line of constant (r,theta,phi)(r, \theta, \phi), will he also be "static" relative to the black hole's asymptotic Lorentz frame (i.e., relative to the "distant stars").
The precise definition of "angular velocity relative to the asymptotic rest frame" -or simply "angular velocity"-is
{:(33.13a)Omega-=(d phi)/(dt)=(d phi//d tau)/(dt//d tau)=(u^(phi))/(u^(t)):}\begin{equation*}
\Omega \equiv \frac{d \phi}{d t}=\frac{d \phi / d \tau}{d t / d \tau}=\frac{u^{\phi}}{u^{t}} \tag{33.13a}
\end{equation*}
(see exercise 33.2). In terms of Omega\Omega, the Killing vectors, and the scalar products of Killing vectors, the 4-velocity of a stationary observer is
A stationary observer is static if and only if Omega\Omega vanishes.
The stationary observers at given r,thetar, \theta cannot have any and every angular velocity. Only those values of Omega\Omega are allowed for which the 4 -velocity u\boldsymbol{u} lies inside the future light cone-i.e., for which
and it is assumed that S//M=a > 0S / M=a>0. The following features of these limits are noteworthy. (1) Far from the black hole, one has rOmega_(min)=-1r \Omega_{\min }=-1 and rOmega_(max)=+1r \Omega_{\max }=+1, corresponding to the standard limits imposed by the speed of light in flat spacetime.
(2) With decreasing radius, Omega_(min)\Omega_{\min } increases ("dragging of inertial frames"). Finally, when g_(tt)g_{t t} reaches zero, i.e., at
Omega_("min ")\Omega_{\text {min }} becomes zero. At and inside this surface, all stationary observers must orbit the black hole with positive angular velocity. Thus, static observers exist outside and only outside r=r_(0)(theta)r=r_{0}(\theta). For this reason r=r_(0)(theta)r=r_{0}(\theta) is called the "static limit"; see Box 33.2. (3) As one moves through the static limit into the "ergosphere," one sees the allowed range of angular velocities become ever more positive (ever more "frame dragging"). At the same time, one sees the allowed range narrow down, until finally, at the horizon
the limits Omega_(min)\Omega_{\min } and Omega_(max)\Omega_{\max } coalesce (omega^(2)=g_(tt)//g_(phi phi))\left(\omega^{2}=g_{t t} / g_{\phi \phi}\right). Thus, at the horizon there are no stationary observers. All timelike world lines point inward. There is no escape from the black hole's "pull."
Further features of stationary observers and "frame dragging" are explored in the exercises.
EXERCISES
Exercise 33.1. KERR DESCRIPTION OF KILLING VECTORS
(a) Use the transformation law from Boyer-Lindquist coordinates to Kerr coordinates [equation (4) of Box 33.2] to show that
in accordance with equations (33.19a,b).
(b) Show that for a stationary observer (world line of constant r,thetar, \theta ), the angular velocity expressed in terms of Kerr coordinates is
Omega-=d phi//dt=d widetilde(phi)//d widetilde(grad)=u^( widetilde(phi))//u^( widetilde(V))\Omega \equiv d \phi / d t=d \widetilde{\phi} / d \widetilde{\nabla}=u^{\widetilde{\phi}} / u^{\widetilde{V}}
so that the entire discussion of stationary observers in terms of Kerr coordinates is identical to the discussion in terms of Boyer-Lindquist coordinates. Differences between the coordinate systems show up only when one moves along world lines of changing rr. Reconcile this fact with the fact that both coordinate systems use the same coordinates (r,theta)(r, \theta) but different time and aximuthal coordinates ( t,phit, \phi versus widetilde(V), widetilde(phi)\widetilde{V}, \widetilde{\phi} ).
Exercise 33.2. OBSERVATIONS OF ANGULAR VELOCITY
An observer, far from a black hole and at rest in the hole's asymptotic Lorentz frame, watches (with his eyes) as a particle moves along a stationary (nongeodesic) orbit near the black hole. Let Omega=d phi//dt\Omega=d \phi / d t be the particle's angular velocity, as defined and discussed above. The distant observer uses his stopwatch to measure the time required for the particle to make one complete circuit around the black hole (one complete circuit relative to the distant observer himself; i.e., relative to the hole's asymptotic Lorentz frame).
(a) Show that the circuit time measured is 2pi//Omega2 \pi / \Omega. Thus, Omega\Omega can be regarded as the particle's "angular velocity as measured from infinity."
(b) Let the observer moving with the particle measure its circuit time relative to the asymptotic Lorentz frame, using his eyes and a stopwatch he carries. Show that his answer for the circuit time must be
(a) Place a rigid, circular mirror ("ring mirror") at fixed (r,theta)(r, \theta) around a black hole. Let an observer at (r,theta)(r, \theta) with angular velocity Omega\Omega emit a flash of light. Some of the photons will get caught by the mirror and will skim along its surface, circumnavigating the black hole in the positive- phi\phi direction. Others will get caught and will skim along in the negative- phi\phi direction. Show that the observer will receive back the photons from both directions simultaneously only if his angular velocity is
Thus in this case, and only in this case, can the observer regard the +phi+\phi and -phi-\phi directions as equivalent in terms of local geometry. Put differently, in this case and only in this case is the observer "nonrotating relative to the local spacetime geometry." Thus, it is appropriate to use the name "locally nonrotating observer" for an observer who moves with the angular velocity Omega=omega(r,theta)\Omega=\omega(r, \theta).
(b) Associated with the axial symmetry of a black hole is a conserved quantity, p_(phi)-=p*xi_((phi))p_{\phi} \equiv \boldsymbol{p} \cdot \boldsymbol{\xi}_{(\phi)}, for geodesic motion. This quantity for any particle-whether it is moving along a geodesic or not-is called the "component of angular momentum along the black hole's spin axis," or simply the particle's "angular momentum." (See §33.5\S 33.5§ below.) Show that of all stationary observers at fixed (r,theta)(r, \theta), only the "locally nonrotating observer" has zero angular momentum. [Note: Bardeen, Press, and Teukolsky (1972) have shown that the "locally nonrotating observer" can be a powerful tool in the analysis of physical processes near a black hole.]
Exercise 33.4. ORTHONORMAL FRAMES OF LOCALLY NONROTATING OBSERVERS
(a) Let spacetime be filled with world lines of locally nonrotating observers, and let each such observer carry an orthonormal frame with himself. Show that the spatial orientations of these frames can be so chosen that their basis 1 -forms are
Show that u=-omega^( hat(t))\boldsymbol{u}=-\boldsymbol{\omega}^{\hat{t}} is a rotation-free field of 1 -forms [domega^( hat(t))^^omega^( hat(t))=0:}\left[\boldsymbol{d} \boldsymbol{\omega}^{\hat{t}} \wedge \boldsymbol{\omega}^{\hat{t}}=0\right.; exercise 4.4].
(b) One sometimes meets the mistaken notion that a "locally nonrotating observer" is in some sense locally inertial. To destroy this false impression, verify that: (i) such an observer has nonzero 4 -acceleration,
(ii) if such an observer carries gyroscopes with himself, applying the necessary accelerations at the gyroscope centers of mass, he sees the gyroscopes precess relative to his orthonormal frame (33.21) with angular velocity
[Hints: See exercise 19.2, equation (13.69), and associated discussions. The calculation of the connection coefficients is performed most easily using the methods of differential forms; see §14.6.]
Exercise 33.5. LOCAL LIGHT CONES
Calculate the shapes of the light cones depicted in the Kerr diagram for an uncharged ( Q=0Q=0 ) Kerr black hole (part II.F of Box 33.2). In particular, introduce a new time coordinate
for which the slices of constant tilde(t)\tilde{t} are horizontal surfaces in the Kerr diagram. Then the Kerr diagram plots tilde(t)\tilde{t} vertically, rr radially, and tilde(phi)\tilde{\phi} azimuthally, while holding theta=pi//2\theta=\pi / 2 ("equatorial slice through black hole").
(a) Show that the light cone emanating from given tilde(t),r, tilde(phi)\tilde{t}, r, \tilde{\phi} has the form
(dr)/(d( widetilde(t)))=a((d( tilde(phi)))/(d( tilde(t))))-(2M//r)/(1+2M//r)+-sqrt((1)/((1+2M//r)^(2))-(r^(2)(d( tilde(phi))//d( tilde(t)))^(2))/(1+2M//r))\frac{d r}{d \widetilde{t}}=a\left(\frac{d \tilde{\phi}}{d \tilde{t}}\right)-\frac{2 M / r}{1+2 M / r} \pm \sqrt{\frac{1}{(1+2 M / r)^{2}}-\frac{r^{2}(d \tilde{\phi} / d \tilde{t})^{2}}{1+2 M / r}}
(b) Show that the light cone slices through the surface of constant radius along the curves
{:(33.26b)dr//d tilde(t)=0","d tilde(phi)//d tilde(t)=Omega_(min)" and "Omega_(max):}\begin{equation*}
d r / d \tilde{t}=0, d \tilde{\phi} / d \tilde{t}=\Omega_{\min } \text { and } \Omega_{\max } \tag{33.26b}
\end{equation*}
where Omega_("min ")\Omega_{\text {min }} and Omega_("max ")\Omega_{\text {max }} are given by expressions (33.15a,b) (minimum and maximum allowed angular velocities for stationary observers).
(c) Show that at the static limit, r=r_(0)(pi//2)r=r_{0}(\pi / 2), the light cone is tangent to a curve of constant r,theta, widetilde(phi)r, \theta, \widetilde{\phi}.
(d) Show that the light cone slices the surface of constant widetilde(phi)\widetilde{\phi} along the curves
{:(33.26c)(d( widetilde(phi)))/(d( widetilde(t)))=0","quad(dr)/(d( widetilde(t)))=-1" and "(1-2M//r)/(1+2M//r).:}\begin{equation*}
\frac{d \widetilde{\phi}}{d \widetilde{t}}=0, \quad \frac{d r}{d \widetilde{t}}=-1 \text { and } \frac{1-2 M / r}{1+2 M / r} . \tag{33.26c}
\end{equation*}
(e) Show that the light cone is tangent to the horizon.
(f) Make pictures of the shapes of the light cone as a function of radius.
(g) Describe qualitatively how the light cone must look near the horizon in Boyer-Lindquist coordinates. (Note: it will look "crazy" because the coordinates are singular at the horizon.)
§33.5. EQUATIONS OF MOTION FOR TEST PARTICLES [Carter (1968a)]
Let a test particle with electric charge ee and rest mass mu\mu move in the external fields of a black hole. Were there no charge down the black hole, the test particle would move along a geodesic (zero 4 -acceleration). But the charge produces an electromagnetic field, which in turn produces a Lorentz force on the particle: mu a=eF*u\mu \boldsymbol{a}=e \boldsymbol{F} \cdot \boldsymbol{u}. (Here u\boldsymbol{u} is the particle's 4 -velocity, and a-=grad_(u)u\boldsymbol{a} \equiv \boldsymbol{\nabla}_{\boldsymbol{u}} \boldsymbol{u} is its 4-acceleration.)
The geodesic equation, a=0\boldsymbol{a}=0, for the uncharged case is equivalent to Hamilton's equations
{:(33.27a)dx^(mu)//d lambda=delK//delp_(mu)","quad dp_(mu)//d lambda=-delK//delx^(mu):}\begin{equation*}
d x^{\mu} / d \lambda=\partial \mathscr{K} / \partial p_{\mu}, \quad d p_{\mu} / d \lambda=-\partial \mathscr{K} / \partial x^{\mu} \tag{33.27a}
\end{equation*}
where lambda\lambda is an affine parameter so normalized that
{:(33.27b)d//d lambda=p=4"-momentum, ":}\begin{equation*}
d / d \lambda=\boldsymbol{p}=4 \text {-momentum, } \tag{33.27b}
\end{equation*}
(see exercise 25.2). Similarly (see exercise 33.6) the Lorentz-force equation, mu a=\mu \boldsymbol{a}=eF*ue \boldsymbol{F} \cdot \boldsymbol{u}, for the charged case is equivalent to Hamilton's equations written in terms of position x^(mu)x^{\mu} and "generalized momentum" pi_(mu)\pi_{\mu} :
{:(33.28a)dx^(mu)//d lambda=delH//delpi_(mu)","quad dpi_(mu)//d lambda=-delH//delx^(mu):}\begin{equation*}
d x^{\mu} / d \lambda=\partial \mathscr{H} / \partial \pi_{\mu}, \quad d \pi_{\mu} / d \lambda=-\partial \mathscr{H} / \partial x^{\mu} \tag{33.28a}
\end{equation*}
The form of the superhamiltonian K\mathscr{K}, in terms of the metric coefficients at the particle's location, g^(mu nu)(x^(alpha))g^{\mu \nu}\left(x^{\alpha}\right), and the particle's charge ee and generalized momentum pi_(mu)\pi_{\mu}, is
[See §7.3\S 7.3§ of Goldstein (1959) for the analogous superhamiltonian in flat spacetime.]
Superhamiltonian for a charged test particle in any electromagnetic field in curved spacetime
The first of Hamilton's equations for this superhamiltonian reduces to
{:(33.29a)p^(mu)-=(4"-momentum ")-=dx^(mu)//d lambda=pi^(mu)-eA^(mu):}\begin{equation*}
p^{\mu} \equiv(4 \text {-momentum }) \equiv d x^{\mu} / d \lambda=\pi^{\mu}-e A^{\mu} \tag{33.29a}
\end{equation*}
(value of pi^(mu)\pi^{\mu} in terms of p^(mu),ep^{\mu}, e, and A^(mu)A^{\mu} ); the second, when combined with the first, reduces to the Lorentz-force equation
{:(33.29b){:[dp^(mu)//d lambda+Gamma^(mu)_(alpha beta)p^(alpha)p^(beta)=eF^(mu nu)p_(v).],[[p_(nu)," not "u_(v):}" because "],[lambda=tau//mu]].:}\left.\begin{array}{c}
d p^{\mu} / d \lambda+\Gamma^{\mu}{ }_{\alpha \beta} p^{\alpha} p^{\beta}=e F^{\mu \nu} p_{v} . \tag{33.29b}\\
{\left[p_{\nu}, \text { not } u_{v}\right. \text { because }} \\
\lambda=\tau / \mu
\end{array}\right] .
For a Kerr-Newman black hole, the vector potential in Boyer-Lindquist coordinates can be put in the form
Vector potential for a charged black hole
"Constants of motion" for a charged test particle moving around a charged black hole:
(1) "energy at infinity" EE
(2) "axial component of angular momentum " L_(z)L_{z}
reduces to the Faraday 2 -form of equation (33.5).
There is good reason for going through all this formalism, rather than tackling head-on the Lorentz-force equation in its most elementary coordinate version,
The Hamiltonian formalism enables one to discover immediately two constants of the motion; the elementary Lorentz-force equation does not. The key point is that the components A_(mu)A_{\mu} of A\boldsymbol{A} [equation (33.30)] and the components g^(mu nu)g^{\mu \nu} of the metric [inverse of g_(mu nu)g_{\mu \nu} of equation (33.2); see (33.35)] are independent of tt and phi\phi (stationarity and axial symmetry of both the electromagnetic field and the spacetime geometry). Consequently, the superhamiltonian is also independent of tt and phi\phi; and therefore Hamilton's equation
dpi_(alpha)//d lambda=-delK//delx^(alpha)d \pi_{\alpha} / d \lambda=-\partial \mathscr{K} / \partial x^{\alpha}
guarantees that pi_(t)\pi_{t} and pi_(phi)\pi_{\phi} are constants of the motion.
Far from the black hole, where the vector potential vanishes and the metric becomes
ds^(2)=-dt^(2)+dr^(2)+r^(2)(dtheta^(2)+sin^(2)theta dphi^(2))d s^{2}=-d t^{2}+d r^{2}+r^{2}\left(d \theta^{2}+\sin ^{2} \theta d \phi^{2}\right)
the constants of the motion become
{:[pi_(t)=p_(t)=-p^(t)=-" energy, "],[pi_(phi)=p_(phi)=rp^( hat(phi))=((" projection of angular momentum ")/(" along black hole's rotation axis ")).]:}\begin{gathered}
\pi_{t}=p_{t}=-p^{t}=- \text { energy, } \\
\pi_{\phi}=p_{\phi}=r p^{\hat{\phi}}=\binom{\text { projection of angular momentum }}{\text { along black hole's rotation axis }} .
\end{gathered}
Thus it is appropriate to adopt the names and notation
{:[(33.31a)E-=(" "energy at infinity" ")-=-pi_(t)=-(p_(t)+eA_(t))],[(33.31b)L_(z)-=([" "axial component of angular "],[" momentum", or simply "],[" "angular momentum" "])-=pi_(phi)=p_(phi)+eA_(phi)]:}\begin{align*}
E & \equiv(\text { "energy at infinity" }) \equiv-\pi_{t}=-\left(p_{t}+e A_{t}\right) \tag{33.31a}\\
L_{z} & \equiv\left(\begin{array}{l}
\text { "axial component of angular } \\
\text { momentum", or simply } \\
\text { "angular momentum" }
\end{array}\right) \equiv \pi_{\phi}=p_{\phi}+e A_{\phi} \tag{33.31b}
\end{align*}
for the constants of the motion -pi_(t)-\pi_{t} and pi_(phi)\pi_{\phi}.
A third constant of the motion is the particle's rest mass
In general, four constants of the motion are needed to determine uniquely the orbit of a particle through four-dimensional spacetime. If the black hole were to possess an additional symmetry-e.g., if it were spherical, rather than merely axially symmetric-then automatically there would be a fourth constant of the motion. But in general, black holes are not spherical; so test-particle motion around a black hole possesses only three obvious constants. It is rather remarkable, then, that a constant turns out to exist. It was discovered by Carter (1968a), using Hamilton-Jacobi methods. As of 1973, nobody has given a cogent geometric explanation of why this fourth constant should exist-although hints of an explanation may be found in Carter (1968c) and Walker and Penrose (1970).
Carter's "fourth constant" of the motion, as derived in exercise 33.7, is
obtained by combining Q,L_(z)\mathscr{Q}, L_{z}, and EE, is often used in place of Q\mathscr{Q}. Whereas Q\mathscr{Q} can be negative, K\mathscr{K} is always nonnegative:
{:[K=p_(theta)^(2)+(L_(z)-aEsin^(2)theta)^(2)//sin^(2)theta+a^(2)mu^(2)cos^(2)theta],[ >= 0" everywhere "],[=0" only for case of photon "(mu=0)" moving along polar axis "(theta=0","pi)]:}\begin{aligned}
\mathscr{K} & =p_{\theta}^{2}+\left(L_{z}-a E \sin ^{2} \theta\right)^{2} / \sin ^{2} \theta+a^{2} \mu^{2} \cos ^{2} \theta \\
& \geq 0 \text { everywhere } \\
& =0 \text { only for case of photon }(\mu=0) \text { moving along polar axis }(\theta=0, \pi)
\end{aligned}
The contravariant components of the test particle's 4-momentum, p^(alpha)=dx^(alpha)//d lambdap^{\alpha}=d x^{\alpha} / d \lambda, are readily expressed in terms of the constants E,L_(z),mu,QE, L_{z}, \mu, \mathscr{Q}, by combining equations (33.31) with the metric coefficients (33.2) and the components of the vector potential (33.30). The result is
{:[(33.32a)rho^(2)d theta//d lambda=sqrtTheta","],[(33.32b)rho^(2)dr//d lambda=sqrtR","],[(33.32c)rho^(2)d phi//d lambda=-(aE-L_(z)//sin^(2)theta)+(a//Delta)P],[(33.32d)rho^(2)dt//d lambda=-a(aEsin^(2)theta-L_(z))+(r^(2)+a^(2))Delta^(-1)P]:}\begin{gather*}
\rho^{2} d \theta / d \lambda=\sqrt{\boldsymbol{\Theta}}, \tag{33.32a}\\
\rho^{2} d r / d \lambda=\sqrt{R}, \tag{33.32b}\\
\rho^{2} d \phi / d \lambda=-\left(a E-L_{z} / \sin ^{2} \theta\right)+(a / \Delta) P \tag{33.32c}\\
\rho^{2} d t / d \lambda=-a\left(a E \sin ^{2} \theta-L_{z}\right)+\left(r^{2}+a^{2}\right) \Delta^{-1} P \tag{33.32d}
\end{gather*}
Equations of motion for charged test particles
Here rho^(2)=r^(2)+a^(2)cos^(2)theta\rho^{2}=r^{2}+a^{2} \cos ^{2} \theta as defined in equation (33.3b), and the functions Theta,R,P\Theta, R, P are defined by
When working in Kerr coordinates (to avoid the coordinate singularity at the horizon), one must replace equations (33.32c) and (33.32d) by
{:[(33.32c')rho^(2)d widetilde(V)//d lambda=-a(aEsin^(2)theta-L_(z))+(r^(2)+a^(2))Delta^(-1)(sqrtR+P)","],[(33.32d')rho^(2)d widetilde(phi)//d lambda=-(aE-L_(z)//sin^(2)theta)+aDelta^(-1)(sqrtR+P)]:}\begin{gather*}
\rho^{2} d \widetilde{V} / d \lambda=-a\left(a E \sin ^{2} \theta-L_{z}\right)+\left(r^{2}+a^{2}\right) \Delta^{-1}(\sqrt{R}+P), \tag{33.32c'}\\
\rho^{2} d \widetilde{\phi} / d \lambda=-\left(a E-L_{z} / \sin ^{2} \theta\right)+a \Delta^{-1}(\sqrt{R}+P) \tag{33.32d'}
\end{gather*}
[These follow from (33.32) and the transformation between the two coordinate systems-see equations (4) of Box 33.2.] In the above equations, the signs of sqrtR\sqrt{R} and sqrtTheta\sqrt{\Theta} can be chosen independently; but once chosen, they must be used consistently everywhere.
Applications of these equations of motion will play a key role in the rest of this chapter.
EXERCISES
Exercise 33.6. SUPERHAMILTONIAN FOR CHARGED-PARTICLE MOTION
Show that Hamilton's equations (33.28a) for the Hamiltonian (33.28b) reduce to equation (33.29a) for the value of the generalized momentum, and to the Lorentz force equation (33.29b). [Hint: Use the relation {:(g^(alpha beta)g_(beta gamma))_(,mu)=0.]\left.\left(g^{\alpha \beta} g_{\beta \gamma}\right)_{, \mu}=0.\right]
Exercise 33.7. HAMILTON-JACOBI DERIVATION OF EQUATIONS OF MOTION [Based on Carter (1968a)]
Derive the first-order equations of motion (33.32) for a charged particle moving in the external fields of a Kerr-Newman black hole. Use the Hamilton-Jacobi method [Boxes 25.3 and 25.4 of this book; also Chapter 9 of Goldstein (1959)], as follows.
(a) Throughout the superhamiltonian K\mathscr{K} of equation (33.28b), replace the generalized momentum pi_(alpha)\pi_{\alpha} by the gradient del S//delx^(alpha)\partial S / \partial x^{\alpha} of the Hamilton-Jacobi function.
(b) Write down the Hamilton-Jacobi equation [generalization of equation (2) of Box 25.4] in the form
(d) Use these metric components and the components (33.30) of the vector potential to bring the Hamilton-Jacobi equation (33.33) into the concrete form
(e) Solve this Hamilton-Jacobi equation by separation of variables. [Hint: Because the equation has no explicit dependence on lambda,phi\lambda, \phi, or tt, the solution must take the form
where the values of the "integration constants" follow from del S//del lambda=-H,del S//del t=pi_(t)\partial S / \partial \lambda=-\mathscr{H}, \partial S / \partial t=\pi_{t}, del S//del phi=pi_(phi)\partial S / \partial \phi=\pi_{\phi}. Insert this assumed form into (33.35) and solve for S_(r)(r)S_{r}(r) and S_(theta)(theta)S_{\theta}(\theta) to obtain
{:(33.36b)S_(r)=intDelta^(-1)sqrtRdr","quadS_(theta)=intsqrtThetad theta:}\begin{equation*}
S_{r}=\int \Delta^{-1} \sqrt{R} d r, \quad S_{\theta}=\int \sqrt{\Theta} d \theta \tag{33.36b}
\end{equation*}
where R(r)R(r) and Theta(theta)\Theta(\theta) are the functions defined in equation (33.33). Notice that the constant Q\mathscr{Q} arises naturally as a "separation-of-variables constant" in this procedure. It was in this way that Carter originally discovered Q\mathcal{Q}, following Misner's suggestion that he seek analogies to a constant in Newtonian dipole fields (Corben and Stehle, 1960, p. 209).]
(f) By successively setting del S//del[Q+(L_(z)-aE)^(2)],del S//delmu^(2),del S//del E\partial S / \partial\left[\mathscr{Q}+\left(L_{z}-a E\right)^{2}\right], \partial S / \partial \mu^{2}, \partial S / \partial E, and del S//delL_(z)\partial S / \partial L_{z} to zero, obtain the following equations describing the test-particle orbits:
{:[(33.37a)int^(theta)(d theta)/(sqrtTheta)=int^(r)(dr)/(sqrtR)],[(33.37b)lambda=int^(theta)(a^(2)cos^(2)theta)/(sqrtTheta)d theta+int^(r)(r^(2))/(sqrtR)dr],[(33.37c)t=int^(theta)(-a(aEsin^(2)theta-L_(z)))/(sqrtTheta)d theta+int^(r)((r^(2)+a^(2))P)/(DeltasqrtR)dr],[(33.37~d)phi=int(-(aEsin^(2)theta-L_(z)))/(sin^(2)thetasqrtTheta)d theta+int(aP)/(DeltasqrtR)dr]:}\begin{align*}
\int^{\theta} \frac{d \theta}{\sqrt{\Theta}} & =\int^{r} \frac{d r}{\sqrt{R}} \tag{33.37a}\\
\lambda & =\int^{\theta} \frac{a^{2} \cos ^{2} \theta}{\sqrt{\Theta}} d \theta+\int^{r} \frac{r^{2}}{\sqrt{R}} d r \tag{33.37b}\\
t & =\int^{\theta} \frac{-a\left(a E \sin ^{2} \theta-L_{z}\right)}{\sqrt{\Theta}} d \theta+\int^{r} \frac{\left(r^{2}+a^{2}\right) P}{\Delta \sqrt{R}} d r \tag{33.37c}\\
\phi & =\int \frac{-\left(a E \sin ^{2} \theta-L_{z}\right)}{\sin ^{2} \theta \sqrt{\Theta}} d \theta+\int \frac{a P}{\Delta \sqrt{R}} d r \tag{33.37~d}
\end{align*}
(g) By differentiating these equations and combining them, obtain the equations of motion (33.32) cited in the text.
(h) Derive equations (33.31) for E,L_(z),muE, L_{z}, \mu, and Q\mathscr{Q} by setting del S//delx^(alpha)=pi_(alpha)=p_(alpha)+eA_(alpha)\partial S / \partial x^{\alpha}=\pi_{\alpha}=p_{\alpha}+e A_{\alpha}.
§33.6. PRINCIPAL NULL CONGRUENCES
Two special families of photon trajectories "mold themselves into" the Kerr-Newman geometry in an especially harmonious way. They are called the "principal null congruences" of the geometry. ("Congruence" is an elegant word that means "space-
Principal null congruences for the spacetime geometry of a black hole
filling family of curves.") These congruences are the solutions to the test-particle equations of motion (33.32) with
{:[(33.38a)mu=0" (zero rest mass; photon) "","],[(33.38b)e=0(" zero charge on photon ")","],[(33.38c)L_(z)=aEsin^(2)thetaquad([" a permissible value for "L_(z)],[" only because "d theta//d lambda" turns "],[" out to be zero "])","],[(33.38d)Q=-(L_(z)-aE)^(2)=-a^(2)E^(2)cos^(4)theta.]:}\begin{gather*}
\mu=0 \text { (zero rest mass; photon) }, \tag{33.38a}\\
e=0(\text { zero charge on photon }), \tag{33.38b}\\
L_{z}=a E \sin ^{2} \theta \quad\left(\begin{array}{l}
\text { a permissible value for } L_{z} \\
\text { only because } d \theta / d \lambda \text { turns } \\
\text { out to be zero }
\end{array}\right), \tag{33.38c}\\
\mathscr{Q}=-\left(L_{z}-a E\right)^{2}=-a^{2} E^{2} \cos ^{4} \theta . \tag{33.38d}
\end{gather*}
For these values of the constants of motion, the equations of motion (33.32) reduce to
{:[(33.39a)k^(theta)-=d theta//d lambda=0],[(33.39b)k^(r)-=dr//d lambda=+-E quad{:[" (" "+" " for outgoing photons, "],[" "-" for ingoing) "]:}],[k^(phi)-=d phi//d lambda=aE//Delta],[(33.39c)k^(t)-=dt//d lambda=(r^(2)+a^(2))E//Delta]:}\begin{gather*}
k^{\theta} \equiv d \theta / d \lambda=0 \tag{33.39a}\\
k^{r} \equiv d r / d \lambda= \pm E \quad \begin{array}{c}
\text { (" }+ \text { " for outgoing photons, } \\
\text { "-" for ingoing) }
\end{array} \tag{33.39b}\\
k^{\phi} \equiv d \phi / d \lambda=a E / \Delta \\
k^{t} \equiv d t / d \lambda=\left(r^{2}+a^{2}\right) E / \Delta \tag{33.39c}
\end{gather*}
Significance of the principal null congruences
In what sense are these photon trajectories more interesting than others? (1) They mold themselves to the spacetime curvature in such a way that, if C_(alpha beta gamma delta)C_{\alpha \beta \gamma \delta} is the Weyl conformal tensor (§13.5), and ^(**)C_(alpha beta gamma delta)=epsilon_(alpha beta mu nu)C^(|mu nu|)_(gamma delta){ }^{*} C_{\alpha \beta \gamma \delta}=\epsilon_{\alpha \beta \mu \nu} C^{|\mu \nu|}{ }_{\gamma \delta} is its dual, then
[This relationship implies that the Kerr-Newman geometry is of "Petrov-Pirani type D" and that these photon trajectories are "doubly degenerate, principal null congruences." For details of the meanings and implications of these terms see, e.g., §8 of Sachs (1964), or Ehlers and Kundt (1962), or the original papers by Petrov (1954, 1969 ) and Pirani (1957).] (2) By suitable changes of coordinates (exercise 33.8), one can bring the Kerr-Newman metric into the form
{:(33.41)ds^(2)=(eta_(alpha beta)+2Hk_(alpha)k_(beta))dx^(alpha)dx^(beta)",":}\begin{equation*}
d s^{2}=\left(\eta_{\alpha \beta}+2 H k_{\alpha} k_{\beta}\right) d x^{\alpha} d x^{\beta}, \tag{33.41}
\end{equation*}
where HH is a scalar field and k_(alpha)k_{\alpha} are the components of the wave vector for one of the principal null congruences (either one; but not both!). [This was the property of the Kerr-Newman metric that led to its original discovery (Kerr, 1963). For further detail on metrics of this form, see Kerr and Schild (1965).] (3) In Kerr coordinates (Box 33.2), the ingoing principal null congruence is
{:[(33.42a)r=-E lambda","quad theta=" const, "quad widetilde(phi)=" const, "quad widetilde(V)=" const. "],[{[" arbitrary normalization "],[" factor; can be removed "],[" by redefinition of "lambda]]]:}\begin{align*}
r=- & E \lambda, \quad \theta=\text { const, } \quad \widetilde{\phi}=\text { const, } \quad \widetilde{V}=\text { const. } \tag{33.42a}\\
& \left\{\begin{array}{l}
\text { arbitrary normalization } \\
\text { factor; can be removed } \\
\text { by redefinition of } \lambda
\end{array}\right]
\end{align*}
These ingoing photon world lines are the generators of the conical surface widetilde(V)=\widetilde{V}= const. in the Kerr diagram of Box 33.2. (4) The only kind of particle that can remain forever at the horizon is a photon with world line in the outgoing principal null congruence (exercise 33.9). Such photon world lines are "generators" of the horizon (dotted curves with a "barber-pole twist" in Kerr diagram of Box 33.2). They have angular velocity
(a) Show that in Kerr coordinates the ingoing null congruence (33.39) has the form (33.42a). Also show that the covariant components of the wave vector-after changing to a new affine parameter lambda_("new ")=lambda_("old ")E\lambda_{\text {new }}=\lambda_{\text {old }} E-are
(b) Introduce new coordinates tilde(t),x,y,z\tilde{t}, x, y, z, defined by
and show that in this "Kerr-Schild coordinate system" the metric takes the form
{:(33.44b)ds^(2)=(eta_(alpha beta)+2Hk_(alpha)^((in))k_(beta)^((in)))dx^(alpha)dx^(beta)",":}\begin{equation*}
d s^{2}=\left(\eta_{\alpha \beta}+2 H k_{\alpha}^{(\mathrm{in})} k_{\beta}^{(\mathrm{in})}\right) d x^{\alpha} d x^{\beta}, \tag{33.44b}
\end{equation*}
where
{:[(33.44c)H=(Mr-(1)/(2)Q^(2))/(r^(2)+a^(2)(z//r)^(2))","],[(33.44d)k_(alpha)^(("in "))dx^(alpha)=-(r(xdx+ydy)-a(xdy-ydx))/(r^(2)+a^(2))-(zdz)/(r)-d bar(t)]:}\begin{gather*}
H=\frac{M r-\frac{1}{2} Q^{2}}{r^{2}+a^{2}(z / r)^{2}}, \tag{33.44c}\\
k_{\alpha}^{(\text {in })} d x^{\alpha}=-\frac{r(x d x+y d y)-a(x d y-y d x)}{r^{2}+a^{2}}-\frac{z d z}{r}-d \bar{t} \tag{33.44d}
\end{gather*}
For the transformation to analogous coordinates in which
ds^(2)=(eta_(alpha beta)+2Hk_(alpha)^((out )k_(beta)^((out) ))dx^(alpha)dx^(beta).d s^{2}=\left(\eta_{\alpha \beta}+2 H k_{\alpha}^{\text {(out }} k_{\beta}^{\text {(out) }}\right) d x^{\alpha} d x^{\beta} .
see, e.g., Boyer and Lindquist (1967).
Exercise 33.9. NULL GENERATORS OF HORIZON
(a) Show that in Kerr coordinates the outgoing principle null congruence is described by the tangent vector
(b) These components of the wave vector become singular at the horizon (Delta=0)(\boldsymbol{\Delta}=0), not because of a singularity in the coordinate system-the coordinates are well-behaved!-but because of poor normalization of the affine parameter. For each outgoing geodesic, let Delta_(0)\Delta_{0}
be a constant, defined as the value of Delta\Delta at the event where the geodesic slices the hypersurface widetilde(V)=0\widetilde{V}=0. Then renormalize the affine parameter for each geodesic
(c) Show that these are the only test-particle trajectories that remain forever on the horizon. [Hint: Examine the light cone.]
§33.7. STORAGE AND REMOVAL OF ENERGY FROM BLACK HOLES [Penrose (1969)]
When an object falls into a black hole, it changes the hole's mass, charge, and intrinsic angular momentum (first law of black-hole dynamics; Box 33.4). If the infalling object is large, its fall produces much gravitational and electromagnetic radiation. To calculate the radiation emitted, and the energy and angular momentum it carries away-which are prerequisites to any calculation of the final state of the black hole-is an enormously difficult task. But if the object is very small (size of object << size of horizon; mass of object ≪\ll mass of hole), and has sufficiently small charge, the radiation it emits in each circuit around the hole is negligible. For example, for gravitational radiation
{:(33.48)((" energy emitted per circuit "))/((" rest mass of object "))∼((" rest mass of object "))/((" mass of hole ")):}\begin{equation*}
\frac{(\text { energy emitted per circuit })}{(\text { rest mass of object })} \sim \frac{(\text { rest mass of object })}{(\text { mass of hole })} \tag{33.48}
\end{equation*}
[see §36.5\S 36.5§; also Bardeen, Press, and Teukolsky (1972)]. Because the energy emitted is negligible, radiation reaction is also negligible, and the object moves very nearly along a test-particle trajectory. In this case, application of the first law of black-hole dynamics is simple and straightforward.
Consider, initially, a small object that falls directly into the black hole from far away. According to the first law, it produces the following changes in the mass, charge, and angular momentum of the black hole:
{:[(33.49a)Delta M=E=(" "energy at infinity" of infalling object ")],[(33.49b)Delta Q=e=" (charge of infalling object) "],[(33.49c)Delta S-=Delta|S|=L_(z)=((" component of object's angular momentum ")/(" on black hole's rotation axis "))]:}\begin{gather*}
\Delta M=E=(\text { "energy at infinity" of infalling object }) \tag{33.49a}\\
\Delta Q=e=\text { (charge of infalling object) } \tag{33.49b}\\
\Delta S \equiv \Delta|\boldsymbol{S}|=L_{z}=\binom{\text { component of object's angular momentum }}{\text { on black hole's rotation axis }} \tag{33.49c}
\end{gather*}
The infalling object will also change the direction of S\boldsymbol{S}. In the black hole's original asymptotic Lorentz frame, its initial angular momentum vector points in the zz-direction,
Consequently, only the zz-component of angular momentum of the infalling object can produce any significant change in the magnitude of S\boldsymbol{S}. But the xx - and yy-components, L_(x)L_{x} and L_(y)L_{y}, can change the direction of S\boldsymbol{S}. If the object has negligible speed at infinity, then it produces the changes (exercise 33.10):
Here a subscript " oo\infty " means the value of a quantity at a point on the orbit far from the black hole (at "infinity").
Consider, next, a more complicated process, first conceived of by Penrose (1969): (1) Shoot a small object AA into the black hole from outside with energy-at-infinity E_(A)E_{A}, charge e_(A)e_{A}, and axial component of angular momentum L_(zA)L_{z A}. (2) When the object is deep down near the horizon, let it explode into two parts, BB and CC, each of which subsequently moves along a new test-particle trajectory, with new constants of the motion e_(B)e_{B} and e_(C),E_(B)e_{C}, E_{B} and E_(C),L_(zB)E_{C}, L_{z B} and L_(zC^(**))L_{z C^{*}} (3) So design the explosion that object BB falls down the hole and gets captured, but object CC escapes back to radial infinity. What will be the changes in mass, charge, and angular momentum of the black hole? According to the first law of black-hole dynamics,
{:[Delta M=([" total energy that distant observers see "],[" fall inward past themselves minus "],[" total energy that they see reemerge "])],[=E_(A)-E_(C).]:}\begin{aligned}
\Delta M & =\left(\begin{array}{l}
\text { total energy that distant observers see } \\
\text { fall inward past themselves minus } \\
\text { total energy that they see reemerge }
\end{array}\right) \\
& =E_{A}-E_{C} .
\end{aligned}
Similarly, Delta Q=e_(A)-e_(C)\Delta Q=e_{A}-e_{C} and Delta S=L_(zA)-L_(zC)\Delta S=L_{z A}-L_{z C}. Not unexpectedly, these changes can be written more simply in terms of the constants of motion for object BB, which went down the hole. View the explosion " A longrightarrow B+CA \longrightarrow B+C " in a local Lorentz frame down near the hole, which is centered on the explosive event. As viewed in that frame, the explosion must satisfy the special relativistic laws of physics (equivalence principle!). In particular, it must obey charge conservation
(p_(A))_("immediately before explosion ")=(p_(B)+p_(C))_("immediately after explosion ")\left(\boldsymbol{p}_{A}\right)_{\text {immediately before explosion }}=\left(\boldsymbol{p}_{B}+\boldsymbol{p}_{C}\right)_{\text {immediately after explosion }}
Moreover, conservation of 4-momentum p\boldsymbol{p} and charge ee implies also conservation of generalized momentum pi-=p-eA\boldsymbol{\pi} \equiv \boldsymbol{p}-e \boldsymbol{A},
and hence also conservation of the components of generalized momentum along the vectors del//del t\partial / \partial t and del//del phi\partial / \partial \phi,
This result restated in words: the changes in mass, charge, and angular momentum are equal to the "energy-at-infinity," charge, and "axial component of angular momentum" that object BB carries inward across the horizon, even though BB may have ended up on a test-particle orbit that does not extend back to radial infinity!
Straightforward extensions of the above thought experiment produce this generalization: In any complicated black-hole process that involves infalling, colliding, and exploding pieces of matter that emit negligible gravitational radiation, the total changes in mass, charge, and angular momentum of the black hole are
{:[(33.52a)Delta M=([" sum of values of energy-at-infinity, "E","],[" for all objects which cross the horizon-with "],[E" evaluated for each object at event of crossing "])","],[(33.52b)Delta Q=((" similar sum, of charges, "e," for ")/(" all objects crossing horizon "))","],[(33.52c)Delta S=((" similar sum of axial components of angular ")/(" momentum, "L_(z)," for all objects crossing horizon ")).]:}\begin{align*}
\Delta M & =\left(\begin{array}{l}
\text { sum of values of energy-at-infinity, } E, \\
\text { for all objects which cross the horizon-with } \\
E \text { evaluated for each object at event of crossing }
\end{array}\right), \tag{33.52a}\\
\Delta Q & =\binom{\text { similar sum, of charges, } e, \text { for }}{\text { all objects crossing horizon }}, \tag{33.52b}\\
\Delta S & =\binom{\text { similar sum of axial components of angular }}{\text { momentum, } L_{z}, \text { for all objects crossing horizon }} . \tag{33.52c}
\end{align*}
This result is not at all surprising. It is precisely what one might expect from the most naive of viewpoints. Not so expected, however, is the following consequence [Penrose (1969)]: By injecting matter into a black hole in a carefully chosen way, one can decrease the total mass-energy of the black hole-i.e., one can extract energy from the hole.
For uncharged infalling objects, the key to energy extraction is the ergosphere [hence its name, coined by Ruffini and Wheeler (1971a) from the Greek word " epsilon rho gamma o nu\epsilon \rho \gamma o \nu " for "work"]. Outside the ergosphere, the Killing vector xi_((t))-=del//del t\xi_{(t)} \equiv \partial / \partial t is timelike, as is the 4-momentum p\boldsymbol{p} of every test particle; and therefore E=-p*xi_((t))E=-\boldsymbol{p} \cdot \xi_{(t)} is necessarily positive. But inside the ergosphere (between the horizon and the static limit), xi_((t))\boldsymbol{\xi}_{(t)}
Changes in M,Q,SM, Q, S for any process
nonradiative black-hole
is spacelike, so for certain choices of timelike momentum vector (certain orbits of uncharged test particles), E=-p*xi_((t))E=-\boldsymbol{p} \cdot \xi_{(t)} is negative, whereas for others it is positive. The orbits of negative EE are confined entirely to the ergosphere. Thus, to inject an uncharged object with negative EE into the black hole-and thereby to extract energy from the hole-one must always change its EE from positive to negative and therefore also change its orbit, after it penetrates into the ergosphere. Of course, this is not difficult in principle-and perhaps not even in practice; see Figure 33.2.
For a charged object, electromagnetic forces alter the region where there exist orbits of negative energy-at-infinity. If the charges of object and hole have opposite sign, then the hole's electromagnetic field pulls inward on the object, giving it more kinetic energy when near the hole than one would otherwise expect. Thus, -p*xi_((t))-\boldsymbol{p} \cdot \boldsymbol{\xi}_{(t)} becomes an overestimate of EE,
{:(33.53)E=-(p-eA)*xi_((t))=-p*xi_((t))+ubrace(eQr//rho^(2)ubrace)_(t[ < 0" if "eQ < 0]):}\begin{equation*}
E=-(\boldsymbol{p}-e \boldsymbol{A}) \cdot \xi_{(t)}=-\boldsymbol{p} \cdot \xi_{(t)}+\underbrace{e Q r / \rho^{2}}_{t[<0 \text { if } e Q<0]} \tag{33.53}
\end{equation*}
and orbits with E < 0E<0 exist in a region somewhat larger than the ergosphere. If, on the other hand, ee and QQ have the same sign, then orbits with E < 0E<0 are confined to a region smaller than the ergosphere. For given values e,Qe, Q, and rest mass mu\mu, the region where there exist orbits with E < 0E<0 is called the "effective ergosphere."
Exercise 33.10. ANGULAR MOMENTUM VECTOR FOR INFALLING PARTICLE
EXERCISE
Derive equations ( 33.49d,e,f33.49 \mathrm{~d}, \mathrm{e}, \mathrm{f} ) for the components L_(x)L_{x} and L_(y)L_{y} of the orbital angular momentum of a particle falling into a black hole. Assume negligible initial speed, E^(2)-mu^(2)~~0E^{2}-\mu^{2} \approx 0.
§33.8. REVERSIBLE AND IRREVERSIBLE TRANSFORMATIONS [Christodoulou (1970), Christodoulou and Ruffini (1971)]
Take a black hole of given mass MM, charge QQ, and angular momentum SS. By injection of small objects, make a variety of changes in M,QM, Q, and SS. Can one pick an arbitrary desired change, Delta M,Delta Q\Delta M, \Delta Q, and Delta S\Delta S, and devise a process that achieves it? Or are there limitations?
The second law of black-hole dynamics (nondecreasing surface area of black hole; Box 33.4; proof in $34.5\$ 34.5 of next chapter) provides a strict limitation.
Then can all values within that limitation be achieved-and can that limitation be discovered by a direct examination of test-particle orbits?
The answer is yes; and, in fact, the limitation was discovered by Christodoulou (1970) and Christodoulou and Ruffini (1971) from an examination of test-particle orbits, independently of and simultaneously with Hawking's (1971) discovery of the second law of black-hole dynamics.
Figure 33.2 .
An advanced civilization has constructed a rigid framework around a black hole, and has built a huge city on that framework. Each day trucks carry one million tons of garbage out of the city to the garbage dump. At the dump the garbage is shoveled into shuttle vehicles which are then, one after another, dropped toward the center of the black hole. Dragging of inertial frames whips each shuttle vehicle into a circling, inward-spiraling orbit near the horizon. When it reaches a certain "ejection point," the vehicle ejects its load of garbage into an orbit of negative energy-at-infinity, E_("garbage ") < 0E_{\text {garbage }}<0. As the garbage flies down the hole, changing the hole's total mass-energy by Delta M=E_("garbage ejected ") < 0\Delta M=E_{\text {garbage ejected }}<0, the shuttle vehicle recoils from the ejection and goes flying back out with more energy-at-infinity than it took down
{:[E_("vehicle out ")=E_("vehicle "+" garbage down ")-E_("garbage ejected ")],[ > E_("vehicle "+" garbage down ")]:}\begin{aligned}
E_{\text {vehicle out }}= & E_{\text {vehicle }+ \text { garbage down }}-E_{\text {garbage ejected }} \\
& >E_{\text {vehicle }+ \text { garbage down }}
\end{aligned}
The vehicle deposits its huge kinetic energy in a giant flywheel adjacent to the garbage dump; and the flywheel turns a generator, producing electricity for the city, while the shuttle vehicle goes back for another load of garbage. The total electrical energy generated with each round trip of the shuttle vehicle is
{:[" (Energy per trip) "=E_("vehicle out ")-(" rest mass of vehicle ")],[=(E_("vehicle "+" garbage down ")-E_("garbage ejected "))-(" rest mass of vehicle ")],[=(" rest mass of vehicle "+" rest mass of garbage "-Delta M)-(" rest mass of vehicle ")],[=" (rest mass of garbage ")+(" amount "","-Delta M","" by which hole's mass decreases ")]:}\begin{aligned}
\text { (Energy per trip) } & =E_{\text {vehicle out }}-(\text { rest mass of vehicle }) \\
& =\left(E_{\text {vehicle }+ \text { garbage down }}-E_{\text {garbage ejected }}\right)-(\text { rest mass of vehicle }) \\
& =(\text { rest mass of vehicle }+ \text { rest mass of garbage }-\Delta M)-(\text { rest mass of vehicle }) \\
& =\text { (rest mass of garbage })+(\text { amount },-\Delta M, \text { by which hole's mass decreases })
\end{aligned}
Thus, not only can the inhabitants of the city use the black hole to convert the entire rest mass of their garbage into kinetic energy of the vehicle, and thence into electrical power, but they can also convert some of the mass of the black hole into electrical power!
To derive the limitation of nondecreasing surface area from properties of testparticle orbits, one must examine what values of energy-at-infinity, EE, are allowed at a given location (r,theta)(r, \theta) outside a black hole. Equations (33.32a,b), when combined, yield the value of EE in terms of a test particle's location (r,theta)(r, \theta), rest mass mu\mu, charge ee, axial component of angular momentum L_(z)L_{z}, and momenta p^(r)=dr//d lambda,p^(theta)=d theta//d lambdap^{r}=d r / d \lambda, p^{\theta}=d \theta / d \lambda in the rr and theta\theta directions:
(One must take the positive square root, +sqrt(beta^(2)-alpha gamma)+\sqrt{\beta^{2}-\alpha \gamma}, rather than the negative square root; positive square root corresponds to 4 -momentum pointing toward future; while negative square root corresponds to past-pointing 4 -momentum; see Figure 33.3.)
Several features of the energy equation (33.54) are noteworthy. (1) For orbits in the equatorial "plane," theta=pi//2\theta=\pi / 2 and p^(theta)-=0p^{\theta} \equiv 0, the energy equation yields an effective potential for radial motion (Box 33.5). (2) Orbits of negative EE must have beta < 0\beta<0 and gamma > 0\gamma>0-which can be achieved only if L_(z)a < 0L_{z} a<0 and/or eQ < 0e Q<0. Thus, one cannot decrease the mass of a black hole without simultaneously decreasing the magnitude of its charge or angular momentum or both. (3) For an orbit at given (r,theta)(r, \theta), with given ee and L_(z),EL_{z}, E is a minimum if p^(r)=p^(theta)=mu=0p^{r}=p^{\theta}=\mu=0. Put differently, the rest mass and the rr - and theta\theta-components of momentum always contribute positively to EE.
By injecting an object into a black hole, produce small changes
in its mass, charge, and angular momentum. For given changes in QQ and SS, what range of changes in MM is possible? Clearly delta M\delta M can be made as large as one wishes by making the rest mass mu\mu sufficiently large. But there will be a lower limit on delta M\delta M. That limit corresponds to the minimum value of EE for given ee and L_(z)L_{z}. The orbit of minimum EE crosses the horizon (otherwise no changes in M,Q,SM, Q, S would occur!), so one can evaluate EE there. At the horizon, as anywhere, a minimum for EE is achieved if mu=p^(r)=p^(theta)=0\mu=p^{r}=p^{\theta}=0. Inserting these values and r=r_(+)r=r_{+}(so Delta=0\Delta=0 ) into equations (33.54), one finds
Properties of test-particle orbits:
(1) EE as function of mu,e,L_(z^('))\mu, e, L_{z^{\prime}}, r,theta,p^(r)r, \theta, p^{r}
(2) effective potential
(3) negative EE requires L_(z)a < 0L_{z} a<0 and/or eQ < 0e Q<0
Changes in black-hole properties due to injection of particles:
Figure 33.3.
Energy-at-infinity EE allowed for a particle of angular momentum L_(z)L_{z} and rest mass mu\mu, which is (1) in the "equatorial plane" theta=pi//2\theta=\pi / 2, (2) at radius r=3M//2r=3 M / 2, (3) of an uncharged ( Q=0Q=0 ) extreme-Kerr ( S=M^(2)S=M^{2} ) black hole. EE is here plotted against L_(z)L_{z}. "Seas" of "positive and negative root" states are shown. The positive root states have energies-at-infinity given by equations (33.54)
and have 4 -momentum vectors pointing into the future light cone. The negative root states (states of Dirac's "negative energy sea") have energies at infinity given by
and have 4 -momentum vectors pointing into the past light cone. In the gap between the "seas" no orbits exist (forbidden region). The gap vanishes at the horizon r=Mr=M (infinite redshift of local energy gap, 2mu2 \mu, gives zero gap in energy-at-infinity). [Figure adapted from Christodoulou (1971).]
corresponding to changes in the black-hole properties of
{:(33.56)delta M >= (a delta S+r_(+)Q delta Q)/(r_(+)^(2)+a^(2))quad((" absolute minimum value of ")/(delta M" for given "delta S" and "delta Q)):}\begin{equation*}
\delta M \geq \frac{a \delta S+r_{+} Q \delta Q}{r_{+}{ }^{2}+a^{2}} \quad\binom{\text { absolute minimum value of }}{\delta M \text { for given } \delta S \text { and } \delta Q} \tag{33.56}
\end{equation*}
Notice an important consequence [Bardeen (1970a)]: if the black hole is initially of the "extreme Kerr-Newman" variety, with M^(2)=a^(2)+Q^(2)M^{2}=a^{2}+Q^{2}, so that one might fear a change which makes M^(2) < a^(2)+Q^(2)M^{2}<a^{2}+Q^{2} and thereby destroys the horizon, one's fears are unfounded. Equation (33.56) then demands (since r_(+)=Mr_{+}=M and S=MaS=M a )
M delta M >= a delta a+Q delta QM \delta M \geq a \delta a+Q \delta Q
so M^(2)M^{2} remains greater than or equal to a^(2)+Q^(2)a^{2}+Q^{2}, and the horizon is preserved.
Box 33.5 ORBITS OF TEST PARTICLE IN "EQUATORIAL PLANE" OF KERR-NEWMAN BLACK HOLE
Radial motion is governed by energy equation (33.54) with theta=p^(theta)=0\theta=p^{\theta}=0 :
{:(3)p^(r)=dr//d lambda:}\begin{equation*}
p^{r}=d r / d \lambda \tag{3}
\end{equation*}
Thus, equation (1) is an ordinary differential equation for dr//d lambdad r / d \lambda.
Qualitative features of the radial motion can be read off an effective-potential diagram. The effective potential V(r)V(r) is the minimum allowed value of EE at radius rr :
As in the Schwarzschild case (Figure 25.2), the allowed regions for a particle of energy-at-infinity EE are the regions with V(r) <= EV(r) \leq E; and the turning points ( p^(r)=p^{r}=dr//d lambda=0d r / d \lambda=0 ) occur where V(r)=EV(r)=E.
Stable circular orbits occur at the minima of V(r)V(r). By examining V(r)V(r) closely, one finds that for uncharged black holes the innermost stable circular orbit (most tightly bound orbit) has the characteristics here tabulated [table adapted from Ruffini and Wheeler (1971b)].
For a charged extreme Kerr-Newman black hole ( M^(2)=Q^(2)+a^(2),Q!=0M^{2}=Q^{2}+a^{2}, Q \neq 0 and a!=0a \neq 0 ) stable circular orbits with 100 per cent binding (E=0)(E=0) are achieved in the limit
(e)/( mu)longrightarrow-oo,quad(Q)/(M)longrightarrow0" (so "a longrightarrow M" ), "quad" and "((e)/( mu))*((Q)/(M))longrightarrow-oo.\frac{e}{\mu} \longrightarrow-\infty, \quad \frac{Q}{M} \longrightarrow 0 \text { (so } a \longrightarrow M \text { ), } \quad \text { and }\left(\frac{e}{\mu}\right) \cdot\left(\frac{Q}{M}\right) \longrightarrow-\infty .
[Christodoulou and Ruffini (1971)].
The effective potential for an uncharged, extreme Kerr black hole (a=M)(a=M) is shown in the figure [figure adapted from Ruffini and Wheeler (1971b)]. For detailed diagrams of orbits in the equatorial plane, see de Felice (1968). For many interesting properties of orbits that are not confined to the equatorial plane, see Wilkins (1972).
The general limit (33.56) on the change in mass can be rewritten in an alternative form [Christodoulou (1970), Christodoulou and Ruffini (1971)]:
is the "irreducible mass" of the black hole. Equation (33.57) states that no black-hole transformation produced by the injection of small lumps of matter can ever reduce the irreducible mass of a black hole. This result is actually a special case of the second law of black-hole dynamics, since the surface area of a black hole is
(Exercise 33.12).
Equation (33.58) can be combined with a=S//Ma=S / M and inverted to yield
{:[(33.60)M^(2)=(M_(ir)+(Q^(2))/(4M_(ir)))^(2)+(S^(2))/(4M_(ir)^(2)).],[[[" irreducible con- "],[" tribution to mass "]]{:[(uarr)/(" electromagnetic con- "){:" tribution to mass ":}],[[[" rotational con- "],[" tribution to mass "]]]:}]:}\begin{align*}
& M^{2}=\left(M_{\mathrm{ir}}+\frac{Q^{2}}{4 M_{\mathrm{ir}}}\right)^{2}+\frac{S^{2}}{4 M_{\mathrm{ir}}{ }^{2}} . \tag{33.60}\\
& {\left[\begin{array}{l}
\text { irreducible con- } \\
\text { tribution to mass }
\end{array}\right] \begin{array}{l}
\frac{\uparrow}{\text { electromagnetic con- }} \begin{array}{l}
\text { tribution to mass }
\end{array} \\
{\left[\begin{array}{l}
\text { rotational con- } \\
\text { tribution to mass }
\end{array}\right]}
\end{array}}
\end{align*}
A black-hole transformation that holds fixed the irreducible mass is reversible; one that increases it is irreversible. The derivation of equation (33.56) revealed that the only injection processes that actually achieve the minimum possible value for delta M\delta M (and thus make deltaM_(ir)=0\delta M_{\mathrm{ir}}=0 ) are those with mu=p^(r)=p^(theta)=0\mu=p^{r}=p^{\theta}=0 at the horizon, r=r_(+)r=r_{+}. Restated in words: To produce a reversible transformation by injecting an object into a black hole, one must (1) give the object a rest mass mu\mu extremely small compared to its charge epsilon\epsilon or axial component of angular momentum L_(z)L_{z},
mu//e≪1" and "//" or "mu^(2)//L_(z)^(2)≪1\mu / e \ll 1 \text { and } / \text { or } \mu^{2} / L_{z}^{2} \ll 1
and (2) set the object down "extremely gently" ( p^(r)=p^(theta)=0p^{r}=p^{\theta}=0 ), extremely close to the horizon (r=r_(+))\left(r=r_{+}\right). This does not sound too difficult until one recalls that objects with p^(r)=p^(theta)=0p^{r}=p^{\theta}=0 at the horizon must be moving outward with the speed of light, and that the nearer one approaches the horizon as one sets down the object, the greater one's danger of "slipping" and getting swallowed!
Clearly, any actual injection process will depart somewhat from irreversibility. Reversibility is an idealized limit, approachable but not attainable.
EXERCISES Exercise 33.11. IRREDUCIBLE MASS IS IRREDUCIBLE
Show that condition (33.56) is equivalent to deltaM_(ir) >= 0\delta M_{\mathrm{ir}} \geq 0.
Exercise 33.12. SURFACE AREA OF A BLACK HOLE
Show that the surface area of the horizon of the Kerr-Newman geometry [area of surface r=r_(+)r=r_{+}and t=t= const (Boyer-Lindquist coordinates) or widetilde(V)=\widetilde{V}= const (Kerr coordinates) ]] is 16 piM_(ir)^(2)16 \pi M_{\mathrm{ir}}{ }^{2}.
Exercise 33.13. ANGULAR VELOCITY OF A BLACK HOLE
A general theorem [Hartle (1970) for relativistic case; Ostriker and Gunn (1969) for nonrelativistic case] says that, if one injects angular momentum into a rotating star while holding fixed all other contributions to its total mass-energy (contributions from entropy and from baryonic rest mass), then the injection produces a change in total mass-energy given by
{:(33.61)delta(" mass-energy ")=((" angular velocity of star ")/(" at point of injection "))delta(" angular momentum "):}\begin{equation*}
\delta(\text { mass-energy })=\binom{\text { angular velocity of star }}{\text { at point of injection }} \delta(\text { angular momentum }) \tag{33.61}
\end{equation*}
By analogy, if one injects an angular momentum delta S\delta S into a rotating black hole while holding fixed all other contributions to its total mass-energy (contributions from irreducible mass and from charge), one identifies the coefficient Omega_(h)\Omega_{h} in the equation
delta M=Omega_(h)delta S\delta M=\Omega_{h} \delta S
Notice that this is precisely the angular velocity of photons that live forever on the horizon [equation (33.42b); "barber-pole twist" of null generators of horizon].
(b) Show that any object falling into a black hole acquires an angular velocity (relative to Boyer-Lindquist coordinates) of Omega=d phi//dt=Omega_(h)\Omega=d \phi / d t=\Omega_{h} in the late stages, as it approaches the horizon. (Recall that the horizon is a singularity of the Boyer-Lindquist coordinates. This is the reason that every object, regardless of its L_(z),E,e,mu,QL_{z}, E, e, \mu, \mathscr{Q}, can approach and does approach Omega=Omega_(h)\Omega=\Omega_{h}.)
Exercise 33.14. SEPARATION OF VARIABLES FOR WAVE EQUATIONS
This chapter has studied extensively the motion of small objects in the external fields of black holes. Of almost equal importance, but not so well-understood yet because of its complexity, is the evolution of weak electromagnetic and gravitational perturbations ("waves") in the Kerr-Newman geometry. Just as one had no à priori reason to expect a "fourth constant" for test-particle motion in the Kerr-Newman geometry, so one had no reason to expect separability for Maxwell's equations, or for the wave equations describing gravitational perturbations-or even for the scalar wave equation ◻psi-=-psi_(alpha)^(alpha)=0\square \psi \equiv-\psi{ }_{\alpha}{ }^{\alpha}=0. Thus it came as a great surprise when Carter (1968c) proved separability for the scalar wave equation, and later when Teukolsky (1972,1973)(1972,1973) separated both Maxwell's equations and the wave equations for gravitational perturbations.
Show that separation of variables for the scalar-wave equation in the (uncharged) Kerr geometry yields solutions of the form
where mm and ℓ\ell are integers with 0 <= |m| <= ℓ;S_(mℓ)0 \leq|m| \leq \ell ; S_{m \ell} is a spheroidal harmonic [see Meixner and Schärfke (1954)]; and u_(ℓm)u_{\ell m} satisfies the differential equation
{:(33.64b)-d^(2)u//dr^(**2)+Vu=0:}\begin{equation*}
-d^{2} u / d r^{* 2}+V u=0 \tag{33.64b}
\end{equation*}
In order to put the equation in this form, define a Regge-Wheeler (1957) "tortoise"-type radial coordinate r^(**)r^{*} by
{:(33.64c)dr^(**)=Delta^(-1)(r^(2)+a^(2))dr:}\begin{equation*}
d r^{*}=\Delta^{-1}\left(r^{2}+a^{2}\right) d r \tag{33.64c}
\end{equation*}
and find an effective potential V(r^(**))V\left(r^{*}\right) given by
In this radial equation Q\mathscr{Q} is a constant (analog of Carter's constant for particle motion), given in terms of mm and ℓ\ell by
{:(33.64e)Q-=lambda_(mℓ)-m^(2);lambda_(mℓ)=[[" eigenfunction of spheroidal harmonic; "],[" see Meixner and Schärfke (1954) "]].:}\mathcal{Q} \equiv \lambda_{m \ell}-m^{2} ; \lambda_{m \ell}=\left[\begin{array}{l}
\text { eigenfunction of spheroidal harmonic; } \tag{33.64e}\\
\text { see Meixner and Schärfke (1954) }
\end{array}\right] .ä
[These details of the separated solution were derived by Brill et al. (1972). For studies of the interaction between fields and Kerr black holes-studies performed using the above solution, and using analogous solutions to the electromagnetic and gravitational wave equations-see Bardeen, Press, and Teukolsky (1972), Misner (1972b), Teukolsky (1972), Ipser (1971), Press and Teukolsky (1973), and Chrzanowski and Misner (1973).]
*This box is based on Misner (1969a).
*The global structure of the Schwarzschild geometry, including the existence of two singularities and two exterior regions, was first discovered by Synge (1950). See Box 31.1.
*On the problem of Plateau see, e.g., Courant (1937), Darboux (1941), or p. 157 of Lipman Bers (1952). †\dagger The uncharged ( Q=0Q=0 ) version was first found as a solution to Einstein's vacuum field equations by Kerr (1963). The charged generalization was first found as a solution to the Einstein-Maxwell field equations by Newman, Couch, Chinnapared, Exton, Prakash, and Torrence (1965). Only later was the connection to black holes discovered; see Box 33.1.
*This is not quite true. Kerr-Newman spacetime possesses, in addition to its two Killing vectors, also a "Killing tensor" which is closely linked to the Boyer-Lindquist coordinates rr and theta\theta. See Walker and Penrose (1970); also §33.5.
Extraction of energy from a
black hole by processes in
Extraction of energy from a
black hole by processes in the ergosphere